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Abstract. This paper is devoted to the study of propagation of chaos and mean-field 
limit for systems of indistinguable particles undergoing collision processes, as formulated 
by M. Kac [42] for a simplified model and extended by H. P. McKean [53] to the Boltz- 
mann equation. We prove quantitative and uniform in time estimates measuring the 
distance between the many-particle system and the limit system. These estimates imply 
<«_^ ■ in particular the propagation of chaos for marginals in weak measure distances but are 

Cn ' more general: they hold for non-chaotic initial data and control the complete many- 

particle distribution. We also prove the propagation of entropic chaos, as defined in [12j . 
answering a question of Kac about the microscopic derivation of the Jif-theorem. We 
finally prove estimates of relaxation to equilibrium (in Wasserstein distance and relative 
>D , entropy) independent of the number of particles. Our results cover the two main Boltz- 

mann physical collision processes with unbounded collision rates: hard spheres and true 

Maxwell molecules interactions. Starting from an inspirative paper of A. Griinbaum |36] 

Q, ■ we develop a new method which reduces the question of propagation of chaos to the one of 

.^fH ' proving a purely functional estimate on some generator operators {consistency estimates) 

. , together with fine differentiability estimates on the flow of the limit non-linear equation 

r^^ • {stability estimates). These results provide the first answer to the question raised by 

jrt \ Kac of relating the long-time behavior of a coUisional particle system with the one of its 

mean-field limit, however using dissipativity at the level of the mean-field limit instead 
of using it at the level of the many-particle Markov process. 
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L. Boltzmann is best known for the equation bearing his name in kinetic theory [71[8]. 
Inspired by Maxwell's discovery [51] of (what is now called) the Boltzmann equation 
and its "maxwellian" (gaussian) equilibrium, Boltzmann [7] discovered the "ff-theorem" 
(increase of the entropy) for this equation which explains how the solutions should be 
driven towards the equilibrium of Maxwell. In the same work he also proposed the deep 
idea of "stosszahlansatz" (molecular chaos) for explaining how the irreversible Boltzmann 
equation could emerge from the Newton laws on the dynamics of particles. Giving a precise 
mathematical meaning to this notion and proving this limit remains a tremendous open 
problem to this date (the best and astonishing result so far [l5] is only valid for very short 
times). 

M. Kac proposed in 1956 [52] the simpler and seemingly more tractable question of de- 
riving the spatially homogeneous Boltzmann equation from a many-particle jump process, 
and he introduced a rigorous notion of molecular chaoqj in this context. The "chaotic- 
ity" of the many-particle equilibrium towards the maxwellian distribution, i.e. the fact 
that the first marginals of the uniform measure on the sphere S {\/N) converges to a 
gaussian as A^ goes to infinity, has been known for a long time, at least since Maxwello 
However in [12] Kac proposed the first proof of the propagation of chaos along time for a 
simplified collision process for which series expansions of the solution are available, and he 
showed how the many-particle limit rigorously follows from this property of propagation 
of chaos. This proof was later extended to a more realistic collision model, the so-called 
cutoff Maxwell molecules, by McKean [53]. 

Since in this setting both the many-particle system and the limit equation are dissipa- 
tive, Kac also raised the natural question of relating their asymptotic behaviors. In his 
mind this program was to be achieved by understanding dissipativity at the level of the lin- 
ear many-particle jump process and he insisted on the importance of estimating the rate of 
relaxation of the many-particle process. This has motivated beautiful works on this "Kac 
spectral gap problem" [iQl [501 HSj [l5l [H] , but so far this strategy has proved unsuccessful 
in relating the many-particle process asymptotic behavior and that of the limit equation, 
see the interesting discussion in [12]. At the time of Kac the study of nonlinear partial 
differential equations was rather young and it was plausible that the study of a linear 
many-dimension Markov process would be easier. However the mathematical developpe- 
ment somehow followed the inverse direction and the theory of existence, uniqueness and 
relaxation to equilibrium for the spatially homogeneous Boltzmann is now well-developed 
(see the many references along this paper). 

This paper is an attempt to develop a quantitative theory of mean-field limit which 
strongly relies on detailed knowledge of the limit nonlinear equation, rather than on detailed 
properties of the many-particle Markov process. As the main outcome of this theory we 
prove uniform in time quantitative propagation of chaos as well as propagation of entropic 
chaos, and we prove relaxation rates independent of the number of particles (measured in 
Wasserstein distance and relative entropy). All this is done for the two important realistic 
and achetypal models of collision, namely hard spheres and true (without cutoff) Maxwell 
molecules. This provides a first complete answer to the question raised by Kac, however 
our answer is an "inverse" answer in the sense that our methodology is "top-down" from 
the limit equation to the many-particle system rather than "bottom-up" as was expecting 
Kac. 



Kac in fact called this notion "Boltzniann's property" in [15] as a clear tribute to the fundamental 
intuition of Boltzmann. 

We refer to [53] for a bibliographic discussion, see also [T5] where [M] is quoted as the first paper 
proving this result. 
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A first version of this work was posted on arXiv as "Quantitative uniform in time chaos 
propagation for Boltzmann colhsion processes" . This new expanded version includes new 
results about entropic chaos and relaxation times. 
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1. Introduction and main results 

1.1. The Boltzmann equation. The Boltzmann equation (Cf. [I9] and [20]) describes 
the behavior of a dilute gas when the only interactions taken into account are binary 
collisions. It writes 

^ + vVJ = Qif,f) 

where Q = Q{f, f) is the bilinear Boltzmann collision operator acting only on the velocity 
variable. 

In the case when the distribution function is assumed to be independent on the position 
X, we obtain the so-called spatially homogeneous Boltzmann equation, which reads 

(1.1) ^(t,v) = Q{fJ){t,v), veR'', t>0, 

where d > 2 is the dimension. 

Let us now focus on the collision operator Q. It is defined by the bilinear symmetrized 
form 

(1.2) Q{gJ)(v) = l f B{\v - v.\,cose) {g'J' + g'f:-9*f- 9 f*)dv, da, 

where we have used the shorthands / = f{v), f = f{v'), g^ = g{v^) and g'^ = g{v'^)- 
Moreover, v' and v^ are parametrized by 

(1-3) v=^— + ^ a, ^* = ^ 2 ^' '^ ■ 

Finally, 6 G [0, vr] is the deviation angle between v' — v'^ and v — v* defined by cos 9 = a -u, 
u = V — v^, u = u/\u\, and B is the Boltzmann collision kernel determined by physics 
(related to the cross-section T,{v — v^:,a) by the formula B = \v — v*\ S). 

Boltzmann's collision operator has the fundamental properties of conserving mass, mo- 
mentum and energy 

' Q{f,f)Hv)dv = 0, (^{v) = l,v,\v\\ 



and satisfying the so-called Boltzmann's H theorem which writes (at the formal level) 

_d_ 
~di 
We shall consider collision kernels 

B = r{\v-v^\)b{cos6) 



[ flogfdv = - [ Qif, f) log(/) dv > 0. 
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with r, b nonnegative functions. In dimension d = 3, here is a short classification of 
possible collision kernels, together with some important examples that we shall consider 
in this paper. 

(1) Short (finite) range interaction are usually modeled by the hard spheres collision 
kernel 

(1.4) (HS) B{\v-v^\,cose) = cst\v-v*\. 

(2) Long-range interactions are usually modeled by collision kernels deriving from in- 
teraction potentials 

V{r) = cstr""*, s > 2. 

They satisfy the formula 

r(z) = \z\^ with 7 = (s - 4)/s 

and 

b{cose) ~e~o CbO"^^'' with 1/ = 2/s 

(b is L^ apart from ^ ~ 0). Such formulas (and others) can be found in J19j . 

This general class of collision kernels includes in particular the true Maxwell 
molecules collision kernel when 7 = and u = 1/2: 

(1.5) (tMM) B{\v -v^\,cos9) = B{cose) ^er^oCbO-^/^. 

It also includes the so-called Grad's cutoff Maxwell molecules when the 
singularity in the 9 variable is removed. Without much restriction we shall consider 
as a typical such model 

(1.6) (GMM) B{\v-v^\,cose) = l. 

1.2. Deriving the Boltzmann equation from many-particle systems. The question 
of deriving the Boltzmann equation from particle systems (interacting via Newton's laws) 
is a famous problem. It is related to the so-called 6-th Hilbert problem proposed by Hilbert 
at the International Congress of Mathematics at Paris in 1900: axiomatize mechanics by 
^^ developing mathematically the lim,iting processes [■ ■ ■] which lead from, the atomistic view 
to the laws of motion of continua" . 

At least at the formal level, the correct limiting procedure has been identified by 
Grad [HI] in the late fourties and a clear mathematical formulation of the open prob- 
lem was proposed in |18j in the early seventies. It is now called the Boltzmann- Grad or 
low density limit. However the original question of Hilbert remains largely open, in spite 
of a striking breakthrough due to Lanford [15], who proved the limit for short times (see 
also Illner and Pulvirenti |39j for a close-to- vacuum result). The tremendous difficulty 
underlying this limit is the irreversibility of the Boltzmann equation, whereas the particle 
system interacting via Newton's laws is a reversible Hamiltonian system. 

In 1954-1955, Kac [52] proposed a simpler and more tractable problem: start from the 
Markov process corresponding to collisions only, and try to prove the limit towards the 
spatially homogeneous Boltzmann equation. Kac's jump process runs as follows: consider 
N particles with velocities fi,. . . , V]\[ S W^. Compute random times for each pair of 
particles (fj, Vj) following an exponential law with parameter T(\vi — Vj\), take the smallest, 
and perform a collision {vi,Vj) — >• {v* ,v*) given by a random choice of a direction parameter 
whose rule is related to b{cos0), then recommence. This process can be considered on 
M ; however it leaves invariant some submanifolds of M (depending on the number of 
conserved quantities during collision) and can be restricted to them. For instance in the 
original simplified model of Kac d = 1 (scalar velocities) and the process can be restricted 
to S^^^{\^£N) the sphere with radius ySN, where £ is the energy. In the more realistic 
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hard spheres of Maxwell molecules models, d = 3 and the process can be restricted to the 
sphere 

gN ._ ^dN-1 (^^]vf ) r\{vi + ---+VN=0}. 

Kac formulated the notion of propagation of chaos that we shall now explain. Consider 
a sequence {f )n>i of probabilities on E , where E is some given Polish space (think to 
E = W^ for most applications): the sequence is said f- chaotic if 

/^~/®^ when A^^oo 

for some given one-particle probability / on E. The meaning of this convergence is the 
following: convergence in the weak measure topology for any marginal depending on a 
finite number of variables. This is a low correlation assumption. 

It was clear since Boltzmann that in the case when the joint probability density / 
of the A^-particle system is tensorized during some time interval into A^ copies /®" of a 
1-particle probability density, then the latter would satisfy the limit nonlinear Boltzmann 
equation during this time interval. But Kac made the key remark that although in general 
interactions between a finite number of particles prevents any possibility of propagation of 
the "tensorization" property, the weaker property of chaoticity can be propagated (hope- 
fully!) in the correct scaling limit. 

The framework set by Kac is our starting point. Let us emphasize that the limit 
performed in this setting is different from the Boltzmann-Grad limit. The Kac's limit is 
in fact a mean-field limit. This limiting procedure is most well-known for deriving Vlasov- 
like equations. In a companion paper [58j we develop systematically our new functional 
approach for Vlasov equations, McKean-Vlasov equations, and granular gases driven by a 
thermal bath. 

1.3. The notion of chaos and how to measure it. Our goal in this paper is to set 
up a general robust method for proving the propagation of chaos with quantitative rate in 
terms of the number of particles N and of the final time of observation T. 

Let us explain briefly what it means. The original formulation of Kac [l2] of chaoticity 



is: a sequence / G Psym{E ) of symmetric probabilities on E is /-chaotic, for a given 
probability / G PiE), if for any i£W and any ip G Cb{E)®^ there holds 

lim (/^^®l^-A = (r^v^ 

N^oo \ / \ 

which amounts to the weak convergence of any marginals. This can be expressed for 
instance in Wasserstein distance: 

hm W, (Ue in , f^') = 

where II^ denotes the marginal on the i first variables. This is Kac's definition of chaos 
that we shall call finite-dimensional chaos. 

We shall deal in this paper with quantitative chaos, in the sense that we measure 
precisely the rate of convergence in the above limit. Namely, we say that /^ is /-chaotic 
with rate e{N), where s{N) — )• when TV — > oo (typically e{N) = N~''' , r > or e{N) = 
{logN)~'^ , r > 0), if for any £ G N* there exists K^ G (0, oo) such that 

Similar statements with other metrics can be also formulated (and shall be used in this 
paper): for some normed space of smooth functions T C Ch{E) (to be specified) and for 
any £ G N* there exists Ki G (0,oo) such that for any ip G J-"®^, \\^\\j^ < 1, there holds 

'u4f^]-f^',^)\<K,e{N). 

The Wasserstein distance Wi is recovered when J^ is the space of Lipschitz functions. 
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Observe that in the latter statements the number of variables I considered in the mar- 
ginal is kept fixed as A^ goes to infinity. A stronger notion of infinite-dimensional chaos 
would be 

lim — = 

N^oo N 

with corresponding quantitative formulations. This amounts to say that one can prove a 
sublinear control on Ki in terms of H. in the previous statements. Variants for other metrics 
could also be considered. 

Finally one can formulate an even stronger notion of (infinite-dimensional) entropic 
chaos (see [12] and definition ()1.7p of the relative entropy below): 

lif(/^|7^)^^^^i/(/|7) 

with obvious quantitative versions. This notion of chaos is particularly interesting since 
it corresponds to the derivation of Boltzmann's entropy from the many-particle system 
entropies, we shall come back to this point. 

Now, considering a sequence of symmetric^ A^-particle densities 

/^eC([0,oo);P,,„,(i?^)) 

and a 1-particle density of the expected mean field limit 

/G C([0,oo); P(i?)), 

we say that there is propagation of chaos on some time interval [0, T] if the /o-chaoticity 
of the initial family f^ implies the /j-chaoticity of the family //^ for any time t G [0,T], 
according to one of definitions of chaoticity above. 

1.4. Kac's program. As already said Kac proposed the problem of deriving the spatially 
homogeneous Boltzmann equation from a many-particle Markov jump process modeling 
the binary collisions, through its master equation (the equation on the law of the process). 
This amounts intuitively to consider the spatial variable as a hidden variable inducing 
ergodicity and markovian properties on the velocity variable. Although the latter point 
has not been proved so far to our knowledge, it is worth pointing out that it is at the 
same a very natural guess and an extremely interesting (and probably also difficult) open 
problem. 

Interestingly enough here is in the words of Kac [42, p. 173] how his approach was 
inspired by Boltzmann's ideas: "This formulation led to the well-known paradoxes which 
were fully discussed in the classical article of P. and T. Ehrenfest. These writers made it 
clear (a) that the "Stosszahlansatz" cannot be strictly derivable from purely dynamic con- 
siderations and (b) that the "Stosszahlansatz" has to be interpreted probabilisticaly. The 
recent attempts of Born and Green, Kirkwood and Bogoliubov to derive Boltzmann's equa- 
tion from Liouville's equation and hence to justify the "Stosszahlansatz" dynamically are, 
in our opinion, incomplete, inasmuch as they do not make it clear at what point statistical 
assumptions are introduced. The "master equation" approach which we have chosen seems 
to us to follow closely the intentions of Boltzmann. " Even we partly disagree with this 
statement and we still believe that a fully satisfying justification of the "Stosszahlansatz" 
(molecular chaos) has to be derived from "purely dynamic consideration" (Newton's laws) 
as was pioneered by Lanford [l5], the latter open problem seems tremendously difficult and 
it is clear that Kac raised a fascinating more tractable question: if we have to introduce 
stochasticity, at least can we keep it under control all along the process of derivation of 
the Boltzmann equation and relate it to the dissipativity of the limit equation. This is one 
of the main questions we attempt to answer in this work. 



i.e. invariant according to permutations of the particles. 
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Let US now discuss more in details the content of Kac's seminal paper |42J . He goes on 
to introduce a jump process on the (A^ — l)-dimensional sphere S (viV) with radius 
yiV corresponding to random rotations among pairs of variables chosen randomly, and 
occuring at random times following exponential laws. He then prove finite-dimensional 
chaos (with no rate) by a beautiful combinatorial argument, based on an infinite series 
"tree" representation of the solution according to the collision history of particles, and 
some Leibniz formula for the iV-particle operator acting on tensor products. 

He then raises several questions that we schematize as follows: 

(1) The first one concerns the restriction of the models as compared to realistic collision 
processes. The fact that the geometry of collision is simplified was relaxed by 
McKean [53] soon afterwards, however for cutoff Maxwell molecules for which the 
collision kernel does not depend on the modulus of the relative velocity. This 
seemingly technical issue is in fact related to deep difficulty for dealing with jump 
process whose jump times follow laws depending on the variables. In the words 
of Kac [421 p. 179]: "The above proof suffers from the defect that it works only 
if the restriction on time is independent of the initial distribution. It is therefore 
inapplicable to the physically significant case of hard spheres because in this case our 
simple estimates yield a time restriction which depends on the initial distribution. 
A general proof that Boltzmann's property propagates in time is still lacking." In 
other words the question is: can one prove propagation of chaos for the 
hard spheres coUision process? 

(2) Following closely the spirit of the previous question about going beyond the limit 
of Kac's original combinatorial insight in order to deal with realistic collision pro- 
cesses, it seems to us very natural to ask whether one can prove propagation of 
chaos for the true Maxwell molecules collision process? The difficulty lies 
now in the fact that the particle system can undergo infinite number of collisions 
in a finite time interval, and no "tree" representation of solutions is available. This 
is related with the physical interesting situation of long-range interactions, as well 
as with the mathematical interesting framework of fractional derivative operator 
and Levy walk. 

(3) Kac then discusses the //-theorem of Boltzmann, which is not surprising as its 
original goal is the derivation of Boltzmann's equation and its irreversible feature. 
He makes the simple observation [421 P- 182] that the "ergodic property of the 
Markov process under consideration" steadily implies that it admits an infinite 
number of Liapunov functions, including the L^ norm and Boltzmann's entropy. 
In contrast with it, the limit equation admits only (in general) Boltzmann's en- 
tropy as a Liapunov function. Kac then heuristically conjectures [121 Eq. (6.39)] 
that (in our notations) H{f^)/N — >• H{ft) along time, which would imply Boltz- 
mann's H-theorem. from the monotonicity of H{f )/N for the Markov process. 
He concludes with: "/f the above steps could be made rigorous we would have a 
thoroughly satisfactory justification of Boltzmann's H -theorem." In our notation 
the question is can one prove propagation of entropic chaos along time? 

(4) He finally discusses the relaxation times, with the goal of deriving relaxation times 
of the limit equation from the many-particle system. This imposes to have esti- 
mates independent of the number of particles on this relaxation times. As a first 
natural step he therefore goes on to consider the LP' spectral gap of the Markov 
process on the sphere and remark [111 p. 187]: "Surprisingly enough this seems 
quite difficult and we have not succeeded in finding a proof. Even for the simplified 
model we have been considering, the question remains unsettled although we are 
able to give a reasonably explicit solution of the master equation." This question 
has triggered many beautiful works (see the next subsection), however it is easy 
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to convince oneself (see the discussion in [12j for instance) that there is no hope 
of passing to the hmit A^ — )■ oo in this spectral gap estimate, even if the spec- 
tral gap is independent of A^. The LP' norm is catastrophic in infinite dimension. 
Therefore following quite closely the intention of Kac, we reframe the question in 
a setting which "tensorizes correctly in the limit N — )• oo" , that is in our notation: 
can one prove relaxation times independent of the number of particles on the 
normalized Wasserstein distance W{f^,j^)/N or on the normalized rel- 
ative entropy H{f^\j^)/N, w^here 7^ denotes the A^-particle invariant 
measure? 

This paper is concerned with solving the four questions outlined above. 

Before concluding this subsection, let us briefly illustrate by some quotation \i'2\ p. 178] 
that Kac was firmly believing in a "bottom-up" approach, which deduces properties on 
the Boltzmann equation from the Markov process: ^^Since the master equation is truly 
descriptive of the physical situation, and since existence and uniqueness of the solution 
of the master equation are almost trivial, the preoccupation with existence and uniqueness 
theorems for the Boltzmann equation appears to be unjustified on grounds of physical in- 
terest and importance." Although as the time Kac was writing, almost no mathematical 
result (apart from [9l [10]) was available for the nonlinear spatially homogeneous Boltz- 
mann equation, mathematical development has not followed this path, and as we have said, 
our approach shall be rather "top-down" , taking advantage of the bynow well-developped 
theory for this equation. 

1.5. Review of previous results. For Boltzmann collision processes, Kac |42j-|43) has 

proved the point (1) in the case of his baby one-dimensional model. The key point in 
his analysis is a clever combinatorial use of a semi-explicit form of the solution (Wild 
sums). It was generalized by McKean [53j to the Boltzmann collision operator but only 
for "Maxwell molecules with cutoff", i.e. roughly when the collision kernel B above is 
constant. In this case the combinatorial argument of Kac can be extended. Kac raised 
in [42 1 the question of proving propagation of chaos in the case of hard spheres and more 
generally unbounded collision kernels, although his method seemed impossible to extend 
(no semi-explicit combinatorial formula of the solution exists in these cases) . 

In the seventies, Griinbaum |36) then proposed in a very compact and abstract paper 
another method for dealing with hard spheres, based on the Trotter-Kato formula for 
semigroups and a clever functional framework (partially remindful of the tools used for 
mean-field limit for McKean- Vlasov equations). Unfortunately this paper was incomplete 
for two reasons: (1) It was based on two "unproved assumptions on the Boltzmann flow" 
(page 328): (a) existence and uniqueness for measure solutions and (b) a smoothness 
assumption. Assumption (a) was indeed recently proved in [32] using Wasserstein metrics 
techniques and in [29] adapting the classical DiBlasio trick [22] , but concerning assumption 
(b), although it was inspired by cutoff maxwell molecules (for which it is true), it fails for 
hard spheres (cf. the counterexample built by Lu and Wennberg in |49] ) and is somehow 
"too rough" in this case. (2) A key part in the proof in this paper is the expansion of the 
"-ffj" function, which is a clever idea of Griinbaum (and the starting point for our idea of 
developing an abstract differential calculus in order to control fluctuations) — however it 
is again too rough and is adapted for cutoff Maxwell molecules but not for hard spheres. 

A completely different approach was undertaken by Sznitman in the eighties [5^ (see 
also Tanaka [68] for partial results concerning non-cutoff Maxwell molecules). Starting 
from the observation that Griinbaum's proof was incomplete, he gave a full proof of prop- 
agation of chaos for hard spheres. His work was based on: (1) a new uniqueness result 
for measures for the hard spheres Boltzmann equation (based on a probabilistic reasoning 
on an enlarged space of "trajectories"); (2) an idea already present in Griinbaum's ap- 
proach: reduce by a combinatorial argument on symmetric probabilities the question of 
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propagation of chaos to a law of large numbers on measures; (3) a new compactness result 
at the level of the empirical measures; (4) the identification of the limit by an "abstract 
test function" construction showing that the (infinite particle) system has trajectories in- 
cluded in the chaotic ones. Hence the method of Sznitman proves convergence but does 
not provide any rate for chaoticity. 

Let us also emphasize that McKean in [52] studied fluctuations around deterministic 
limit for 2-speed Maxwellian gas and for the usual hard balls gas. Graham and Meleard 
in [3_5j have obtained a rate of convergence (of order 1/A^ for the ^-th marginal) on any 
bounded finite interval of the A^-particle system to the deterministic Boltzmann dynamic 
in the case of Maxwell molecules under Grad's cut-off hypothesis. Finally Fournier and 
Meleard in |30tl31j have obtained the convergence of the Monte-Carlo approximation (with 
numerical cutoff) of the Boltzmann equation for true Maxwell molecules with a rate of 
convergence (depending on the numerical cutoff and on the number N of particles) . 

After we had finished writing our paper, we were told by I. Bailleul about the recent book 
|44j by V. N. Kolokoltsov. This book is focused on fluctuation estimates of central limit 
theorem type. It does not prove quantitative propagation of chaos but weaker estimates 
(and on finite time intervals), however the comparison of generators for the many-particle 
and the limit semigroup is reminiscent of our work. 

1.6. The method. The main inspiration at first came from Griinbaum's paper [36]. Our 
original goal was to construct a general and robust method able to deal with mixture 
of jump and diffusion processes, as it occurs for granular gases (see for this point the 
companion paper [SE]). It turns out that it leads us to develop a new theory, inspiring 
from more recent tools such as the course of Lions on "Mean-field games" at College 
de France, and the master courses of Meleard [55j and Villani [70] on mean-field limits. 
One of the byproduct of our paper is that we make fully rigorous the original intuition 
of Griinbaum in order to prove propagation of chaos for the Boltzmann velocities jump 
process associated to hard spheres contact interactions. 

As Griinbaum [36j we shall use a duality argument. We introduce S^ the semigroup 
associated to the flow of the A'^-particle system and T^ its "dual" semigroup. We also 
introduce S^^ the (nonlinear) semigroup associated to the meanfield dynamic (the expo- 
nent "NL" recalling that the limit semigroup is nonlinear in the most physics interesting 
cases) as well as T^ the associated (linear) "pushforward" semigroup (see below for the 
definition). Then we will prove the above kind of convergence on the linear semigroups 
Tf and Tt°°. 

The first step consists in defining a common functional framework in which the A^- 
particle dynamic and the limit dynamic make sense so that we can compare them. Hence 
we work at the level of the "full" limit space P{P(E)) (see below). Then we shall identify 
the regularity required in order to prove the "consistency estimate" between the genera- 
tors G^ and G°° of the dual semigroups T^ and T^ , and finally prove a corresponding 
"stability estimate" at the level of the limiting semigroup S^^ . 

The latter crucial step shall lead us to introduce an abstract differential calculus for 
functions acting on measures endowed with various metrics. More precisely, we shall 
define functions of class C^' on a probability space by working on affine subspaces of the 
probability space for which the tangent space has a Banach space structure. This notion is 
related but different from the notion of derivability developed in the theory of gradient flow 
by Ambrosio, Otto, Villani and co-authors in [21 ffll [61] as well as to the differentiability 
notion introduced by Lions in |46j . 

Another viewpoint on this method is to consider it as some kind of accurate version (in 
the sense that it establishes a rate of convergence) of the BBGKY hierarchy method for 
proving propagation of chaos and mean- field limit on statistical solutions. This viewpoint is 
extensively explored and made rigorous in Section [8] where we revisit the BBGKY method 
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for the spatially homogeneous equation as developped in [3j. The proof of uniqueness for 
statistical solutions to the hierarchy becomes straightforward within our framework by 
using differentiability of the limit semigroup as a function acting on probabilities. 

This general method is the purpose of Theorem 13.11 It is, we hope, interesting per se for 
several reasons: (1) it is fully quantitative, (2) it is highly flexible in terms of the functional 
spaces used in the proof, (3) it requires a minimal amount of informations on the A^-particle 
systems but more stability information on the limiting PDE (we intentionally presented 
the assumption as for the proof of the convergence of a numerical scheme, which was our 
"methodological model"), (4) the "differential stability" conditions that are required on 
the limiting PDE seem (to our knowledge) new, at least at the level of Boltzmann or more 
generally transport equations. 

1.7. Main results. Let us give some simplified versions of the main results in this paper. 
All the abstract objects shall be fully introduced in the next sections. 

Theorem 1.1 (Summary of the main results). Consider some initial distribution /o G 
P{R'^)nL°°{W^) with compact support or polynomial moment bounds. Consider the corre- 
sponding solution ft to the spatially homogeneous Boltzmann equation for hard spheres or 
Maxwell molecules, and the solution f^ of the corresponding Kac 's jump process starting 
either from the N-fold tensorization of fo or the latter conditionned to S . 
The results can be classified into three main statements: 

(1) Quantitative uniform in time propagation of chaos, finite or infinite 
dimensional, in weak measure distance (cf. Theorems \5.1^5.3A6.R6.2^) : 

w, (n,/f , (ft 

ViV> 1, Vl<£<iV, sup ^ -^ '-^<a{N) 

for some a{N) —)• as A^ —)• oo. In the case f^ = /® one has moreover explicit 
power law (for Maxwell molecules) or logarithmic rate (for hard spheres) estimates 
on a, 

(2) Propagation of entropic chaos (cf. Theorem 7.10\ -(i)): Consider the case where 



the initial datum of the many-particle system is restricted to S^ . If the initial 
datum is entropy- chaotic in the sense 



^H{f,^\j^)^^^±^Hifoh) 



N 



with 



1.7) F(/o^|7^):= / log %f^{dV) and H (foh) := [ folog^dv 

and where 7 is the gaussian equilibrium with energy £ and 7^ is the uniform 
probability measure on S^ , then the solution is also entropy chaotic for any later 
time: 

Vt>0, ^^(/f|7^)^^^^±^^(/t|7). 

This proves the derivation of the H -theorem in this context. 
(3) Quantitative estimates on relaxation times, independent of the number 
of particles (cf. Theorems \5.3\^6^ and Theorem \7.10\ - (ii)) : Consider the case 
where the initial datum of the many-particle system is restricted to S . Then we 
have 

ViV> 1, VI <£<iV, Vi>0, ^ ^ ^' '' </?(t) 
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iu 



for some (3{t) — t- as t —)■ oo. Moreover in the case of Maxwell molecules, and 
assuming moreover that the Fisher information of the initial datum /o is finite: 

Hfo) ■■= I — 7 dv < +00, 



/o 



the following estimate also holds: 

1 



ViV>l, -if(/f|7^)</3(t) 



for some function I3{t) -^ as t 



oo. 



l.i 



Some open questions. Here are a few questions among those raised by this work: 

(1) What about the optimal rate in the chaoticity estimates along time? Our method 
reduces this question to the chaoticity estimates at initial time, and therefore to 
the optimal rate in the quantitative law of large numbers for measures according 
to various weak measure distances. 

(2) What about the optimal rate in the relaxation times (uniformly in the number 
of particles)? Spectral gap studies predict exponential rates, both for the many- 
particle system and for the limit system, however our rates are far from it! 

(3) Can uniform in time propagation of chaos be proved for non-reversible jump pro- 
cesses (such as inelastic collision processes) for which the invariants measures 7^ 
and 7 are not explicitely known (e.g. granular gases)? 

1.9. Plan of the paper. In Section [2] we set the abstract functional framework together 
with the general assumption and in Section[3]we state and prove the abstract Theorem l3.1l 
In SectionUlwe present some tools and results on weak measure distances, the construction 
of initial data restricted to submanifolds for the many-particle system, and the rate for 
sampling the limit distribution by empirical measures. In Section we apply the method 
to (true) Maxwell molecules (Theorems 15.11 and 15. 3|) . In Section we apply the method 
to hard spheres (Theorems 16.11 and 16. 2p . Section [7] is devoted to the study of entropic 
chaos. Finally in Section [5] we revisit the BBGKY hierarchy method for the spatially 
homogeneous Boltzmann equation in the light of our framework. 



2. The abstract setting 

In this section we shall state and prove the key abstract result. This will motivate the 
introduction of a general functional framework. 

2.1. The general functional framework of the duality approach. Let us set up the 

framework. Here is a diagram which sums up the duality approach (norms and duality 
brackets shall be specified in Subsections [ 



N itzxN 



E'ye 



-^=/^^ 



Pn{E) C PiE) 



observables 



Kolmogorov 



Kolmogorov 



-P.,^(i?^)- 



P{PiE)) 



duality 



duality 



-^Ck{E^) 
\ 



-c 



R^ 



^C,{P{E)) 



In this diagram: 
E denotes a Polish space: 



12 S. MISCHLER AND C. MOUHOT 

This is a separable completely metrizable topological space. We shall de- 
note the distance on this space by dE in the sequel. 

- & denotes the N -permutation group. 

- Psym{E^) denotes the set of symmetric probabilities on E^: 

Given a permutation a G &^ , a vector 

V = {v,,...,VN)eE'', 

a function 99 G Ch{E ) and a probability / G P{E ), we successively 
define 

Va = K(l),---,'^a(Ar)) G E^, 

and 

ipa G Cb{E^) by setting ^^(F) = ip{V^) 

and finally 

/i^ G P{E^) by setting (/j^, ^> = (/^, 9..) • 

We then say that a probability /^ on i?^ is symmetric if it is invariant 
under permutations: 

VaG6^, /i^ = /^. 

- The probability measure fiy denotes the empirical measure: 



l'v--=^Y.^-^^ V = ivi,.. .,vn)€E^ 



where 6y. denotes the Dirac mass on E at point Vi G E. 

Pn{E) denotes the subset {fJ.y , V G E^} of empirical measures of P{E). 

P{P{E)) denotes the set of probabilities on the polish space P{E). 

Ch {P{E)) denotes the space of continuous and bounded functions on P{E): 

This space shall be endowed with either the weak or strong topologies (see 
Subsection 12. 4p . and later with some metric differential structure. 

The map tt^ from E^ /G^ to Pn{E) is defined by 

The map tt^ from Cb{P{E)) to Cb{E^) is defined by 

v$ G c,{P{E)) , vy G S^, {^^^) {V) := $ {^i^) . 

The map R^ from Cb{E^) to Cb{P{E)) is defined by: 

V(^ G Cb{E^), yp G P{E), <(p) := (p«^, vp) . 
The map vr^ from P,y^{E^) to PiPiE)) is defined by: 

for any /^ G Psym(-E^) and any <I> G Cb{P{E)), where the first bracket 
means (•, ■)p(p(E)),Cb{P(E)) ^nd the second bracket means (•, •)p{E^),C'i,(E^)- 
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Let US now discuss the "horizontal" arrows: 

- The arrows pointing from the first column to the second one consists in writing the 
Kolmogorov equation associated with the many-particle stochastic Markov process. 

- The arrows pointing from the second column to the third column consists in writing the 
dual evolution semigroup (note that the A^-particle dynamics is linear). As we shall discuss 
later the dual spaces of the spaces of probabilities on the phase space can be interpreted 
as the spaces of observables on the original systems. 

Remark 2.1. Consider a random variable V on E^ with law f^ S P{E^). It is common 
notation in probability to simply write /Xy for the random variable on P{E) with law 
TTp/^ G P{P{E)). Our notation is slightly less compact and intuitive, but at the same 
time more accurate. 

Remark 2.2. Our functional framework shall be applied to weighted probability spaces 
rather than directly in P{E). More precisely, for a given weight function m : E ^ M^- we 
shall use affine subsets of the weighted probability space 

{f€P{Ey, M„(/):=(/,m)<oo} 

as our basis functional spaces. Typical examples are 

m{v) := m{dE{v,vo)) for some fixed vq ^ E with m{z) = z or m{z) = 
e»^*, a,k>0. 

We shall sometimes abuse notation by writing M^ for Mm when m{z) = z^ in the above 
example. 

2.2. The A^-particle semigroups. Let us introduce the mathematical semigroups de- 
scribing the evolution of objects living in these spaces, for any A^ > 1. 

Step 1. Consider the trajectories V/^ G E^ , t > 0, of the particles (Markov process 
viewpoint). We make the further assumption that this flow commutes with permutations: 

For any a G (3^, the solution at time t starting from (V(^) is (V/^) . 

This reflects mathematically the fact that particles are indistinguishable. 

Step 2. This flow on E^ yields a corresponding semigroup S^^ acting on Psy„i{E^) for 
the probability density of particles in the phase space E^ (statistical viewpoint), defined 
through the formula 

V/^GP.,„.(ii;^), ^eC,{E^l (5f(/^),^>=E(^(Vf)) 

where the bracket obviously denotes the duality bracket between P{E^) and Cb{E^) and 
E denotes the expectation associated to the probability space in which the process V/^ 
is built. Since the flow (V^ ) commutes with permutation, the semigroup 5^ acts on 
Psym{E^). In other words, if the law /q^ of Vq belongs to PsymiE^), then for later times 
the law //^ of V/^ also belongs to Psym{E^). And one can associate to this semigroup a 
linear evolution equation with generator denoted by Ajy: 



dtr = AV'\ r ePsy^E' 

which can be interpreted as the forward Kolmogorov equation on the law of Vf . 

Step 3. We assume that there exists a semigroup T/^ acting on the functions space Cb{E^) 
of observables on the evolution system (V/^) on E^ (see the discussion in the next remark) 
which is dual to the semigroup 5^. More precisely, we assume 

yf'^ePiE^ ^eCb{E^), (/^^,^(^)> = (5f(n,^). 
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We associate to this semigroup the following linear evolution equation with generator 
denoted by G^: 

5i^ = G^M, ^eCb{E^). 



2.3. The mean-field limiting semigroup. We now define the evolution of the limiting 
mean-field equation. 

Step 1. Consider a semigroup S^^ acting on P{E) associated with an evolution equation 
and some operator Q: 

For any /q G P{E) (assuming possibly some additional moment bounds), 
then 5f^(/o) := ft where ft G C(M+,P(S)) is the solution to 

(2.1) dtft = Q{ft), f\t=o = fo- 

This semigroup and the operator Q are typically nonlinear for mean-field 
models, namely bilinear in case of Boltzmann's collisions interactions. 

Step 2. Consider then the associated pushforward semigroup T^ acting on Ch{P{E)): 

v/ G p{E), $ G c,{p{E)), rrmu) ■■= ^ {s^if)) ■ 

(Again additional moment bounds can be required on / in order to make this definition 
rigorous.) Note carefully that Tj°° is always linear as a function of $, although of course 
T^°° [$](/) is not linear in general as a function of /. We shall associate (when possible) 
the following linear evolution equation on Cb{P{E)) with some generator denoted by G°°: 

at$ = G°°($). 

Remark 2.3. The semigroup T^ can be interpreted physically as the semigroup of the 
evolution of observables of the nonlinear equation (|2.ip . Let us give a short heuristic 
explanation. Consider a nonlinear ordinary differential equation 

-_ = F{V) on M°' with V^, • F = and V\t=o = v 

with divergence-free vector field for simplicity. One can then define formally the linear 
Liouville transport partial differential equation 

5t/ + V,-(F/) = 0, 

where / = ft{v) is a time-dependent probability density over the phase space M"^, whose 
solution is given (at least formally) by following the characteristics backward ft{v) = 
fo{y-t{v))u Now, instead of the Liouville viewpoint, one can adopt the viewpoint of ob- 
servables, that is functions depending on the position of the system in the phase space (e.g. 
energy, momentum, etc.) For some observable function ipo defined on W^, the evolution of 
the value of this observable along the trajectory is given by (pt{v) = ipo{Vt{v)) and ^pt is 
solution to the following dual linear PDF 

dtif-F- Vyif = 0. 

Now let us consider a nonlinear evolution system 

— = Q{f) in an abstract space f £ Ti. 



The formula would be slightly more complicated in the case where F is not divergence-free and would 
involve the jacobian of F, but this is not relevant for the present heuristic discussion. 
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By analogy we define two linear evolution systems on the larger functional spaces -P('H) 
and Cb(^): first the abstract Liouville equation 



dtn+--iQ{f)7r) 



0, vr G P{n) 



and second the abstract equation for the evolution of observables 

of 

However in order to make sense of this heuristic, the scalar product have to be defined 
correctly as duality brackets, and, most importantly, a differential calculus on % has to 
be defined rigorously. Taking % = P{E), this provides an intuition for our functional 
construction, as well as for the formula of the generator G°° below (compare the previous 
equation with formula (|2.8|) 1. Be careful that when Ti = P{E), the abstract Liouville and 
observable equations refers to trajectories in the space of probabilities P{E) (i.e. solutions 
to the nonlinear equation (|2.1|) 1. and not trajectories of a particle in E. Note also that 
for a dissipative equation at the level of Ti (such as the Boltzmann equation), it seems 
more convenient to use the observable equation rather than the Liouville equation since 
"forward characteristics" can be readily used in order to construct the solutions to this 
observable equation. 

Summing up we obtain the following picture for the semigroups: 



Pf on E^/&^ 



observables 



Pn{E) C PiE) 



observables 



- T/v on CbiE^) 
f 



R^ 



\ 



rr on C{P(E)) 



observables 



5f^ on P{E) 



Hence a key point of our construction is that, through the evolution of observables, we 
shall "interface" the two evolution systems (the nonlinear limiting equation and the A^- 



particle system) via the applications vr^ and R 



N 



From now on we shall denote n 



N 



TT, 



2.4. The metric issue. C{P{E)) is our fundamental space in which we shall compare 
(through their observables) the semigroups of the A^-particle system and the limiting 
mean-field equation. Let us make the topological and metric structures used on P[E) 
more precise. At the topological level there are two canonical choices (which determine 
two different sets C{P{E))): 

(1) The strong topology which is associated to the total variation norm, denoted by 
II • ||a,/i; the corresponding set shall be denoted by Ci,(P{E),TV). 

(2) The weak topology, i.e. the trace on P{E) of the weak topology on M^{E) (the 
space of Radon measures on E with finite mass) induced by Ch{E)] the correspond- 
ing set shall be denoted by Ch{P{E),w). 
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It is clear that 

ChiPiE),w)cCt{P{E),TV). 
The supremum norm ||<I>||j;^oo(p(£;)) does not depend on the choice of topology on P{E), 
and induces a Banach topology on the space Cb{P{E)). The transformations vr and R 
satisfy: 

(2.2) lk^^||Loo(£;iV) < WHi^iPiE)) and WR^Wloc^^p^^e)) - W^h^iE^^)- 

The transformation tt^ is well defined from Ch{P{E),w) to Cii{E^), but in general, it 
does not map Cb{P{E),TV) into Cb{E^) since the map 

V e E^ ^ fi^ e {P{E),TV) 

is not continuous. 

In the other way round, the transformation R is well defined from Cb{E ) to Cb{P{E),w) 
and therefore also from Cb{E^) to Ch{P{E),TV): for any G Ch{E^) and for any se- 
quence fk^ f weakly, we have f®^ -^ f®^ weakly, and then R^Mfk) ^ R^Mf)- 

There are many different possible metric structures inducing the weak topology on 
Cb{P{E),w). The mere notion of continuity does not require discussing these metrics, 
but any subspace of Cb{P{E),w) with differential regularity shall strongly depend on this 
choice, which motivates the following definitions. 

Definition 2.4. For a given weight function m : E ^ M4., we define the subspaces of 
probabilities: 

Pn, := {f G PiE); (/,m)<oo}. 

We also define the corresponding bounded subsets for a > 

BPm,a ■■= {f G Pn.; {f,m)<a}. 

For a given constraint function m : E ^ R^ such that the components m are controlled 
by m, we define the corresponding constrained subsets 

^m,m,r := {/ e P„; (/, m) = r}, r G R^ . 

We also define their corresponding bounded subsets for a > 

SPm,m,a,r ■= {/ G SPm,a; (/, m) = r}, 

and the corresponding vectorial space of "increments" 

IPm,m := {/2 - /i; 3r € M^ S.t. /i,/2 G P^.m.r} 

and 

TPm,rr.,a := {/2 - fl\ 3 r G M^ S.t. /i,/2 G BPm,rr.,r,a] ■ 

Let us now define the notion of distances over probabilities that we shall consider. 

Definition 2.5. Consider a weight function mg and a constraint function mg. We shall 
use for the associated spaces of the previous definition the simplified notation Pg for P^ , 

Pg^r for Pm,m,r, BPg^a for BPm,a, SPg^r,a for BPm,m,a,r, I'Pg for XPm,m and XPg,a for 
XP 

We shall consider that a distance dg which 

(1) either is defined the whole space Pg (i.e. whatever the values of the 
constraints) , 

(2) or such that there is a Banach space Q D TPg endowed with a norm 
II • II g such that dg is defined on any Pg^n r G M^, by setting 

y f, 9 G Pg,r, dg{f,g) :=\\g- f\\g. 
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Let us finally define a quantitative Holder notion of equivalence for the distances over 
probabilities that we shall consider. 

Definition 2.6. Consider some weight and constraint functions mg, mg. We say that 
two metrics do ^-nd di defined on Pg are Holder equivalent on bounded sets if there exists 
K £ (0, oo) and, for any a G (0, oo), there exists Ca G (0, oo) such that 

yf,geBPg,a do{f,g)<Ca[diif,g)r, diif,g)<Ca[doif,g)r 

for some constant Ca depending on o > 0. If do and di are resulting from some normed 
spaces Go and Gi we generalize this definition as 

VrGM^, V/,5G^Pg,r,a, do{f,g)<Ca[di{f,g)r, di{f,g)<Ca[do{f,g)r 
for some k £ (0, oo) and some constant Ca depending on a > 0. 



Example 2.7. The choice 

me=mg:=l, mg := 0, || • |b := || • ||mi 
recovers Pg{E) = P{E). More generally on can choose 

me = 1, mg^^iv) := dE{v,vo)'', mg^ := 0, || • Ib^ := || • d^lv, wo)^||mi- 

For fes > and ki > k2, the spaces Pg^, and Pg^, are topologically uniformly equivalent 
on bounded sets of Pg^ . 

Example 2.8. There are many distances on P{E) which induce the weak topology, see for 
instance [64J . In the next section, we present some of them which have a practical interest 
for us, and which are all topologically uniformly equivalent on bounded sets of P{E) in 
the sense of the previous definition, with the choice of a convenient (strong enough) weight 
function. 



2.5. Distances on probabilities. Let us discuss some well-known distances on P( 
(or defined on subsets of P{M.'^)) which shall be useful for the sequel. These distances are 
all topologically equivalent to the weak topology a{P{E), Cb{E)) on the sets BPk^a{E) for 
k large enough and for any a G (0, oo), and they are all uniformly topologically equivalent 
(see [BHl [T7] and section I2.5.6P . We refer to [6ll Ell [T7] and the references therein for more 
details on these distances. 

2.5.1. Dual-Holder -or Zolotarev's- distances. Denote by dE a distance on E and let us 
iix vo £ E (e.g. uo = when £' = M*^ in the sequel). Denote by Cq'*(£^), s G (0, 1) (resp. 
LipQ(ii^)) the set of s-Holder functions (resp. Lipschitz functions) on E vanishing at one 
arbitrary point vq £ E endowed with the norm 

b], := sup \v(y)-^(f\ ^,^ (0,1], [^]^,^ := [^],. 

We then define the dual norm: take mg := 1, mg := and Pg{E) endowed with 

"^Lg^Pg, [g -/]:■■= sup ^^ff^. 
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2.5.2. Wasserstein distances. Given q G (0,cxd), define Wq on 

Pg{E) = Pg{E) := {/ G PiE); {f,dE{-,Vof) < 00} 
by 

yf,gePg{E), Wg{f,g):= inf / dE{x,yyU{dx,dy), 
where II{f,g) denote the set of probability measures 11 G P(E x E) with marginals / and 

g- 

U{A, E) = f{A) and U{E, A) = g{A) for any Borel set Ac E. 
Note that for Vi, V2 G E^ and any g G [1, 00), one has 

\ 1/9 



(i))" 



( \ ^ 
(2.3) VF,K,0 =d£,(£;iv/6,)(Fi,y2) := min ^ d^ ((^i)., (^^2). 

\ 4 = 1 

and that 

V/,5GA(i^), VFi(/,5) = [/-5]i= sup (/-5,95>. 

¥'GLipo{-B) 

2.5.3. Fourier-based norms. Given E = M , ?tt,cj^ := \v\, mgj := 0, let us define 
V/GXPg,, ll/lle^ = I/I, := sup ^, SG(0,1], 



where / denotes the Fourier transform of / defined through the expression (when say 

m--= I f{x)e-'-<dx. 
Similarly, given E = W^, mg^ := |up, mg^ := v, we define 

V/GXPg„ ||/||g^ = I/I, := sup ^, sG(1,2]. 

Obviously high-order versions of this norm could be defined similarly by increasing the 
number of constraints. However we shall see in the next subsubsection how to extend this 
notion of distance without constraints. 

2.5.4. More Fourier-based norms. More generally, given E = W^ and /c G N*, we set 

mg:=\v\^, mg := (v")^gp^dj^|<fc_i 



with |q| = «! + ••• + arf and 

^"=(<S...,^;'/), a = («!,..., ad)GN' 
and we define 



V/GXPg, ||/||g = |/|, :=sup^, sG{0,k]. 



In fact, we may extend the above norm to M^{W^) in the following way. We first define 
for 

/ G mI_^{R'^) and a G N"*, |a| = ai + h Od < A: - 1 

the following moment: 

Ma[f]-= I V^fidv). 
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Consider a fixed (once for all) function x ^ C^(M ) (compact support), such that x = 1 
on the set {v E W^, \v\ < 1}. This implies in particular 



T-'{x){v)dv = x{0) = l. 
Then we define A4fc[/] through its Fourier transform 



|a|<fc-l 



a! 



a! := ail . . . a^! 



Note that this is a mollified version of the {k — l)-Taylor expansion of / at ^ = 0. Then 
we may define the semi-norms 



\k ■■= sup 



fio-Mk[fm 



and the norms 



lel' 



:=|/U+ Yl l^"[/]l 

aGN'i,\a\<k-l 



2.5.5. Negative Soholev norms. Given s £ {d/2, d/2 -\- 1/2) and 

E = 'R'^, mg^ := \v\, mg^ := 
we define the following negative homogeneous Sobolev norm 



V/GXPg,, 



H- 



m) 



i^h 



L2 



Observe that probabilities are included in the corresponding non homogeneous negative 
Sobolev space H~'^{M.'^) as soon as s > d/2. 
Similarly, given s € [d/2 + 1/2, d/2 + 1) and 



we define 



E 



yf^IPg,, 



, mg2 ■■= \v\ , mg2 



52 



H- 



fiO 



lei^ 



L2 



2.5.6. Comparison of distances when E = M. . All the previous distances are Holder equiv- 
alent on bounded sets in the sense of Definition 12.61 Precise quantitative statements of 
these equivalences are given in Lemma 14.11 in Section 14.11 

2.6. Differential calculus in probability spaces. We start with a purely metric defi- 
nition in the case of usual Holder regularity. 

Definition 2.9. Given some metric spaces Gi and Q2, some weight function 

A 1^1^ [l,+oo), 
we denote by 

ucA{gi,g2) 

the weighted space of uniformly continuous functions from Qi to ^2 7 that is the set of 
functions S : Qi ^ Q2 such that there exists a modulus of continuity uj so that 

(2.4) V/i, /2 G §1 dg, (5(/i),5(/2)) < A(/i,/2)^(dg, (/l,/2)), 

with 

A(/i,/2):=max{A(/i),A(/2)} 
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and where dg^, denotes the metric of Qk- Note that the tilde sign in the notation of the 
distance has been removed in order to present unified notation with the next definition. 
For any rj G (0, 1], we denote by 

the weighted space of functions from Qi to G2 with r^-Holder regularity, that are the 
uniformly continuous functions for which the modulus of continuity satisfies u}{s) < C s'^ 
for some constant C > 0. We then define the semi-norm 

['^]c°-''{gi,g2) ^°^ S (^C^/'{Gi,g2) 
as the infimum of the constants C > such that ()2.4p holds with oo{s) = C s^ . 

We now define a first order differential calculus, for which we require a norm structure 
on the functional spaces. 

Definition 2.10. Given some Banach spaces ^1, ^2 and some m,etric sets Gi, Q2 such that 

Qi —Qi <Z Gi, i = 1,2, 

some weight function 

A: ^1 ^ [1,00), 
we define 

UC/^ iGi,Gi;G2,Q2 

(later simply denoted by UCj^{Gi',G2)), the space of continuously differentiable functions 
from Gi to G2, whose derivative satisfies some weighted uniform continuity. 
More explicitely, this is the set of uniformly continuous functions 

S : Gi ^ G2 
such that there exists a uniformly continuous function 

DS:Gi^CiGi,G2) 

(where £(^1,^2) denotes the space of bounded linear applications from Gi to G2 endowed 
with the usual operator norm), some constants Cj > 0, i = 1, 2, and some modulus of 
differentiability fi (that is a function il : M-|- — ?• M-|- such that il(s)/s — )• when s — )• 0) so 
that for any fi, f2 ^ Gi'- 

(2.5) l|5(/2)-5(/i)||5^ < C7iA(/i,/2) 11/2 -/lib, 

(2.6) \\{DS[h],f2-fi)\\g^ < C2A(/i,/2) 11/2 -/lib, 

(2.7) ||5(/2)-5(/i)-(i^5[/i],/2-/i)|b^ < C3A(/i,/2)f](||/2-/i|bJ. 
For any ry G (0, 1], we also denote by 

C'a'' (^i>^i;^2,^2 

(later simply denoted by C^'''(^i; ^2)), the space of continuously differentiable functions 
from Gi to G2 J whose derivative satisfies some weighted ?7-H61der regularity, which is the set 
of continuously differentiable functions such that the modulus of differentiability satisfies 
r2(s) < C s^^^ for some constant C > 0. We define respectively Cf , Cf, Cf, as the 
infimum of the constant Ci , C2 , C3 > such that respectively ()2.5p , (12. 6p , ()2.7p with 
Q{s) = C3 s^"'"'' hold. We then define the sem,i-norm,s 

[S]^oi := Ci , [S]^io := C2 , [S]c\'V := C3 , ||5||^i,^ := Q + C2 + C3 . 

^A ^A ^A ^A 

In the sequel we omit the subscript A or we replace it by the subscript b in the case when 
A = l. 
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Remark 2.11. In the sequel, we shall apply this abstract differential calculus with some 
suitable subspaces Gi C P{E). This choice of subspaces is crucial in order to make rigorous 
the intuition of Griinbaum [36] (see the — unjustified — expansion of Hf in [36]). It is 
worth emphasizing that our differential calculus is based on the idea of considering P{E) 
(or subsets of P{E)) as "plunged sub-manifolds" of some larger normed spaces Qi. Our 
approach thus differs from the approach of P.-L. Lions recently developed in his course 
at College de France [46] or the one developed by L. Ambrosio et al in order to deal 
with gradient flows in probability measures spaces, see for instance [2]. We develop a 
differential calculus in probability measures spaces into a simple and robust framework, 
well suited to deal with the different objects we have to manipulate (1-particle semigroup, 
polynomial, generators. . . ). And one of the main innovations of our work is the use of 
this differential calculus to state some "differential" stability conditions on the limiting 
semigroup. Roughly speaking the latter estimates measure how this limiting semigroup 
handles fluctuations departing from chaoticity. They are the corner stone of our analysis. 



This differential calculus behaves well for composition in the sense that for any given 

4'(^~i>^~2) and VGCl'Jl 



UGcl:;j{gug2) and vgcI;'^ (02,^3) 



there holds 



s:=vouecl:^{gi,g3) 



for some appropriate weight function Ag. We conclude the section by stating a precise 
result well adapted to our applications. 

Lemma 2.12. For any given 

UecY'{Gi,G2) and VeC^''^{G2,G3) 
there holds 

S:=VoUe Cl;,\,{gi,g3) and DS[f] = DV[U{f)] o DU[f]. 
More precisely, there holds 



and 



[5]^o.i < [V]co,i Mco^i, [5]^i,o < [V]ci,o [^/l^i.o 



[5]^!,. < [v]ci.o [uy, + [vic-i.. m'll- 



When further V € C^'^{G2,Q3), we also have 

S:=VoUeC];^igi,g3) 

with 

[S]c\v < [V]ci.o [UU„ + [V]cM [^^]'o,(i+.)/2. 

a2 "-"A O^ 

Proof of Lemma \2.1IA For /i , /2 G Qi we have 

UU2) = UUi) + [DUifi], (/2 - /i)) + TZu (/i, /2) 
with 

\\{Dl([pi],if2-fl))\\g^ < [^]^l,oA(/i,/2) ||/2-/l||g, 

and 

\\nu{fl,f2)\\g,<[liy^vAifl,f2) \\f2-fl\\lt\ 

and a similar Taylor expansion holds for V: for gi,g2 S 02, 

V(52) = V(<7i) + {DV[gi], (52 - 5i)) + 7^v (51,92) 
with 

\\{DV[gi],ig2-gi))\\g.^ < [VJcm \\g2 - gi\\g. 
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and 

||7^v(ffl,<?2)lb3<[v]cl,. 1152-51 list'- 

We then write 

\\S{f2) - 5(/i)lb3 ^ Mco^ \mf2)-U{h)\\g^ < Mco. [U]^o.iA{h,h) 11/2 - /2||g, 

which imphes 

[5]^o,i < [V]co,i [U]^o,i 

and 

s{f2) = {Vou){f2) = v{u{h) + {DU[hi{f2-h))+nu{fuf2) 
= v{u{h)) + nv{u{f2)M{h)) 

+ (^DV[U{h)l{{DU[h],{f2- h)) +'Ru{hj2))\ 
which imphes 

{DS[h], ih - h)) = {DV[U{h)], ( {DU[h], (/2 - /i)> 
and 

Finahy we estimate the remaining term: 

\S{f2)-S{h)-{DS[h],{f2-h))\ 

= \TZv{U{f2)M{h)) + {DV[U{h)inu{h,f2))\ 

< (a(/i,/2)^+'' [V]ci,.[^^]Ji!+A(/i,/2) [V]ci,o [Z^]^!,.) ||/2-/i||Jt'- 
We hence conclude that 

[sy.. < [v]ci,o [uy„ + [v]ci.. [u]]+i . 

The variant 

a2 a O^ 

is easily obtained by estimating instead 



|2 



\nv{U{f2),U{fl))\ < [V]cM l|Z^(/2)-Z^(/lJ||g, 

< [V]cmA(/i,/2)' [^]^o,(i+.)/2 ll^-Mlit'- 

D 

2.7. The pushforward generator. As a first application of this differential calculus, let 
us compute the generator of the pushforward limiting semigroup. 

Lemma 2.13. Given some Banach space Q and some probability space Pg{E) (see Defini- 
tions \2.4^2.5\} associated to a weight function m and constraint function m, and endowed 
with the metric induced from Q, then for some 5 G (0, 1] and some a G (0, 00) we assume 
that for any a G (a, 00); 

(i) The equation (|2.ip generates a semigroup 

which is 5-Hoder continuous locally uniformly in time, in the sense that for any 
T G (0, 00) there exists Cr G (0, 00) such that 

yf,gGBPg,,a, sup ||5fV-Sf^<7L, <Cr \\f-g\\g,- 

te[o,r] 
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(ii) The application Q is hounded and 6-Hdlder continuous from BPg^a into Q. 
Then for any a G (a, oo) the pushforward semigroup T^ defined by 

V/ G BPg,a{E), CD G Ua{BPg,aiE)), Trmif) ■= ^ {8^(1)) 

is a Co-semigroup of contractions on the Banach space U Ch{BPg ^a{E)) . 

Its generator G°° is an unbounded linear operator on U Cb{BPg ,a{E)) with domain 
Dom{G°°) containing UC^{BPg^a{E)). On the latter space, it is defined by the formula 



(2.8) 



VcD G UCl{BPg,a{E)), V/ G BPg,a{E), (G^^) (/) := {D<!>[f],Q{f)) . 



Remark 2.14. Note that the restriction to uniformly continuous functions ^ on probabihty 
spaces win be harmless in the sequel for two reasons: first in most cases our choice of 
weight, constraints and distance yields a compact space BPg^a{E), and second and most 
importantly we shall only manipulate this pushforward semigroup for functions $ having 
at least Holder regularity. 

Proof of Lemma \2.1'J[ The proof is split in several steps. 

Step 1. We claim that for any /o G BPg^aiE) and r > the application 

Sifo):[0,T)^Pg, f^5f^(/o) 

is right differentiable in t = with 

5(/o)'(0+) = Q(/o). 

Denote ft := S^^fo. First, since ft G BPg^a for any t G [0,r] and Q is bounded on 
BPg^a{E) (assumption (ii)), we deduce that uniformly on /o G BPg,a[E) we have 

ft 



(2. 



Wft-fi 



Ollg 



Q{fs) ds 



< Kt. 



We then use the previous inequality together with the fact that Q is ^-Holder continuous 
(assumption (ii) again), to get 

ft 



\\ft-fo-tQ{fo)\\g 



{Q{fs)-Q{fo))ds 







= ^ Wfs- foWg ds 
Jo 

ft 4 

< L / {KsYds = LK^- 
Jo 1 



1+5 



+ 5' 



which implies the claim. 

Then the semigroup property of (S^^) implies that t i— t- /< is continuous from R+ into 
Pg{E) and right differentiable at any point. 

Step 2. We claim that (T^°°) is a Co-semigroup of contractions on UCb{BPg^a{E)). 

First for any $ G UCb{BPg^a{E)), we denote by a;$ the modulus of continuity of $. We 
have thanks to the assumption (i): 

vtG[o,T], \{Tr^){g)-{Tr^){f)\ = |i>(sf^(5))-^(5f^(/))| 

< u^(cA\g-f\\gX 
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SO that T^^ G UCb{13Pg^^a{E)) for any t G [0,r], and then, by iteration, for any t > 0. 
Next, we have 

iirni = sup \\Tm = sup sup \^ {sru))\ < 1 

since 

||$|| = sup |<^(/)|- 

feBPg,aiE) 

FinaUy, from 1^, for any $ G UCb{BPg^a{E)), we have 

l|rt°°$-$||= sup |$(sf^(/))-$(/)| <w$(i^t)^o. 

fel3Pg,,{E) 

As a consequence (Tj°°) has a closed generator G°° with dense domain 

Dom(G°°) C UCb{BPg,a{E)) 
(see for instance [62l Chapter 1, Corohary 2.5]). 

Step 3. We shall now identify this generator, at least on a part of its domain. Let us 
construct a natural candidate provided by the heuristic of Remark 12 .31 Let us define G°°$ 

by 

V$ G UCl{BPg,a{E)), V/ G BPg^aiE), (G^^W := {D^f],Q{f)) . 
The right hand side is well defined since 

D^f]£C{g,R)=g' and Q(/) G g. 
Moreover, since both applications 

f^D^if] and f^Qif) 
are uniformly continuous on BPg^a{E), so is the application 

/^(G'-$)(/). 
Hence (5°^$ G UCb{BPg,a{E)). 

Step 4. Finally, by composition, for any fixed $ G f/Cj^(i3Pg,a(ii;)) and / G BPg^a{E), 
the map 

t^TrHp) = <^os^'^{p) 

is right differentiable in t = and 

|(rr<I>)(p)k=o := |(<I>o5(/.)(t))|,=o 

Z)ci>[5(/,)(0)],|5(p)(0) 



= {D<i>[p],Qip)) = [G^<^){p), 
which precisely means that $ G Dom(G°°) and that (12. Sh holds. D 
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2.8. Duality inequalities. Our transformations vr^ and R^ behave nicely for the supre- 
mum norm on Cb{P{E),TV), see ()2.2p . More generahy we shall consider ''^duality pairs'^ 
of metric spaces as follows: 

Definition 2.15. We say that a pair {J^,Pg) of a normed vectorial space J-" C Ch{E) 
endowed with the norm || • \\jr and a probability space Pg C P{E) endowed with a metric 
dg satisfy a duality inequality if 

(2.10) v/,gGPg, V<^G^, \{g-f,ip)\<Cdg{f,g)M^, 

where here (•, •) stands fot the usual duality brackets between probabilities and continuous 
functions. In the case where the distance dg is associated with a normed vector space G, 
this amounts to the usual duality inequality \{h,ip)\ < \\h\\g \\(p\\j^. 

The "compatibility" of the transformation R for any such pair follows from the mul- 
tilinearity: if T and Q are in duality, T C Ch{E) and Pg is endowed with the metric 
associated to || • \\g, then for any 

LP = ^piX ■■■ X Lp]^ £ jr8)Af^ 

the polynomial function R^ in Cb{P{E)) is C^'^(Pg,M). Indeed, given /i,/2 G Pg^, we 
define 

g ^ M, h^ DR'^ [/i] (h) := f; ( n ^/i' '^^•) I ^^' '^^) ' 

and we have 

<(/2)-<(/i) = e( n ^/2,¥^fc)) (/2-/i,v'.) ( n ^/I'V'fc) 

j=i yi<fc<j y \i<k<e 

and 

<(/2) - < (/l) - I)< [/l] (/2 - /l) = 

= E I n (/a'V'fe) I {f2-fi,^j) I n (/I'^'i) I (/2-/i,^.) ( n ^/i'^^^) 

l<j<j<7V \l<fc<j J \j<k<i J \i<k<e 

Hence for instance R^ G C^'^(Pg;R) with 

|<(/2)-<(/l)| < A^II</'IU(L-)-- 11/2 -/lib, 
\DR^[h]ih)\ < N ||(/.||_^^(^^)^_, \\h\\g, 

and 

(2.11) \R^ih) - <(/l) - I)<[/l](/2 - /l)| < ^^^^~ ^ V ll.F2^(L°^r-^ 11/2 - flWl , 

where we have defined 

ll¥'llj-fccg)(L°°)Jv-fc := . max H^Jj, ||j- . . . ||v3ij|j- TT ||c^j||ioo(£;) 

^ ' ll,...,lk distincts in [lljA^I] \ "^ "^ 

Remarks 2.16. (1) It is easily seen in this computation that the assumption that if is 
tensor product is not necessary. In fact it is likely that this assumption could be 
relaxed all along our proof. We do not pursue this line of research. 
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(2) The assumption T C C\){E) could also be relaxed. For instance, when 

T := Lipo(^) 

is the space of Lipschitz function which vanishes in some fixed point xq £ E, Q is 
its dual space, and 

Pg ■■= {/ G Pi{E); {f,distE{;Xo)) < a} 
for some fixed a > 0, we have i?<^ G C^''^{Pi{E);R) with 

fP^l ^AT^N-2\\,„\\ fO^l N{N - 1) jy_]^ 

N-2 



or equivalents i?^ € C'^\Pi{E)-^) with A(/) := ||/||^;;t 



In the other way round, for the projection tt it is clear that if the empirical measure 
map 

V (^E^ ^^l^£ P{E) 

belongs to C^'^{E^ , Pg) for some norm structure Q, then by composition one has 

(2-12) lk^(^)|lc"=>'7(£;iV;K) < C'tt \Mc'^,v{Pg) ■ 

However the regularity of the empirical measure of course heavily depends on the choice 
of the metric Q. 

Example 2.17. In the case T = {Cb{E),L°°) and G = {M^{E),TV), ([232]) trivially holds 
with C'^''' replaced by C^. 

Example 2.18. When J^ = LipQ(£^) (Lipschitz function vanishing at some given point vq) 
endowed with the norm ||(/>||i,jp and Pg{E) (constructed in Subsubsection [23^ is endowed 
with the Wasserstein distance Wi with linear cost, one has (I2.12p with /c = 0, r/ = 1: 

|$(/.^) -$(^^)| < ||$||co,i(P,)TVi(/if,/iy) < ||$||c,o,i(P,)||X-y||,i, 

where we use ()2.3p . which proves that 

IK {^)\\co^{EN) < ll^llco>i(^6)' 
when E is endowed with the £^ distance defined in (12. 3p . 

3. The abstract theorem 

3.1. Assumptions of the abstract theorem. Let us list the assumptions that we need 
for our main abstract theorem. 

(Al) Assumptions on the A^-particle system. 

G and T^ are well defined on Cb{E ) and invariant under permuta- 
tion, and they satisfy the following moment conditions: 
(i) Propagation of integral moment bound: There exists a weight function 
mg-^^, a time T E (0, oo] and a constant Ci^t £ (0,oo), possibly de- 
pending on T and mg-^^, but not on the number of particles A^, such 
that 

ViV>l, sup (/f,M4 )<C™,^,T. 

0<t<T \ ^ / 

(ii) Support moment bound at initial time: There exists a weight function 
mg^ and a constant C^ G (0, +oo), possibly depending on the number 
of particles A^, such that 



Supp/o^ c{ve i?^;M4^(y) < c4^,o} ■ 
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Note that the name of the weights mg^ and mg^ in (Al) above are chosen for coherence 
with the functional spaces in the other assumptions. 

(A2) Assumptions for the existence of the pushforward 

semigroup. 

Given some Banach space Qi and some probabihty space Pg^{E) (see 
Definitions I2.4||T^ associated to a weight function mg^ (as in (Al)-(i)) 
and a constraint function mg^ , and endowed with the metric induced from 
Qi, then for some 6 G (0, 1] and some a G (0, oo) we assume that for any 
a G (a, oo): 
(i) The equation (|2.ip generates a semigroup 

5f^ : BPg,,a ^ l3Pg,,a 

which is 5-Hoder continuous locahy uniformly in time, in the sense 
that for any r G (0, oo) there exists Cr G (0, oo) such that 

yf,gGBPg,,a, sup ||SfV-5f^5L, < Cr \\f - 9\\g,- 

te[o,r] 

(ii) The application Q is bounded and 5-Holder continuous from BPg^^a 
into Qi. 

So in particular the semigroups S^ , T^ , S^^ and T^°° are well defined as well as the 
generators G^ and G°°. 

We then need the key following consistency assumption. It intuitively states that the A^- 
particle approximation of the limiting mean- field equation is consistent. More rigorously 
this means a convergence of the generators of the A^-particle approximation towards the 
generator of the limiting pushforward semigroup within the abstract functional framework 
we have introduced. 

(A3) Convergence of the generators. 

In the probability metrized set Pg^ introduced in (A2) (associated to 
the weight function mg^ and constraint function mg^ ) we define 

Rg, := {r G M^; 3 / G Pg, s.t. mg, (/) = r} . 

We also define a weight function 

m'g^ < C mgi 

possibly weaker than mg^ and we define the associated weight on the dis- 
tribution: 

Ai(/):=(/,m'g^>. 
Then we assume that for some function 

e{N) ^0 as iV ^ oo, 
the generators G^ and G°° satisfy 
V<I>G fl Clf(Pg„r;M) 

reRg^ 



MZJ (G^7r;v-vr^G-)<i> < e{N) sup [^U.,, ., 
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where MI^„ is defined by 



N 



Note the following aspect, which shall be a crucial source of difficulty in the application 
of the following abstract theorem: the loss of weight in the consistency estimate 
(A3) has to be matched by the moment bounds propagated on the A^-particle 
system in the assumption (Al)-(i). (In fact the loss of weight in the consistency 
estimate mg^ can be even slighlty higher than the one from the stability estimate Ai, see 
the above relation.) 

More specifically, the best we are able to prove uniformly in N on the A^-particle system 
are polynomial moment bounds. This thus constraints the kind of loss of weight we can 
afford in the following stability estimate. 

We now state the second key stability assumption. Intuitively this corresponds to the 
abstract regularity that needs to be transported along the flow of the limiting mean-field 
equation so that the fluctuations around chaoticity can be controlled. More rigorously 
this means some differential regularity on the pushforward limiting semigroup, which cor- 
responds to some differential regularity on the limiting nonlinear semigroup according to 
the initial data and in probability spaces. 

(A4) Differential stability of the limiting semigroup. 

We consider some Banach space Q2 D Gi (where Qi was defined in (A2)) 
and the corresponding probability space Pg^ (E) (see Definitions 12. 4112. 5p 
with the weight function mg^ and the constraint function mg^, and en- 
dowed with the metric induced from Q2- 

We assume that the fiow S^'^ is C^'''(Pg^^ri Pg2) ^^^ ^^J ^ ^ ^Gi i^ the 
sense that there exists C^ > such that 

where rj E (0, 1) is the same as in (A3), where 

l + V 



{v,v") = {r]A) or {r]',r]")= il, 



2 

and with the weight (recall that Ai was defined in (A3)) 

A2=A,^. 

We finally state a weaker stability assumption on the limiting semigroup. It shall be used 
intuitively for proving that initial error in the law of large number when approximating 
a probability by empirical measures is propagated by the limiting semigroup. The reason 
for dissociating this assumption from the previous one is because we need some fiexibility 
in the independent choice of distances for these two assumptions. 

(A5) Weak stability of the limiting semigroup. 

We assume that, for some probabilistic space Pg^{E) associated to the 
weight function mg^ (as in (Al)-(ii)), a constraint function m^g and some 
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metric structure dg^ , for any a,T > there exists a concave modulus of 
continuity @a,T such that we have 

yfl,f2eBPg,,a{E) 

sup dg, (5f^(/l), 5f^(/2)) < ea,T {dg, (/i, /a)) . 

[o,r) 



Observe that in the latter assumption, we require that /i, /2 € BPg^^aiE) which implies 
in particular that we require the bounds 



Mmo, (/l) 



/ /i "1^3 du < a, Mmg (/2) = / /2 "iga ^^ < «• 



When applying the assumption to some empirical measure for one of the argument /i or 
/2, this shall induces the requirement of controlling pointwise terms like M^ . This is 
the reason for the assumption (Al)-(ii). 

3.2. Statement of the result. We are now in position to state the main abstract result. 
This result can be considered intuitively as a convergence in approximation theory, in 
the sense of proving that approximation errors between the A^-particle system and the 
limiting mean-field system are propagated along time without instability amplification 
mechanism. More specifically the approximation error means in the present context some 
kind of distance between the discrete A^-particle system and the limiting mean-field system, 
within our abstract functional framework. This result implies in particular the propagation 
of chaos. 

Theorem 3.1. Consider a family of N -particle initial conditions 



N\ 



f^' ePsyUE''), N>1 



and the associated solutions 

ff = 5f (/o^) . 
Consider a 1-particle initial condition /q G P{E) and the associated solution 

ft = Sr{fo) 

of the limiting mean-field equation. 

Assume that the assumptions (A1)-(A2)-(A3)-(A4)-(A5) hold for some spaces Pg^, 
Qk and Tk, fc = 1, 2, 3 with Tk C Cb{E), and where Tk and Qk are in duality (that is (I2.10|) 
holds). 

Assume also that the 1-particle distribution satisfies the moment bound 



M, 



CXD. 



mea (/o) = if o,mg^) < + 

Then there is an explicit absolute constant C € (0,oo) such that for any N,i gW , with 
N > 21, and for any 



ip = ipi (^ ip2 



^i G (J-in-FanJ-s)^ 



we have 
(3.1) 



sup 

[0,T) 



<C 



5f(/o^)-(5f^(/o))^^),V.) 



.2 \m 



N 



+ Cmg^^TCT£{N)f \\f\\jrli^(L'^)e- 



+^ \M:Fs^{L^y-^ Qa'^,T [y^dg,^{T^pfi',Sfo)j 
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where a^ > depends on Cg, ^ and Mg^{fQ), and where Wdg stands for an abstract 
Monge-Kantorovich distance in P{Pg^{E)) (see the third point in the next remark) 



(3.2) 



m,^{^^f^.5f,)= I dg,{^Ji'^,fo)f^'{dV). 



Remarks 3.2. (1) The spirit of this method is to first treat the A^-particle system as a 

perturbation (in a very degenerated sense) of the hmiting problem, and second to 
minimize assumptions on the many-particle systems in order to avoid complications 
of many dimensions dynamics. 

(2) In the applications the worst decay rate in the right-hand side of (j3.ip is always 
the last one. This last term controls two kind of errors: (1) the chaoticity of the 
initial data, that is how well /g ~ /q , (2) the rate of convergence in the law of 
large numbers for measures in the distance dg^ . 

(3) Let us discuss more the meaning of this last term and the related issue of sampling 
by empirical measures in statistics (see also Section H]). 

Following the abstract definition of the optimal transport Wasserstein distance 
we define 



V ^1 , /i2 G -P [Pg-i ) , Wdg., ifJ-i , fJ-2) 



inf 



dg-i {pi,P2) TT{dpi,dp2 



where n(/ii,/i2) denotes the probability measures on the product space Pg.^ x Pg^ 
with first marginal /Ui and second marginal fi2- In the case when fi2 = <5/o then 

n(/^i,5/o) = {w«"^/o} 
has only one element, and therefore 

W.,3(7r^/o^,5/o) 



inf / dg^{pi,p2) TT {dpi , dp2 ) 

^6n(^^/,f,5^j7Pg3XPg3 



EN 



dgAf^vJo) fo'^idV) 



which explains the notation (|3.2p . We simply denote in the tensorized case: 

Comparisons of the W^ functionals and estimates on the rate 

^dU) ^0 as N ^oo 
depending on the choice of the distance d are discussed in Subsection 14.21 
3.3. Proof of the abstract theorem. For a given function 

VPG(J-inJ-2nJ-3)®^ 

we break up the term to be estimated into three parts: 



(5f(/o^)-(5r(/o))^''),^^ 



(^N-e 



< 



:!NffN 



(^N-e 



N 



< (^6r(/^'),v9®i-— ^^-^5f(/^^),i?^o/i{) 

+ |(/o^,r,^(4 o p^)) - (/„^, (rr<) o p^ 

We deal separately with each part step by step: 



(=: Ti) 

=--T2) 
{=--T,). 
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7i is controlled by a purely combinatorial arguments introduced in |36J. Roughly 
speaking it is the combinatorial price we have to pay when we use the injection 
TT^ based on empirical measures. 

Ts! is controlled thanks to the consistency estimate (A3) on the generators, the dif- 
ferential stability assumption (A4) on the limiting semigroup and the propagation 
of integral moment bounds (Al)-(i). 

Ts is controlled in terms of the chaoticity of the initial data thanks to the weak 
stability assumption (A5) on the limiting semigroup and the support moment 
bounds at initial time (Al)-(ii). 



Step 1: Estimate of the first term 7i. Let us prove that for any t > and any N > 2i 

there holds 

e2 



(3.3) 



Ti 



5f(/o),^®l 



(giN-i 



qN f fN\ ut ^ ,,N 



< 



2^ I|9'IIl°°(£;0 



N 



Since S^ (f^) is a symmetric probability measure, estimate (13. 3p is a direct consequence 
of the following lemma. 



Lemma 3.3. For any cp G Ci,{E^) we have 



(3.4) 



yN>2£, 



f(^l 



(gN~e 



ttatR 



Sy 771 



< 



2^ I|¥'IIl°°(£;«) 



N 



where for a function (p G C\,{E ), we define its symmetrized version (psym cls: 



■' sym 



T©^ ^ "^^ 



aee^ 



where we recall that 6 is the set of N -permutations. 
As a consequence for any symmetric measure we have 



(3.5) 



V/^ e PsUE""), if'^R'M)) - (/^,^> 



< 



2^ IIv5|Il°°(£«) 



N 



Proof of Lemma \3.3[ . This lemma is a simple and classical combinatorial computation. We 
briefly sketch the proof for the sake of completeness. 
For a given i < N/2 we introduce 



AN,e:={{ii,...,ii)e[\l,N\]' : ^ k ^ k' e [\l,i\], ik + iy } 



and 



Sat,, :=A^,, = {(zi,...,i£)E[|l,iV|]^}\A^,£. 
Since there are N\/{N — 1)1 ways of choosing i distinct indices among [|1, A''|] we get 
\Bni\ , N\ 



m 



(N -i)\N^ 



e-i 



N 



'i-i 



l-exp(^log(l--^ 
\i=o ^ 



< 1-exp -2j;- <2^ 



j=0 



i=0 



N - N' 



where we have used 

Vj;G [0,1/2], log(l-x)>-2x and Vx G M, e""" > 1 - x. 
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Then we compute 



N 



^Um^) = ^ E v^(-- 



V 



^Viij 



(«l,...,i£)GA]v,f (ii,...,if)e-Bjv/ 



o-e6]v 
We now use the same estimate 



1 



AT! 



< 



{N -ey.N^ - N 



as above to get 



R. 



i (..N 






ai^&N 



ip(E> 1 



(^N-e 



+ 



sym 



2£2 

AT 



m\L° 



and the proof of (13. 4p is complete. 

Next for any /^ G PsymiE^) we have 

and (13. 5p trivially follows from ([37 



n 



Step 2: Estimate of the second term 72- Let us prove that for any t £ [0, T) and 

any N >2i there holds 



(3.6) 



Ta := 



^Af rpN I Tji ^ , N\\ fN (rpOO Tji \ ^ , N 

Jo . ^t [Ru:>°f^Vl) - {Jo A^t Ru,} °fJ'V 



< Cmg^^T C^ e{N) I \W\\jr2^(^L-^Y-2. 

We start with the following algebraic identity 

tSn - TTNTr = - /* ^ {Tf^s ^N Tr ) ds = I' Ti[, [G^vr^ - 7rr,G^] T^ ds. 
We then use assumptions (Al)-(i) and (A3) and we get for any t £ [0,T) 



< 



Mi^,^ 5il, (/o^) , (M^^^ ) [G^vr^ - nr,G-] (t^R^ 



ds 



< I sup (/f ,M„^ 

,0<i<T 



M, 



A 



(3.7) 



<Cmc,T^iN) sup 



T ^ 



[G vttv - vrivG°°] \T^R^ 
ds. 



ds 



L°°(Ejv) 



TOO D-^ 



^i'l^C^SLr) 
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Now, let us fix r G Hg^ . Since 



■<i.»?/' 



Tr{%) = % o 5^ with 5^ e C^A^iPgur; Pg.) 

and i?^ E ^'^'^(Pgj) because C/9 G Tf (see subsection 12. Sp . we can apply Lemma [2. 121 and 
use assumption (A4) to obtain 



l,Vf 



TriRf,) G C^f (Pg„ 



with 



T^iK 



cj::(p,.,r)-V^^'^^^A^'~^o,.,Po,) 



+ [s, 



NLl l+V 



-^0,v 



CZ' iP5,.r,Pg^) 



r: 



C^'^iPg^) 



and 



A2 = Ai 



^,1/(1+'?') 

Note that the two different sets of indices 

(r?',r/") = (r?,l) or (r/',r/") = (1, (1 + ??)/2) 

in Lemma 12.121 correspond to the two cases discussed in Lemma 12.121 
We then deduce thanks to (12. lip and assumption (A4): 

(3-8) / [T^{Rl,)]^i,r,i^p^^_^'j ds < Cfl^ (^||(/?||j-255(icx,)f-2 + \W\\r^<^(L^Y~-ij 

Then we go back to the computation (13. 7p . and plugging (13. 8p we deduce (13. 6p . 



Step 3: Estimate of the third term Tjj. Let us prove that for any t > and N > (. 

we have 



T.:-- 



< 



where &a,T was introduced in assumption (A5), and a = a^ is defined by 



,A^ 



max 



{Mg,{fo) , C4^,o} 



where C^ q was introduced in assumption (Al)-(ii). 
Assumption (Al)-(ii) indeed implies that 



Supp/o^ C /C := I F G M'^^ s. t. M^^ {^^) = 1 f2^g.Av^) < C^,} • 



Hence we are in position to apply (A5) for the functions /o and fiy on the support of 
f^ since Mmg (/o) is bounded by assumption, and M 
restricting to V G IC thanks to the previous equation. 

Let us also rec; 

We then write 



f^ since Mmg (/o) is bounded by assumption, and Mg^{^y) is bounded by C^ q when 
Let us also recah that i?^ G C°'^(Pg3,M) because c/? G J"f^. 



% = |(/o^<(^f"(/^^)))-(/o^4(^f'(/o))) 

We now apply (A5) to get 

Vt G [0,r], dg, (5f^(/o),5f^(/xC^)) < G,^,r (dSa (/o,/^y)) 
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on the support of /q^, and therefore 

We can further rewrite the last equation in the following way by using the concavity of 
the Q^N X function: 

which concludes the proof of this step. 

The proof of Theorem 13.11 is complete by piling the previous steps. 

4. The iV-PARTiCLE approximation at initial time 

4.1. Comparison of distances on probabilities. In the following lemma we compare 
the different metrics and norms defined Subsection 12.51 Let us denote 



Mk{f,g):=me,x{{f,{v)'); {g, {v)')} ■ 



Lemma 4.1. Let f,gG P(M ) and k G (0,oo), then the following estimates hold: 
(i) For any q £ (l,+oo) and any k £ [q — l,oo).' 



^-- .- ^^-.- ^Mi-^) 



(4.1) Wi{f,g) < Wgif,g) < 2— Mk+i{f , g)— Wi{f , g)^ 
(ii) For any s £ (0, 1], 

(4.2) 1/ - gU < 2(1-^) Ws{f,g) < 2'-'^'^ Wi{f,gr. 
(iii) For any s G {d/2,d/2 + 1), 

^^■'^> 11^ ^^^H~s^ (^2s-d)\A{d + 2-2s)J '-^ ^'i • 

(iv) For any s > and k > we have 



d 



k 



(4.4) [/ - g]l < C{d, s, k) Mk+iif, g)^^TWi) 1/ _ ^i^^+'^c^^^) 

for some constant C{d, s,k) > depending on d, s and k. 
(v) For any 

d ] d 



and k > we have 



s £ [ max "i 2 ' ^ r 2 "^ ^ 



(4.5) [f - g]*i < C{d,s,k) M,+,if,g)^mi^ \\f - gr^_'^^ 

for some constant C{d,s,k) > depending on d, s and k. 
Proof of Lemma \4-l\ Let us consider each inequality one by one. 

Point (i). The inequality (|4.ip is well-known in optimal transport theory, we refer for 
instance to [^21 [T?] . 

Point (ii). Let us prove inequality (|4.2|) . Let vr £ Il{f,g). We write 



ifio-mi 



< 

~«xll 



'ir{dv, dw) 



g-iv-S, _ ^-iw-e, 



< 2(1-") / \v-wY\iYT^{dv,dw), 
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which yields the first inequahty in (j4.2p by taking the supremum in ^ G M"^ and the infimuni 
in TT G Il{f,g). The second inequality is then immediate from the concavity estimate 

W,if,g)<(Wi{f,g)r. 

Point (hi). Let us prove the inequality (I4.3p . Consider R> and the ball 

Bii ■= ixGM.'^ ; \x\ <r\, 
and write 

ll/-9ll^-. = 



< \f-9\l[ 
JE 



2s di + 1^^ ^j^^ d^ 

R 



B'i, \i?' 



+ 4 ' 



|^|2(s-l) 



< 



5d-l 



R 



Br,,, .^^ 

i?2 



BZ le- 



d-2s 



\f-9\i + 



We optimize this estimate by choosing 

'4{d + 2-2s) 



R 



(2s - d) 



{d + 2-2s)' " {2s -d) 

\f-9\i' 



which yields 



||;_,||?,_^,it!lf,j£^iL,V-* ,;_#-. 



'^'^ - {2s-d) \4:{d + 2-2s)J 
which concludes the proof of (j4.3p . 

Point (iv). Let us now prove inequality (j4.4p . We introduce a truncation function 

where 

X G C,°°(M^), Ml < 1, and < X < 1, X = 1 on B{0, 1) 
and a mollifier 



7£(x) = e '^Tf")' e>0 with 7(x) 



e 2 



.£/ ■ ■ (27r)'^/2' 

In particular we have an explicit formula for the Fourier transform of this mollifier: 

%{0 =7(e0 =e = 2 . 
Fix (p G W^'°°{'R'^) such that [ip]i < 1, c^(0) = and define 

•^R ■= ^ XR, ^R,e = V'R*7e 
and write 

ip {df -dg) = j ipR^s {df - dg)+ j {^r - ipa^e) {df - dg) + j {p - ipn) {df - dg). 
For the last term, we have 



(4.6) 



{^R - V>) {df - dg) 



< j{l-XR)W\{df + dg) 

< j M 

Jb'l 



^ Rk 



{df + dg) < 



Mk+i{f,g) 
Rk 
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Concerning the second term, we observe that 

\V^r{x)\ < X (I) + Hx)\ \V{xr){x)\ < X (I) + ^ 

so that 

VgG[l,oo], \\VipR\\L,<C{q,d,x)R~^, 
for some constant dependmg on q, d and x- Next, using that 






\V>R -V>R,£\\oo < \\Vip 



R\\oo 



7e(x) |x| dx 



e\\V^ 



R\\oo 



-f{x)\x\dx < C{q,d,x)£:, 



we find 

(4.7) 



{V>R - V>R,e) {df - dg) 



< C{q,d,x)£- 



Concerning the first term, using Parseval's identity, we have (the "hat" denotes as usual 
the Fourier transform) 

1 



VR,e if - g) 



< 



1 



'PR % if - g) di 
f-g 



l|Vv^«|| 



Li 



l?h 



< C{d,x)R''\f-g\se-'^''+'-'^ 

(4.8) < C R'^ e-'^'^+'-^^ \f - g\s. 

Gathering ([O]) . (fiTT]) and (gS]), we get 

'Mfc+i(/,5) 



2 l«|- 



\^r'e-' ^di 



[f-g]l<C{q,d,x) 



j^iV4^(M + , + ^^,-(^+^-i)|/_,|^). 



This yields (14. 4p by optimizing the parameters e and R with 



d+s-l 



R = Mk+i{f,g)^ \f - g\s '^' e^^ 



and then 



e = Mfc+i(/,5)d+M^+-) \f -g\ 



d+k{d+3) 



Point (v). Let us now prove inequality ()4.5p . 

Let us consider some smooth 99 such that [(p\i < 1, ip{0) = and let us perform the 
same decomposition as before: 



(p (df -dg) = / (pR^s (df -dg)+ I (ipR - ipR^e) {df - dg) + / (c/? - ipr) [df - dg). 
The first term is controlled by 

(/ - g) 



fR,e {df - dg) 



V>R,e lel' 



ICh 



< W'PrAh^ 11/ -5 



\H- 



with 



\^R,s\\hs 



l\c\'\rR\'\^\'^'-'^m'd^^ 



1/2 



< \\'^{^R)\\L^mr'%{Ohoo 

< C{s) \\V{^r)\\^, e'-' <C{d,s)Ri e-^-^-'^ 
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with the same arguments as above. 

The second term and the last term are controlled exactly as in ()4.6p and (|4.7p . which 
yields 

[/ -S-Ji < C* ( e+ j^ + R2e'- ' \\f - gWn-s j ■ 

This yields (|4.5p by optimizing the parameters e and R with 

R = Mk+i{f,g)^^ 11/ - g\Q'^' e'-& 
and then 



e = Mfc+i(/,<7)^+^||/-5ll^+_r. 



D 



4.2. Quantitative law of large numbers for measures. Let us recall and extend the 
definition of the functional W^(/) which was introduced in (j3.2p . For any function 

D:Pk{E)xPk{E)^R+, {f,g)^D{f,g) 

(where A; > is the index of a polynomial weight possibly required for the correct definition 
of D) such that 

D{f, g) =0 if and only \i f = g 
it is legitimate to define 

V/GP,(M'^), >V^(/):=(r^,I)(vrf,/)>= / D {^,^ j) f^^ {dV) 

{D should be thought as some function of a proper distance). 

For well chosen function D, the goal of the next lemma is to quantity the rate of 
convergence 

Wi^(/) ^"'^°°) in the case E = M.'^. 

Lemma 4.2. We have the following rates for the W function: 
(i) Let us consider 



|2 

Then for any s G {d/2, d/2 + 1) and N > 1 there holds 



V/,5GP2(ffi^), Z)i(/,<7):=||/-<?|,^-.. 



(4.9) v/ e P2{M''), wSAf) = I hv - /111-. <^(^) < 



c_ 

N 

for some constant C depending on the second moment of f . 
(ii) Let us consider 

yf,gGP2iM.'), D2{f,g):=\\f-gfH-s. 
Then for any s > d/2 and N > 1 there holds 

(4.10) V/GP2(M'^), WS,if)= [ hv-f\\l-sdf^''iV)<^ 

for some constant C depending on the second moment of f . 
(iii) Let us consider 

V/,5ePi(M^), Ds{f,g):=Wi{f,g). 

Then for any r/ > there exists k > 1 such that for any N > 1 

(4.11) V/GPfc(M'^), wS,{f)= [ Wi{fi^,f)df^''{V)< ^ 



1 



for some constant C depending on r] and the k-th moment of f. 
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(iv) Let us consider 

V/,5eP2(K'^), Di{f,g):={W2{f\g)f. 
Then for any r^ > there exists k > 2 such that for any N > 1 

(4.12) yfePkiM.''), WE,if)= [ {W2{fi^,f))'df'^''iV)< 

for some constant C depending on r] and the k-th moment of f. 



C 



_^max{d,2} + rj 



Remarks 4.3. (1) Estimate ()4.12p has to be compared with the fohowing classical es- 

timate (see e.g. [64j): for any and any A*" > 1 there holds 

(4.13) yfGPd+.iR"), W^2(/)<-^ 

where the constant C > depends on the {d + 5)-th moment of /. It is worth 
mentioning that our estimate (14.120 improves on (14.130 when d < 3 and k is large 
enough. 

Similarly, if we try to deduce from (14.130 some estimate on the rate for the Wi 
distance, by using a Holder inequality we get that for any A^ > 1 

V/GPd+5(K'^), W^,(/)<^ 

iVd+4 

where the constant C > depends on the {d+5)-th moment of /. For any d this is 
worst than (|4.11|) as soon as when the probability / G P^iM.'^) with k large enough. 
(2) When f,g ^ P{W^) are compactly supported, observe that the estimate (|4.5|) 
improves into 

Vs>l [f-9]l<C\\f-g\\f_^, 

for a constant C depending on s and on a common bound R of the support of / 
and g. 

If furthermore d = 1, we can take s = 1 in order to apply (14. 9p in the proof of 
(14.110 below and we obtain the "optimal rate" of convergence in the functional law 
of large numbers in Wasser stein distance Wi: 



C 



In higher dimension d > 2, the restriction s > d/2 means that we do not produce 
a better estimate than (|4.1ip by this line of argument. 
(3) As was kindly pointed out by M. Hauray, estimate (|4.12p should also be compared 
with some estimates in [53] where the related quantity 



Z^{f) := f W, K,^{yj df^^{V,)df^''{V,) 

is considered. When / S P(]R'^) has compact support and d > 3, they prove that 

Z^(/) < 



tN/ - - C* 



where the constant depends on the support of /. 
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Since for any f,g^ Pi(IR ) and for any ip G Lip]^(M ) we have 

jRdN jR<i(N+l) 



N 

/ ip{v) dg{v) - ^yZ 



^{v,)dr\v) 



ip{v){dg -df){v), 
we deduce by minimizing in ip that 

W,if,g)< f Wi{g,fj.^)df'^^iV), 

jRdN 

and therefore 

W^^(/)<^^(/). 

As a consequence, when / S P(M ) has compact support and d > 3 we obtain 
from this Une of argument the stronger estimate 

^wSf) < ^,- 

It is hkely that one could obtain similar estimates to (|4.11|) by tracking the formula 
for the constants in the results of [2S] and combining them with moment bounds 
and some interpolation. 

On the other hand, observe that our estimate (j4.1ip is almost optimal in the 
sense that we can not expect a better convergence rate than (|4.11|) with ?/ = 0, as 
it is stressed in [HHl Appendix]. 



Proof of Lem,m,a \4-^ We split the proof into two steps. 



Step 1. Let us prove ()4.9p (note that (I4.10p is then readily implied by (14. 9p ). 
Let us fix / G ^2(1^ )• We write in Fourier transform 



1 ^ 



which implies 



f4-f' 



Wr. if) = \ U12S d^ df'^ (V) 

TH E L_.^ ^s U^df--iv). 



Observe then that 

f (e— r«_/(^))d/(«^.) = o, j = l,...,d, 

jRd 

which implies that 

/ ^ rnh ^cie<^y) = o 

jR(N+l)d \^\^^ 
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as soon as ji / J2 , and 



-iti-g 



/(e) df{v) 



1 _ e— ? /(e) - e-« /(O + \fm'\ df{v) 
1-|/(0P- 



We deduce that 












II iij^- 


(/) 

s 


1 ^ /• 


e--rC _ /(e) 


2 

7t f'^^(rlV) 


|^|2. 


"^ / Iffll^J 






1 /■ e— «-/(e) 


2 




iv L. |ep« 


ai; j{av) 






1 /■ (i-i/(e)P) ,, 


Finally, denoting 




M2:= / (t;)2d/(7;) 


we have 




/(e) = i + i(/,^)-c + o(M2iep), 


and therefore 













l/(6P 



which impUes 



(1 + i (/,«)• e + o(M2 lei')) (i - z (/, ^) • e + o(M2 ie|2) 
i + o(M2iep), 



wV (/) 

!l II u — s 



N 



N 



(1-1/(01^ 



I5l<i 



lei^ 



■(ie + 



:i-i/(e)p 



I5l>i 



|{|- *, 



I5l<i 



M2 ,, f 1 ,A 



< 



c_ 

N 



from which (14.91) follows. 



Step 2. Let us now prove ()4.1ip . 
We use first (14. 5p in order to get 



W, 



jRdN 



< C (M,+i(/) + A4+i(^^))''+-^ ||/.{)^-/ll^ 

J^dN \ 

We then perform a Holder inequality with exponents 



N ^||2 \ d+2ks jf^N fY\ 



P 



d + 2ks 



P 



d + 2ks 
d + k{2s-l) 



and get 



W^i(/)<C 



R-^^ 



(Mfc+i(/) + M,+i (^{)^)) ^+'=(^-1) <^(y) X 



d+fc(23-l) 
d+2fca 



^diV 



l^-fB-sdr^'iv) 



d+2k3 
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Since 



we deduce by using (|4.9p that 

_/Yd+2fes 

where the constant C(/, A;) depends on the {k + l)-th moment of /. 
This ahows to conclude the proof of ()4.1ip since 

• if d = 1 we can take s = 1 in (14. 5p and then k large enough so that k/{d + 2ks) = 
2 + 1] with some r/ > as small as wanted, 

• if d > 2 we take s close to d/2 and then k large enough so that k/{d+2ks) = l/d+rj 
with some ry > as small as wanted. 

Then the estimate ()4.12p follows from ()4.1ip with the help of ()4.ip and a Holder in- 
equality. 

D 

4.3. Chaotic initial data w^ith prescribed energy and momentum. In many as- 
pects, the simplest A^-particle initial data is the sequence of tensorized initial data /® , 
N > 1, where / is a 1-particle distribution. This means perfect chaoticity. On the other 
hand it has a drawback: since in all applications we shall use pointwise bounds on the 
energy of the A^-particle system (and also sometimes pointwise higher moment bound as in 
(Al)-(ii)), this implies for this kind of initial data that / has to be compactly supported. 
There is another "natural" choice of initial data, by restricting to one of the subspaces left 
invariant by the dynamics: 



1 ^ 1 ^ 



5- = <(yGM-^ s.t. j;^}^h\'=£, j;^}^v, = M 

■ =1 i=l 



Without loss of generality we shall often set A4 = and iS = 1 in this formula in the 
sequel. 

The drawback is now that an initial data on S^ cannot be perfectly tensorized, and 
some additional chaoticity error is paid at initial time. However an advantage of this 
viewpoint is that it is simpler to study the asymptotic behavior of both the A^-particle 
and the limiting mean-field equation in this setting. Moreover it has historical value since 
this approach was introduced by Kac (see the discussion in [42, Section 5 "Distributions 
having Boltzmann's property"]), although in his case there was only one conservation 
law, namely the energy, and therefore S was replaced by S (yiV). We shall present 
some results on the construction of chaotic initial data on S , whose proofs are mostly 
extensions of the precise statements and estimates recently established in |12j on this issue 
in the setting of Kac on S^^^. We refer to the work in progress [16] where an extensive 
study and precise computations of rates shall be performed. 

Lemma 4.4. Consider an initial data 

/o G Pi (m^) n L'^ 
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which fulfills some moment conditions 

Mmg^ if) = if, mg^) < +00, Mmg.^ (/) = (/, mg^) < +00 

for some positive radially symmetric increasing weight functions mg^ and mg^, and let us 
denote 

[ \v\^foidv)=££{0,oo). 

Let us define a non- decreasing sequence (aAf)Af>i «s follows: 

• // /o has compact support 

(4.14) Supp/o C {?; G M^ \v\ < a| 

for some ^ > 0, then 

\fN>l, ON ■=mg,^{A). 

• If fo ^^^ ''T'Ot compact support, then (aAr)Ar>i can be chosen as any non- decreasing 
sequence such that 

lim ajy = +00 

Af-5-oo 
(note in particular that this sequence can grow as slow as wanted). 

Then there exists 

/o^GP(M^^), iV>l, 

such that 

(i) The sequence {f^)N>i is fo-chaotic. 



(ii) Its support satisfies 



Supp/o^ C 5^. 



/o 
(iii) It satisfies the following integral moment hound based on mi : 



ViV>0, (/o^M4^)<C(/o,me,) 



where the constant C > depends on M^^ . 
(iv) It satisfies the following "support moment bound": 

(v) It satisfies a uniform relative entropy bound 

^^(/o^l7^) < C, 

for some constant C > (see (|1.7|) for notations). 
(vi) If furthermore the Fisher information associated to /o is bounded, that is /(/o) < 00 
(see (jl.Sp for notations), then /q^ can be built in such a way that it satisfies a 
uniform relative Fisher information bound 



^(/o It ) •= 1^7 

iV J^N 



Vlog 



d/o^ ' 



d'-f 



N 



pN 



f^' < c, 



for some constant C > 0, where the gradient in this formula stands for the Rie- 
mannian gradient on the manifold S^ . 
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Proof of Lemma \4-4\ For the sake of simplicity, we assume with no loss of generality that 
the energy £ = 1 and that /o is centered. 

We aim at defining our initial data /q by conditioning the tensorized initial data f^ 

to 5^: 

iN 



f-(v)= fr 

with 



^" ■ \f^ (Vn' 



SN 



N 



Such a construction obviously satisfies (ii). 

It is proved similarly as in ^12j (see for instance Theorem 9 in this reference) that this 
conditioned measure is well-defined, and that it is /o-chaotic, which proves (i). 

Rem,ark 4.5. Among many interesting intermediate steps and other results, it is also proved 
in [12] the following estimate: assume for simplicity that d = 1 and that /o has energy 1, 
then the function 

^'- 7^(r) 
is asymptotically divergent except for r = yN, for which 



with 




S = J/ {v^-lfdf{v) 



(In fact this result was sketched by Kac ^42j but the proof is made more precise in [12] . 

This shows in particular that the sequence of chaotic initial data f^ , A^ > 1, as 
considered many times in the sequel, asymptotically concentrates on the energy sphere 
with the energy of /q. This manifestation of the central limit theorem explains why 
the construction of Kac (to condition to a given energy sphere) is very natural. It also 
enlightens why it is possible to expect the kind of uniform in time propagation of chaos 
results that we shall prove in the next sections for such chaotic initial data. 

Point (iii) is just a consequence of the chaoticity with the test function M^ (actually 
an easy truncation and passage to the limit proceedure is needed in full rigor). 
Concerning point (iv), first if /q is compactly supported ()4.14p we deduce that 

Supp/o^ C {y G M^^, M4^(y) < mg,{A)] 

and (iv) holds. 

In the non compactly supported case, for any increasing sequence (^fc)fc>i of positive 
reals (with Aq big enough for the following to be well-defined) we define 

^ /ol|f|<Afc 

^°''" fo{{\v\<Ak}y 
Using the previous we know that 



/, 



N 
0,k 



®N 



•IO,k 



(conditioning to the sphere) is /o^fc-chaotic. Conditions (ii) and (iii) will therefore be 
immediately satisfied. 
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We now want to choose a sequence kjy ^ oo such that (iv) is satisfied and at the same 
proving chaoticity towards Jq. It is clear that 

For any given sequence (a at) which tends to infinity, we define kisf in such a way that 
"T-aC^fciv) — '^N so that /cat — )■ c« when A^ — ;• oo. The chaos property is equivalent to the 
weak convergence of the 2-niarginal, which can be expressed in Wasserstein distance for 
instance: 

w, (te), , r ) < w, ((/^,), , /«,^) + w, (/g, r ) . 

The last term of the RHS converges to zero only depending on A; — >• 0, while the first term 
in the RHS converges to zero for fixed A; as A^ — )• from the previous part of the proof. 
Therefore, maybe at the price of a slower increasing sequence k^ we can have both the 
support moment condition (iv) and 



which shows the chaoticity and concludes the proof. 

For the proof of (v) and (vi) we refer to [Kj . D 

Remarks 4.6. (1) We note that if one only wants to get rid of the compact support 

requirement in /o (used for deriving the support bounds on f^ on the energy and 
mg^ ) , and not necessarily to prescribe a given energy, another strategy could have 
been to simply perform the cutoff in the end of the previous proof. In principle 
it could allow to get better information on the rate of convergence. However a 
drawback of this approach is that, in the absence of conditioning to an energy 
sphere, the bound on the support of the energy of f^ shall grow with A^. In our 
applications it induces a growth in N of the moment bounds that we prove along 
time on the A^-particle system. This growth should be matched by the decay of the 
scheme and a precise optimized balance could be searched for. We do not pursue 
this line of research. 
(2) Observe that the process of conditioning on the energy sphere is obviously com- 
patible with the equilibrium states in the following sense: if one denotes by 7 a 
centered gaussian equilibrium of the limiting equation with energy normalized to 
1, then 

7^(y) := [^^%. 
is the uniform measure on 5^, i.e. an equilibrium of the A^-particle system. 

Let us also state a refinement of the previous lemma which is needed for the applications. 



Lemma 4.7. We use the same setting and assumptions as in Lemma \4-4\ We consider 
some function ©a(a^) such that 

Va>0, @a{x) ^— > 0. 

Then the sequence (/q^), N > 1 of Lemma \4-4\ satisfies the more precise chaoticity 
estimate: 



(4.15) Ww, (vr^ (/o^) , /o) = / W, (^^, /o) df^{V) ^^^^ 0. 
as well as 

(4.16) G.. {Ww, (vr^ (/cf ) Jo))) ^'^^ 0. 



with 



a =max{aAr; (/cmgg)}. 
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Remarks 4.8. (1) Using Lemma 14.11 it would be immediate to extend the previous 

statement to the other weak measure distances we have discussed so far. 
(2) Some polynomial rates could be obtained by refining the techniques from |12] . see 

LMJ. 

Proof of Lemma \4.7\ Let us prove that (|4.15p holds. First, thanks to [66t Proposition 2.2] 
and the fact that the sequence /g^ constructed in Lemma H^ is /o-chaotic, we deduce that 

N rN 



7^'^f^' ^5f, in P[P 

(which means convergence when testing against functions in C{P{W^))). 

Next, thanks to [72] Theorem 7.12], (|4.15p boils down to prove the tightness estimate 

(4.17) hm sup / Wi{p,f^)dT^^f^{p) = Q. 

R^Qo TveN* JWi{pJo)>R 

Let us prove that it follows easily from the following bound 

(4.18) sup / {Wl{^l^Jo)f df^{V)<<^. 
Indeed (j4.18p implies that uniformly in A^ > 1 



/ 



W, ip, /o) d7r^/o^(p) = / W, (/.^, /o) d/o^(y) 

Wi{p,fo)>R Jv^EN s.t. Wi{p^,fQ)>R 



<i/^^(«--,w./„)r<;/o"(v)4 



i?— 5>oo „ 
> 



which concludes the proof of (I4.17|) . 

In order to show ()4.18p . we infer that from j72l Theorem 7.10] 

{W,{p^,fo)f < ||/z{:^-/o||',. <2||^i:^||^,,+2||/o||i,i 

< 2 (Mf(y))' + 2 ||/o||^,i <2M2^(F) + 2 ll/oll^i, 
which implies 

/ {w^ {p^jo))' df^{v) < 2 m\l,i + Hf^,M^) , 

JE^ ^ 

which, together with (ii) in Lemma 14.41 ends the proof of ()4.18p and then of (14.150 . 

Then it is an easy diagonal process exercise to build a sequence {a^) such that (|4.16p 
holds by using the assumption on the function Ba together with (j4.15p . D 

5. True Maxwell molecules 

5.1. The model. Let us consider E = W^, d > 2, and a A^-particle system undergoing 
space homogeneous random Boltzmann type collisions according to a collision kernel 

B = r{z)b{cos9) 

(see Subsection ll.ip . More precisely, given a pre-collisional system of velocity variables 

the stochastic process is: 

(i) for any i' ^ j' , draw a random time Ty(\v-,-v-,\) of collision accordingly to an 
exponential law of parameter T(\vii — Vji\), and then choose the collision time Ti 
and the colliding pair {vi,Vj) (which is a.s. well-defined) in such a way that 



T: =Triu„-„.\\ ■■= min Tj 



1 = ^r(|.,-.,|) := mm ir(K,-.^,|); 
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(ii) then draw a G S'^~^ according to the law b{cos9ij), where 

cosOij = a ■ {vj - Vi)/\vj - Vi\; 
(iii) the new state after colhsion at time Ti becomes 

(5.1) V*j = {vi,...,v*,...,v*,...,vn), 

where only velocities labelled i and j have changed, according to the rotation 

(5.2) ^^=^-^ + ^' ^.=^ ^- 

The associated Markov process 

{Vt)t>o on {rX 

is then built by iterating the above construction. 

After rescaling time t — )• t/N in order that the number of interactions is of order 0(1) 
on finite time interval (see [65]) we denote by /^ the law of Vt and S^ the associated 
semigroup. We recall the notation G^ and T^ respectively for the dual generator and 
dual semigroup, as in the previous abstract construction. 

The so-called Master equation on the law f^ is given in dual form by 

(5.3) 5,</,^^> = (/,^G^c^> 

with 

(5.4) {G^^) {V) = ^ Y. r (b. - ^,1) f 6(cos%) [<^:, - y,] da 

iV /K!d— 1 

l<i<j<N -^^ 

where ip*j = if {V*j) and if = ip{V) G Gb (m^^) . 

This collision process is invariant under velocities permutations and satisfies the micro- 
scopic conservations of momentum and energy at any collision time 

N N N N 

^.* = ^., and \V*\^ = Y\v*/ = Y\v,\^ = \V\^ 
j=i j=i j=i j=i 

As a consequence, for any symmetric initial law /q G Psym(lK ) the law f^ at later times 
is also a symmetric probability, and it conserves momentum and energy: 

Va = l,...,d, / if2''j,a] f['{dV)= [ {J2^j,a] f^idV), 

where {vj^a)i<a<d denote the components of Vj G M.'^, and 

(5.5) V0:M+^M+, / c^i\V\^)f,^{dV)= [ cl>i\V\^) f^ (dV) 

(equality between possibly infinite non- negative quantities). 

The (expected) limiting nonlinear homogeneous Boltzmann equation is defined by (II. ip . 



p.2p . (II. 3p . The equation generates a nonlinear semigroup 

S^'-ifo) := ft for any /o G P2 (m'^ 

where ^2(1^) denotes the space of probabilities with bounded second moment. 
Concerning the Cauchy theory for the limiting Boltzmann equation: 

• In the case (GMM) (Maxwell molecules with angular cutoff), see equation (jl.6 
in Subsection ll.lt we refer to [67] : 
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• In the case (tMM) (true Mawell molecules without angular cutoff), see equa- 
tion ()1.5p in Subsection I l.H we refer to [69] : 

• In the case (HS) of hard spheres, see equation ()1.4p in Subsection II. H we refer to 
[59] (Li theory) and [291 [321 HH] (-P2(K'^) theory). 

For these solutions, one has the conservation of momentum and energy 



Vt>0, / vft{dv)= vUdv), / \v\'ft{dv)= \v\'Udv). 

JRd JRd Jm.'i JRd 

Observe that the change of variable 

maps the domain 

6* e [-vr, 7r/2] n [7r/2, vr] in 6* G [-7r/2, 7r/2]. 

Therefore without restriction we can consider, for the limiting equation as well as the 
iV-particle system, kernel function h such that Supp5 C [0,1]. We still denote by h the 
symmetrized version of 6 by a slight abuse of notation. 

In this section we aim at considering the case of the Maxwell molecules kernel. We shall 
indeed make the following general assumption: 

' r = i, 6eL~([o,i)) 

(5.6) I r 

Va>0, Caib):= b{cos 9) {I - cos e)^^"" da < oo . 

Let us show that Maxwell molecules model ()1.5p enters this general framework. Indeed 
for any positive real function t/j and any given vector u G M we have 

d-2\ / ,;,/„„„ Q\ „;„d-2 



V'(n ■a)da = \S'^'^\ / V(cose) sin^^^^ OdO. 
S'*-! Jo 

Therefore the model (|1.5p satisfies (in dimension d = 3) 

b{z) ^ K {1 - z)-^/^ as z^l, 



which hence fulfills ()5.6p . This assumption also trivially includes the Grad's cutoff Maxwell 
molecules model (11.61). 



5.2. Statement of the results. Our main propagation of chaos estimate result for this 
model then states as follows: 

Theorem 5.1 (Maxwell molecules detailed chaos estimates). Let us consider a 1-particle 
initial distribution /q G P{M.'^). and a hierarchy of N -particle distributions 

rN _ oN ( r®N 



as well the 1-particle of the limiting semigroup 



ft = 5f^ (/o 



where we assume that the collision satisfies ()5.6p . 
Then for any rj > 0, for any ^ G N* and for any 

(f = ipi (g) ip2 f^ ■ ■ ■ f^ 9?£ G 7"*^^ ipi^ J^, i=l,...,£, 
where J- shall be specified below, we have: 

(i) Cases (GMM) and (tMM).' Consider a tensorized initial datum /q = /|^ for 
the N -particle system, and assume that /o has compact support, and take 



T := {if-.W ^M.; \\ip\\jr 



■■= I {l + \C\')mO\dC<oo]. 
Jr<* J 
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(5.7) 



Then we have 

\/N>2£, sup 



< C/n 



.2 m 



N 



+ 



(5f(/o^)-(Sr(/o))^^),V. 



^l_^ ll'/^llj-2(g)(L°°)<-2 +^ ||99||,yi,oog,(ioo)<?-i >Vvy2(/o) 



for some constant Crj > (possibly blowing up as r] ^ 0) depending only r], on the 
collision kernel, and on the size of the support and some moments on /q. 

We deduce the following rate of convergence as N goes infinity by using (I4.12p - 



sup 

t>0 



S^{fi')-{Sr{fo)f''),^) 



< ^'n^' II II 



]\lK(d,ri) 



with 



r i - r? tf d<2, 



K{d,rj) := 



V if d = 3, 



77TI ^/^>4. 



d+4 



The constant C' > may blow up when r/ — t- 0, and depend on b and on /o through 

the size of its support and its moments. 

Case (GMM) with optimal rate for finite time: On a finite time interval [0, T], 

the following variant is available: consider tensorized initial data /g = f® for 
the N -particle system and assume that /o has compact support, and take J- = H'^ 
with s > d/2 high enough. Then we have 



sup 

0<t<T 



< C 



5f(/o^)-(5f^(/o))^^),^) 



,2 ll<^l 






N 






By using (|4.9p this proves 



sup 

te[o,T] 



(5f(/o^)-(Sr(/o))^^),V. 



^ /)2 Cr),T II II 



(iii) 



with T = K^ , with the optimal rate of the law of large number. 

Cases (GMM) and (tMM) conditoned to the sphere: Finally consider T as in 

(i), some initial data 



/oGA 



nL" 



and the N -particle initial data the sequence (/o^)Af>i constructed in Lemma 4-4 
and \4-7\ by conditioning to the energy sphere. Then we have 



sup 

i>0 



C„ 



rj2W 



5f (/o^) - {srifo)) 



®N 



.95 



N 



+ 



ii,„ii I I) ii,„ii ia;jv I „n fN X \ 



which goes to zero as N goes infinity thanks to Lemma \4^ and hence proves the 
propagation of chaos, uniformly in time. 
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Remark 5.2. Observe that by using the point (iii) below and taking t — )■ oo we deduce 
that 

for any fixed £ > 1, where 7^ is the ^-marginal of the uniform measure on the energy 
sphere 7^, and 7 is the gaussian equilibrium with energy £. However in our proof we 
need to build a sequence of initial data f^ which is supported on S^ and chaotic. And 
the only proof we know makes use of the fact that 7 is 7-chaotic. Nevertheless, this 
line of argument is not really based on an explicit computation nor on variational or 
entropy optimization characterization of the equilibrium. It only relies on the action of 
the N-particle and limiting semigroups. As a consequence, it can be applied to other 
situations where neither the steady states for the A^-particle system nor the steady state 
for the mean-field equation are known explicitely and where a direct computation cannot 
be made (which is typical of open systems). We refer to the work in progress \57\ for a 
application to some dissipative Boltzmann equation related to granular gases. 

We now state the key Wasserstein version of the propagation of chaos estimate, which is 
valid for any number of marginals, although with a possibly worse (but still constructive) 
rate. Combined with previous results on the relaxation of the A^-particle system we also 
deduce some estimate of relaxation to equilibrium independent of N and, again, for any 
number of marginals. 

Theorem 5.3 (Maxwell molecules Wasserstein chaos). Under the same setting as in The- 
orem [^n[ either 

(a) with /o compactly supported and /q^ = f^ or 

(b) with /o G P4 n L°^ and f^ constructed by Lemma 4-4' 
we have 



(5.8) ViV> 1, VI <^< A^, sup ^ -^ '-^<a{N) 

for some a{N) — t- as N ^ 00 (in the case (a) one has moreover explicit power law rate 
estimates on a). 

Then in the case (b) we have 

(5.9) VAf>l, Vl<^<Af, Vt>0, ^ ^ * ' ^ ^^ ^^ < P{t) 

for some /3(t) —)• as t —)• 00, where 7 is the gaussian equilibrium with energy £ and 7^ 
is the uniform probability measure on S^ {£). 

In order to prove Theorem 15. H we have to establish the assumptions (A1)-(A2)-(A3)- 
(A4)-(A5) of Theorem 13.11 with T = 00. The application of the latter theorem then 
exactly yields Theorem 15.11 by following carefully each constant computed below. Then 
the proof of Theorem 15.31 will be done in Subsection 15.91 it is deduced from Theorem 15.11 
by using Lemma 14.11 together a result from [37j . 

5.3. Proof of (Al). When the collision kernel B is bounded the operator G is a linear 
bounded operator on C{Br) with Br := {V G M'^^; \V\ < R} for any R e (0,oo) with an 
operator norm independent of R. As a consequence, G^ is also well defined and bounded 
on 

C7\o(IK"^):={^GC'(R'^'^) s. t. ^^0 as |y| ^ ooj 

endowed with the norm 

\ML-^:=sup\^iV)\{V)-^ 
V 
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for any A; G M. It is also easy (and classical) to verify that G^ is dissipative in the sense 
that 

V(^ G CLkfl (k'"^) , VA > ||(A - G^) <A|L^^ < A UWl^^. 

From the Hille-Yosida theory we deduce that G^ generates a Markov type semigroup TJ^ 
on C£^q(M'^^) and we may also define S^ by duality as a semigroup on P]^{W^^). The 



nonlinear semigroup S^^^ is also well defined on Pfc(R'^), see for instance [69 



semigroup u^ 

For the true Maxwell molecules model, the operator G^ is not bounded anymore and 
some additional explanations are needed. The fastest way to argue is just to say this in 
that case B can be approximated by a sequence of bounded collision kernels 

Be := bsicose) with b^ G L°° and b^ /- b. 

We may then define the associated generator G^'^, the associated semigroups T^'^ on 
continuous functions and S^ '^ on probabilities and the nonlinear semigroup S^ '"^ on 
probabilities. We first write estimate (j5.7p for any fixes e > 0. Then since (1) the right- 
hand side term in ()5.7p does not depend on e > (as a consequence of the estimates 
established in the proof below) and (2) S^^'%fo) -^ 5^^(/o) weakly in P{R'^) and (3) 
S^'^ifo^) -^ ft^ weakly in P(M^'^), we can conclude that (j5.7p holds for the true Maxwell 
molecules model by letting e go to 0. 

Possible other direct arguments (without using approximations) could be (1) to establish 
and use Wasserstein Wi stability of the many-particle equation, or (2) use the following 
core C := Wj^^ and prove that ip G Wf^^ implies G'^tp G G^fl (this follows from an easy 
decomposition between singular and non-singular angles in the formula for G ). 

Hence the semigroups T^^ and 



are well defined on C*^^ q(M ). Moreover since for tp G -L^(5 ) we have 

_/v 
{G^^,ip)^,^^^^ = _^ ^ /■ f 6(cos%) [v^:, - ip]' da^^'idV) < 

it is easily seen by arguing similarly as above that they are Co-semigroups of contractions 
on this space L^(5^). 

Then it remains to prove bounds on the polynomial moments of the A'^-particle system. 
We shall prove the following more general lemma: 

Lemma 5.4. Consider the collision kernel 

B = \v — Vi.l'^ b{cos 9) with 7 = or 1 

and 6 > such that 



/ b{z) (1 - zfdz < +00. 
Jo 



This covers the three cases (HS), (tMM) and (GMM). 

Assume that the initial datum of the N -particle system satisfies: 

N 



Supp/o^ C {y G M^'^; M^{V) < 8^] where M^ = ^^1 






and 



(/o^, Mf ) < Co,fc < 00 where M^ = -Y. I^^l'' ^ ^ 2. 
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Then we have 

sup(//^,M^) < max{Co,fc; Ofc} 
t>o 

where a/^ G (0, oo) depends on k and £q. 

Proof of Lemma \5.4\ By using ()5.5p with the function of the energy 

(p{z) := lz>Neo 
and the assumption on f^ we deduce 

(5.10) Vt > 0, Supp//^ civ e R^'^; Mi^{V) < £o\ . 

Next, we write the differential equality on the k-th moment: 

with 

1 /■ r , , , .1 

da. 






\v* I'^ + b* I'' 



We then apply the so-called Povner's Lemma proved in |59[ Lemma 2.2] (valid for 
singular collision kernel as in our case) which implies 

for some constants Ci,C2 E (0, oo) depending only on k and h. 

By using the inequalities \vj^ — Vj^\ > \vj^\ — \vj^\ and \vj^ — Vj^\ > \vj^\ — \vj^\ in order 
to estimate the last term when 7 = 1, we then deduce 



\vn - VnVnvn^v,,) < C, [(1 + \v,,\'+^-') (1 + 



''32\ 



(1 + K p) (1 + K,|^+^-i)] - C2 (b,j'=+^ + 



^h 



fc+7^ 



for a constant C3 depending on Ci and C2. 
Using (symmetry hypothesis) that 

VA:>0, {f!',\v,\') = {f!',M^), 

and (|5.1Up we get 

I (//^, W\^) <2C, {f,^, (1 + M,V,) (1 + M,^)) -2C, {f,^M,^,) 

<2C3{l+£)(l + (ff, Ivil"^''-')) - 2 C2 (/,^, l^il^+A 



Using finally Holder's inequality 

we conclude that y{t) = {ft , \vi\ ) satisfies a differential inequality of the following kind 

with ^1 > 1 and 82 < 0i, and for some constants Ki, K2, K3 > 0, which concludes the 
proof of the lemma. D 

Lemma 15.41 proves (Al)-(i) with 

mg^{v) := {vf. 
Moreover we do not need (Al)-(ii) in the present case since we may take mg^ = 0. 
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5.4. Proof of (A2). Let us define 

Pg. ■■= {/ e n^''); if, me) < S and (/,mg,) < +00} 
endowed with the distance induced by | -12, as weU as 

BPg,,a-={f^Pg,; {f,mg,)<a}. 

Let us recall the following result proved in |671 [33l EH EH] . We briefly outline its proof 
for the the sake of completeness and most importantly because we shall need to modify it 
in order to adapt it to our purpose in the next sections. 

Lemma 5.5. For any /o,5o £ P2{^- ), the associated solutions ft and gt to the Boltzmann 
equation for Maxwell molecules satisfy 

(5.11) 



sup|/t-5tl2 < 1/0-5012 

i>0 



Moreover, there exists a G (0, 00) such that 

VaG[a,oo), S^^:BPg^^ 



BPgu 



Proof of Lemma 15.51 We only prove (j5.1ip and we refer to the above quoted references for 
the moment estimates in the statement of the Lemma (which is nothing but a variation 
on the arguments presented in the proof of Lemma I5.4p . 

We recall Bobylev's identity for maxwellian collision kernel (cf. [6]) 



with 



and 






da, 



f, G = g, F±=F(e±), G± = G(C^), I 



^+ = ^(e + ick), r = ^(e-ick). 



lei 



With the shorthand notation D = g — f, S = g + f, the following equation holds 



(5.12) 



dtD = Q{S,D) 



b(a-^ 



D+ S- D- S+ 
+ 



D 



da. 



2 2 

We perform the following cutoff decomposition on the angular collision kernel: 



b = bx + b^; with 
for some well-chosen 5{K). As in [69J observe that 



bK [o' ■ S,] da = K, bx = bl 



satisfies 



Rk{0 



vee 



^K 



a-i 



+ 



D- S+ 



D 



\e\>5(K) 



da 



\RKm<rk\C\^ where 



Tfc ^0 



and rx depends on moments of order 2 on (i and s (hence bounded by the energy) . 
Using that H^Hoo < 2, we deduce in distributional sense 



d \D\ ,, \D\ 
— - — - ^K - — - < 



sup 



\D\ 
W 



sup / bx (cr ■( 



e^ 



+ 



r 



da \ +rK 



with 



V2 



1 + a-e 



1/2 



V2 



1-a-^ 



A 1/2 
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+ 



d \D\ ^^ \D\ ^^ / \D\ \ 



By using 



we deduce 



which imphes 

for any value of the cutoff parameter K. Therefore by relaxing i^ — ?• oo, we deduce 

Hence we deduce that S'f^ e C°'^{Pgi, Pgi) and (A2)-(i) is proved. 
Lemma 5.6. For any f,g€ P(M'^) with finite second moment 



f \v\ dv < oo, 



g \v\ dv < oo 



and same momentum 

we have 

(5.13) 

and 



f Vi dv 



\QifJ)\2< 



gvidv, i = l,...,d, 



c(/^^(i + MV/w)^ 



(5.14) |(3(/ + 9,/-9)l2<C ( /^(l + kl) {dfiv)+dg{i.))\ [\f - gh + Hf - g)v\,). 



Proof of Lemma 15.51 We prove the second inequality (I5.14|) . The first inequalities (I5.13P 
then follows immediately by writing 

Qif, f) = Qif, f) - Q{M, M) = Q{f-MJ + M) 

where M is the maxwellian distribution with same momentum and energy as /, and then 
applying (fSTTil) with f - M and / + M. 
We write in Fourier: 



J'{Q{f + 9J-9)) = Q{D,S) 



1 



b{^ ■ i) {S{t) D{C) + S{C) D{t) - 2 D{0) 



where Q is the Fourier form the symmetrized collision operator Q, which we can rewrite 

Q{D,S) 



l?l^ 



<Ti + T2 + Ts. 



We have 



ri< 



--|2 



bia-0 \S{e)\ ^^^U ^47PrdcT<C\D\2 



^-1 ■ ■ le-P \e 

for some constant C > 0, where we have used 

|2 

\2 



\1_ 



(l-cos6l)^ 



which permits to control the angular singularity of b. 
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Similarly we compute 

T< f h(n ,A i-p(e+)i |gr)-2i in , .^ipi 

Js<*-i 14^1 I? I l?l 



[l + \v\){df{v)+dg{v)) 



for some constant C > 0, and 
T3<2 



6(..|)M:tM^.. 



lei 



for some constant C > 0. This concludes the proof of (j5.14p by piling these estimates. D 

The proof of (A2)-(ii) is a consequence of (|5.14|) . Indeed, from (|4.2|) we have 
yf,geBPg^,, \{f-g)v\,<C[{f-g)v]l. 
Moreover for f,g £ BPg^^a, 

[if-g)v]l < M{[(f-g)vxn]l + [{f-g)vil-XR)]l} 

< Ca'/'ilf-gTif 

where we have used the bound on the fourth moment. 
Finally, from (|4.4|) . we conclude that 

yf,g€BPg,,a, [f-g]l<C\f-gr^, 

for some constants C > and i^ £ (0, 1) depending on d and a. 
Gathering all these estimates as well as (|5.14p we deduce that 



\Q{f + 9j-9)\, 



< C 



{l + \v\) {df{v) + dg{v))^ (^|/_5|2 + |/-5|^). 



5.5. Proof of (A3). Let us define m'g and 

Ai(/) :=(/,<,) = (/,(^>'> 

for any f e Pg^. 

Let us prove that for any 



i-*?, 



rGRg, and $ S C^f (Pg, 



we have 



^4(^)) (G^vr^-^^G-)$ 



L°°(£jv) 



Ci £" I- I 



for some constant Ci > 0. 

First, consider velocities v,v^,w,w^, G M'^ such that 



V+V^ \V — vj V + V^, \V — Vj ^j_-| 

w = — 1 a, w^: = — a, cr E s> 



KAC'S PROGRAM IN KINETIC THEORY 



55 



Then 5y + 5^^ — Sw — 5w, G ^^-Pgi • Performing Taylor expansions, we get 



e™-? _|_ e™*'^ — e™'^ — e™*''^ 



+i K - V,) e (e^^'« + O (It; - t;*| lel) + O {\w, - v,\^ \^\^) 
= 0{\v-v,\'^\C\'^ sme/2) 

thanks to the impulsion conservation and the fact that 



\w — v\ = \w^ — f*| = \v — v^\ sm -. 



We hence deduce 



Ig™-? _|_ £«■"*■? _ giw-S, _ gJi«*-C| 
\Sv + Sv, - Sw - 5wJ2 = S^P "^ |7T2 ^ < Clf - V*|^ (1 - COS0). 



As an immediate consequence, for any V £ E and V*; defined by (j5.2p . we have 



Hy. fly 



c 

2 - N 



< ^ \Vi -Vj\^ (1 -COS %) 



Consider V G E^r and define 

rv := ((A^y,2i),---,(/^{/,^d),(i^y,kP)) G ^S 
Then for a given <I> G C'A"'(^Gi.rv)' ^^ ^^^ 

4> := D$ [/xy] and Uij = {vi - Vj) 
and we compute: 



1 



N 



TV 



cDU^:^, -cD(^^) 6(cose,j)dcT 



^1.7 — ^ 



2N 






o 






l+r; 



da 



=:h{V) + l2{V). 
Concerning the first term /i(I^), thanks to Lemma 12.131 we have 

1 ^ r 



6(cos 0) [(t>{v*) + ^{w*) - (t>{v) - 4>{w)] Hy{dv) fiyidw) da 
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For the second term l2{V), using that 

1 ^ 



v^^ij 



2 



3 / ^ 



^ ^ (^1 + 77 (^^'"^''jj -^^<i K) =CM^,^(y), 



we deduce 






We finally use that 



y oil b {cos Oij) (1 + \vi\^ + \vj\^) (1 



^M^^^(F)oQ^_^6(cos%) (l + |^;,|2 + |^^.|2) (i-a.n,^.)d^) 
(recall for the last line that mg^iv) = (v)^ and m'g (v) = (v)'^), which implies 



cm:^^(v) 



mg 



and concludes the proof. 

5.6. Proof of (A4) uniformly in time. Let us consider some 1-particle initial data 
/o)5o £ -P4(M'^) (probability with bounded fourth moment). And let us define the associ- 
ated solutions ft and gt to the nonlinear Boltzmann equation (|1.1|) under the assumption 
(|5.6|) . We then define 

ht := £r [/o] (50 - /o) 
the solution to the linearized Boltzmann equation around ft, given by 

dtft = Qift, ft), f\t=o = h 

dtgt = Q{gt,gt), g\t=o = 9o 

, dtht = 2 Q{ht, ft), h\t=o = ^o := 5o - /o- 

We shall prove assumption (A4) with the choice of indices 

which means that the weight which has to be used is 

We shall now expand the limiting nonlinear semigroup in terms of the initial data, 
around /q. 
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Lemma 5.7. There exists A G (0,oo) and, for any r] G (0,1), there exists C^ > such 
that for any 

r E Rg^ and fo,go^Pgi,r 
we have 

(5.15) |/t-9tl2<C,e-(i^^)^* (max{M4(/o),M4(5o)})^ l/o-ffoll, 
and 

(5.16) |/ii|2 < C^e-(i-'')^* (max{M4(/o),M4(9o)})^ l/o-5oll 

where we recall that 

V/GP(M'^), M4(/):=(/,H^>. 

^s a consequence, the operator C^^[fo]{-) is bouned from Qi^r to itself. 

Remark 5.8. Observe that the loss of weight exactly matches the need of assumption (A4) 
since 

As = Ap 

(recall that we have made the choice the indices {ri',rj") = (??, 1) in (A4) in the present 
application). 



Proof of Lemma 5.7. We shall proceed in several steps. 

Step 1. Estimate in \ ■ {4. We closely follow ideas in [671 E]- We shall use the notation 

M = Ma, M = Mi, 
introduced in Example 12.5.41 as well as 

d:=f-g, s:=f + g 

and 

d:=d-M[d], D := F{d), S := F{s) and D := F{d) = D - M[d\. 

The equation satisfied by D is 

(5.17) dtD = Q{D,S)-dtM[d\ 

= Q{b,S) + (Q{M[d],S)-M[Q{d,s)]]. 



We infer from [67] that for any a G N , there exists some absolute coefficients (a^^^), /3 < a 
(which means /3j < Oj for any 1 < i < d), depending on the collision kernel b through 

(5.18) / ^ 6(cos0) [iv^y + {v% - (v^) - (z;")J da = Y^ a,,^ (v^) (v^-^^ 

where a, (3 G N are coordinates indices and 

V := Vi V2 . . .V^ . 

These multi-indices are compared through the usual lexicographical order, and we use the 

standart notation 

d 

\a\ := y^Qfc- 

fc=i 

We deduce that 

V|q| < 3, V'^M[Q{d,s)M^^^ = M^[Q{d,s)] = ^ a^,^ M^[d]M^_p[s] 

/3,/3<a 
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together with 

V|a| <3, V'^Q{M[d],S)^^^^ = M^[Q{M[d],s)] 

since 

M^[M[d]] =M^[d] 
for any |a| < 3 by construction. As a consequence, we get 



Y, a^^pMp[d]M^_p[s] 

/3,/3<a 



(5.19) \J^eM.\ M[Q{d,s)]-Q{M[d],S) <C\i\^ [Y.\^'^^f^-9t\\ 

\H<3 

On the other hand, from |67l Theorem 8.1] and its corohary, we know that there exists 
some constants C, A G (0, oo) such that 



(5.20) Vt>0, iY. Wa[ft-gt]\\ <Ce-^M Y l^"[/o-5o] 

\l"l<3 / \|«|<3 

We perform the same decomposition on the angular colhsion kernel 



b = bK + bK with / bK [cr ■ n da = K, bx = b he\>5(K) 
as in the proof of Lemma 15.51 and use the straightforward estimate 

RKio--=Qbi.{D,sm 

satisfies 



X->oo 



>0 



VeeM^ |i?i^(e)|<rj^|e|' where r^ 
where Qb<= denotes the collision operator associated with the part 6^ of the decomposition 
of the angular collision kernel, and where r^ depends on moments of order 4 on d and s. 
Then we gather Jjf^TTt} . (f^T^ and (f^:^ and we have 



^m)\^K\^J,^,m)\ 



dt leh 



\^\' 



\^\' 



sup / bx (a -^ 



^' 



+ 



r 



da 



+ Ce-^* E l^^[/o-9oJ| I +rK. 
\H<3 

Let us compute (the supremum has been droped thanks to the spherical invariance) 



^K ■= I bx [a -S, 



C' 



+ 



r 



da = bx [a ■ $, 



l+ia-C 



A 2' 



da, 



so that 



Xk-K 



Igd-i 



bK (cr -^ 



l-ia-C 



K^oo J gd 






da 



1- a-f 
i ) I ^7^ — I da := -A G (-oo, 0) 



A 2' 
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where in the last step we used the factor 



■ ~ C 6|2 as 6* ~ 
2 I 

V 



in order to control the singularity of b. 
Then, thanks to Gronwall lemma, we get 



sup — -—7 — S e^ '^ ' sup 






g -At ^(XK—K)t > ^ 



J«l<3 



Therefore, passing to the limit K — )• 00 and choosing (without restriction) A G (0, A), 
we obtain 

or equivalently (and with the notations of Example 12. 5. 4p . 

|||dt|||4<Ce-^*|||do|||4. 



Step 2. From | • I4 to | • I2 on the difference. Prom the preceding step and a straightforward 
interpolation argument, we have 

\f-9\2<\f-9-M[f-g]\, + C (^ \MM-9]\ 

\hl<3 



<\\f-g-M[f-g]\\'j;^\f-g-M[f-gf/' + C Y. l^aif - g]\ 

< C {1 + AUifo) + M^igo)) e-(^/2)*. 
Then by writing 

\f-9\2<\f-9\l\f-9\l~\ 

using Lemma fS.lll for the first term of the right hand side and the previous decay estimates 
for the second term, we obtain 

\ft-9t\2<Cr,e-^'-^^^' {max{M,{fo),M,{go)}f-^^ |/o-5o|^ 
This concludes the proof of (15.150 by using (1 — 77) <!/(! + ry). 

Step 3. From the difference to the hnearized semigroup. The same computations imply 
exactly the same estimate on ht as on the difference {ft — gt), that is inequality ()5.16p in 
Lemma 15.71 D 

We can now consider the second-order term in the expansion of the semigroup. Let us 
recall that the crucial point here is to prove that this second-order term is controlled in 
terms of some power strictly greater than 1 of the initial difference. 
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Lemma 5.9. There exists X G (0,oo) and, for any r/ G (0,1), there exists C^ such that 
for any 

r G Rg^ and fo,go^Pg^,r, 

we have 

li^iU < C e ^ 1 50 - /0I2 

where 

cot ■■= gt-ft~ht = s^Hgo) - srUo) - /:r[/o](5o - /o). 

Remark 5.10. As proved below ujt always has vanishing moments up to order 3, which 
implies that the norm \ujt\4, is well-defined. Moreover observe that in this estimate there 
is no loss of weight. 

Proof of Lemma \5.9[ . We consider the angular cutoff decomposition as in Lemma [5.5[ Con- 
sider the error term 

to := g — f — h, [7 := d). 
which satisfies the evolution equation 

dtUt = Q{^t,f + g)-Q^{hJ-g), Wo = 

and (in the Fourier side) 

dtn = Q{n,s)-Q+{H,D). 

Let us prove that 

V|a|<3, Vt>0, Ma[oJt]-= v"du;t{v) = 0. 

We shall use again the fact that, for maxwell molecules, the a-th moment of Q{fi,f2) 
is a sum of terms given by product of moments of /i and /2 whose orders sum to \a\, see 
equation (jS.lSp . 

We obtain 

V \a\ < 3, j^Ma[uJt] = Y^ aa,/3 Mp[u;t] M^_p[ft + gt] + Yl «°-/3 ^/^N M^~p[ft - gt] 

and since 

V|a|<l, M^[ht] = M^lft - gt] = 0, 
we deduce 

V |a| < 3, -^MaiuJt] = Y^ aa,i3 Mfs[iOt] M^-plft + gt]- 

This concludes the proof of the claim about the moments of ut since ujq = 0. 
We now consider the equation in Fourier form 

dtn = Q{n,s)-Q+{H,D) 

and we deduce in distributional sense 
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(depending on some moments of order 1 of oj, h, d), and 



Ti 



sup 



b(a-^ 



\C\' 



m^)s{c 



+ 



m-)si^+) 



da 



< sup 



< 



b{a-C 



mm ic+i' , Mm \rf 



\t 



sup 






sup 



b(a-i 



+ 



\t 



ler 



+ 



da 



da 






where Xk was defined in Lemma 15.51 and 
1 r b(a-i 



T2 



77 sup 



1 



lel^ 



< - sup / 6(cr-C 



^(e+)i i^(r)i irp , mm? mm? le 



le+p le-p \i? 

6 ( a • Co ) f 1 - (7 • lo) da 



+ 



le 



\i' 



\i? 



da 



< \ht\r, \dt\ 



= Ce-^i-")^* I/10I2 |do|^ < Ce^^i-")^* |dol2^^ 

by using the estimates of Lemma 15.51 
Hence we obtain 



\m[\ 



mA 



l^^-^jS>.(^^S^j^C.-.-.,oir'.^.. 



We then deduce from the Gronwah inequahty, relaxing the cutoff parameter K as in 
Lemma 15.71 and assuming without restriction (1 — i])X < A, that 



sup 






<Ce-(i-^)^* 190-/01^+". 



This concludes the proof. 



D 



5.7. Proof of (A5) uniformly in time in Wasserstein distance. We know from ^67j 
that 

sup W2 (5f ^/o, 5f ^50) < W2 (/o, 50) • 



t>0 



As a consequence, by using 

Hi = Wi< W2, 

we deduce that (A5) holds with 

G(x) = x, J3 = Lip(M'^) and Pg^^ = P2{R'^) 

endowed with the distance dg^ = W2- 

By using Theorem 13.11 whose assumptions have been proved above, this proves point (i) 
in Theorem 1 5 . 1 1 and the rate follows from the estimate on W!t2{f) from Lemma |4.2[ 

By using Lemma WTH in order to relate yVw'^i'^pifo^)^ fo) with WvKi(7rp (/(^), /o) and 
then Lemma 14.71 in order to estimate 



7V-5>oo 



y^wA^pifo)Jo)^^^^o 
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for the sequence of initial data conditioned on the energy sphere constructed in Lemma 
we then deduce point (iii) in Theorem 15.11 

5.8. Proof of (A5) with time grooving bounds in negative Sobolev norms. It is 

also possible (and in fact easier) to prove, in the cutoff case, that the weak stability holds 
in Sobolev space on finite time: 

Lemma 5.11. For any T > and s > d/2 there exists Ct,s such that for any ft, gt 

solutions of the Boltzmann equation for Maxwell molecules ()5.6p and initial data /q and 
go, there holds 

snpWft- gt\\H-^ <Ct,s Wfo-goWH-"- 

[0,T] 

Sketch the proof of Lemma \5.11[ We integrate (|5.12p against D/{1 + |Cp)*: 
d„ „n Iff /A \D-S+D + D+S- D- 2101"^] 

and we use Young's inequality together with the bounds 

\\S^\\ . Il^^ll <II/ + 5|Imi<2 
to conclude. D 

This proves (A5) with the alternate choice 

e(x) = X, J"3 = H'iR'^) and Pg^ = P2{R'^) 

endowed with the distance of the normed space Q^ = i7*(M'^). Then point (ii) in Theo- 
rem [5T] follows from the abstract theorem 13.11 where the rate is provided by the estimate 
on W/^_s-,2(/) from Lemma 14.21 

5.9. Proof of infinite-dimensional Wasserstein chaos. We shall prove Theorem 15.31 

in this subsection. Let us proceed in several steps. Let us emphasize that we do not search 
for optimality on the rate functions given by our argument. 

Step 1: Finite-dimensional Wasserstein chaos. It is immediate that Theorem 15.11 
implies that, under one of the two possible assumptions on the initial data, for any given 
i > 1, one has 



sup 

t>0 



H-" 



for some power law rate function ao{i, N) — t- as A^ — t- 0. 
Then by using Lemma 14.11 we deduce that 

supW,(Ue[ff],f?')<aii,N) 

for some power law rate function a(^, A^) — )• as A^ — ?• 0. 

Note carefully that at this point our rate function still depends on I and in fact a quick 
look at Theorem 15.11 shows that they scale like (^ , therefore making impossible to choose 

Step 2: Infinite-dimensional Wasserstein chaos. We shall use here the following 

result obtained in ^STJ, see also [56, Theoreme 2.1]: for any / G ^(M ) and sequence 
/^ G Ps.vm(M'^) we have 

for some constructive constant C, ai, 02 > 0. 
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By combining this estimate with the previous step we immediately obtain 

sup supVFi(n4/,^],/f^)<a(iV) 

l<i<N t>0 ^ ' 

for some power law rate function a(A^) — t- as A^ — t- 0. This concludes the proof of ()5.8p . 

Step 3: Relaxation in Wasserstein distance. We shall prove (15. 9p and we shall 
consider here initial data j^ constructed by conditioning /® to the energy sphere. We 
first write 

w,^!^,^^) ^ "^^[fi^fr) ^i(/r,7^^) 1^1(7^^,7^) 

N - N ^ N '^ N ' 

Since /j — )■ 7 in L^ and /j — t- 7 in L^ as t ^ +00, one can pass to the limit in the 
Wasserstein distance and get from the previous step 

N 
Moreover it is immediate that 



< a{N). 



Wiift,7). 



N 
Finally it was proved in [HU] that under our assumptions on /g one has 

\\{ft-l){v)h.<Ce-'-' 

for some constants C > and Ai > which implies 

W,{ftn)<\\{ft-l){v)\\L^<Ce-^^'. 

Hence, gathering these three estimates, we deduce that 

(5.21) Wi{fj^,j^) ^2a{N) + Ce-^''. 
It was proved in [T5] that there exists A2 > such that 

where h = df /d'y is the Radon-Nikodym derivative of / with respect to the measure 
7^ so that /^ = h^-i^. When f^ = [f^^]s,M-.^^-) with /o G Pi{W^) we easily upper 
bound the right hand side term by 

\\h^ - ill < 4^ 

II \\l'^{S^ ,'y'^) — ' 

where A = A(/o) > 1- Thanks to Cauchy- Schwartz inequality and the control of the 
Wasserstein distance in terms of a weighted total variation norm (see [72^ Proposition 
7.10]) we also have 

II"- ^llL2(§diV-i(^^^^]V-) ^ II"- ^llLi(5'«'V-i(y7v),7^) 

1=1 

and we deduce 

(5.22) ViV>l, Vt>0, VFi(/f,7^) <^^e-^2t_ 
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Finally by combining (|5.2ip and (|5.22p and taking the better of these two estimates 
depending on N and t, we easily obtain 

ViV>l, vt>o, Wi{ft^,^'')<m 

for some /3(t) — )• as t — )• +00, which concludes the proof of (j5.9p . 

6. Hard spheres 

6.1. The model. The limiting equation was introduced in Subsection 11.11 and the sto- 
chastic model has been already discussed Subsection 15.11 

We consider here the case of the Master equation ()5.3p . (15. 4p and the limit nonlinear 
homogeneous Boltzmann equation (II. ip . ()1.2p . ()1.3p with 

(6.1) B{z,cos6) = r{z)b{cose) = r{z) = \z\. 

6.2. Statement of the result. Our fluctuations estimate result for this model then states 
as follows: 

Theorem 6.1 (Hard spheres detailed chaos estimates). Let us consider a 1 -particle initial 
distribution /o G P(M'^) with energy less than £ 

M2(/o) < £, 
and a hierarchy of N -particle distributions 

j-iV _ oN ( r^N 

issued from the initial data /g , as well the 1-particle of the limiting semigroup 

ft = 5f^ (/o) 

where we assum,e that the collision satisfies (j6.ip . 
Let us finally fix some 6 G (0, 1). 
Then 
(i) Suppose that /q has compact support 

Suppfo C jw eR'^, \v\ <a\ 

and that the N -particle initial data is tensorized 

viv>i, f^ = fr- 

Then there are 

— some constants ki >2 depending on 5 and £; 

— some constant Cs > depending on 5 and £ , and blowing up as 6 ^ 1; 

— some constant Cf, > depending on the collision kernel, 
such that for any i gW , and for any 

(p = ipi^iP2Ci)---(^ fee wi'°°(M'^)®^ 

we have 

(6.2) yN>2i, sup ((5f(/o^)-(5f^(/o))n,^ 

[0,T] ^ ^ ^ 



< ||93||vi/i,°°(IRd)S 



2f Csnlfohn 



The last term of the right hand side (which is also the dominant error term as 
N goes to infinity in our estimate) is given by 

C 

^^^^ " (l + |logA^|)" 
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for some constants C, a > 0. 
(ii) Under the same setting but assuming instead for the limiting initial data 

/o G L°° (rA s. t. [ e^l^'l dfo{v) < +00 

for some z > 0, and taking for the N-particle initial data the sequence (/q )Ar>i 
constructed in Lemma \4^\ and \4-'T\ we have the same estimate ()6.2p however with 
no information on the rate 

e{N) ^^ 0. 

This nevertheless proves the propagation of chaos, uniformly in time. 

We now state the key Wasserstein version of the propagation of chaos estimate, which is 
vahd for any number of marginals, although with a possibly worse (but still constructive) 
rate. Combined with previous results on the relaxation of the A^-particle system we also 
deduce some estimate of relaxation to equilibrium independent of N and, again, for any 
number of marginals. 

Theorem 6.2 (Hard spheres Wasserstein chaos). Under the same setting as in Theo- 
rem\6.1i either 



(a) with /o compactly supported and f^ = /® or 

(b) with 



/o G L~ (rA s. t. I e^l^l dfo{v) < +00 



and /q constructed by Lemma \4-4 
we have 



w, (n,//^, (ft 

(6.3) ViV> 1, VI <^< A^, sup ^ -^ '-^<a{N) 

t>o " 

for some a{N) -^ as N ^ 00 (in the case (a) one has moreover explicit estimates on 
the rate a like a power of a logarithm). 
Then in case (b) we have 

, , Tyi(n£/f,n^(7^)) 

(6.4) VAf > 1, VI <^< Af, Vt> 0, ^ ^ ^ '' </3{t) 

for some /3(t) — )• as t — )• 00, where 7 is the gaussian equilibrium with energy £ and 7^ 
is the uniform probability measure on S {£). 



In order to prove Theorem 16 -H we shall prove assumptions (A1)-(A2)-(A3)-(A4)- 
(A5) of Theorem 13.11 with T = 00. The application of the latter theorem then exactly 
yields Theorem 16.11 by following carefully each constant computed below. We fix 

Ti=T2 = C6(M°') and T3 = Lip(R'^). 

Then the proof of Theorem 16.21 is deduced from Theorem 16.11 in a similar way as Theo- 
rem 15.31 was deduced from Theorem 15.11 see Subsection 16. 9[ 
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6.3. Proof of (Al). From the discussion made in section [5^ we easily see that for the 
Hard sphere model the operator G^ is bounded from C-k+i,oi^'^'^) onto C^kfli^'^^) for 
any fc G M. Since G^ is close, dissipative and C^k+iflO^'^'^) is dense in C_fc,o(IK'^^), the 
Hille-Yosida theory implies that G^ generates a Markov type semigroup T^ on C-kflO^'^^) 
and we may also define S^^ by duality as a semigroup on Pfc(R'^^). The nonlinear semigroup 
S'f^ is also well defined on Pk(^'^), k>2, see for instance [32| [29| liH] . 

Lemma 15.41 was proved both for Maxwell molecules and hard spheres. It first shows 
that 

\ft > 0, Supp//^ C En := {v G R'^^; Mi' {V) < sl . 

It also proves that for any k > 2, 

sup(//^,Mf><Cf 
t>o 

where C^ depends on k, £, on the collision kernel and on the initial value 

(/o^Mf) 

which is uniformly bounded in A^ in terms of k and A. This shows that (Al)-(i) holds 
with rni^v) := \v\^'^ for any ki > 2. The precise value of ki shall be chosen in Section f6. 71 
As for (Al)-(ii), we remark that for a given A^-particle velocity 



we have 



V E Supp/o®^ ^^ Vi = 1, . . . ,iV, t;^ e Supp/o 



which implies 

\/i = l,...,N, mg^{vi) < mg^{A) with mg^{v) := e''^''^ 

for any constant a > 0, which shall chosen later on. 
We conclude that 

Supp/o«^ C {y G M'^^; M^JF) < mg.M)} , 
and (Al)-(ii) holds for the exponential growing weight mg.^. 
6.4. Proof of (A2). We define 

Pg, ■■= {f € PiR''); {f,{v)''^)<+oo] 

that we endow with the total variation norm || • \\g-^ := || • ||j\//i. Note that since ki > 2, 
elements of Pg-^ have in particular finite energy. 

This assertion (A2)-(i) reads: 

For some constant a^j > and for any a > a^j and for any i > 0, the 



application 
maps 



/o ^ Sf^/o 



BPgua ■■= {/ e Pg„ MkAf) = (/, {v)"') < a} 
continuously into itself. 
It is postponed to section 16.61 where we prove in (|6.6|) a Holder continuity of the flow in 
iSPgua- 

Let us consider assertion (A2)-(ii), that is the fact that the operator Q is bounded and 
Holder continuous from Pg^ to Qi. 
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For any f,g£ Pg^ we have 

\\Q{9,9)-Q{fJ)\\M^ = \\Q{9-f,9 + f)hn 

< 2/ / / b{9)\v — v^:\\f — g\\f^: + g^:\dadv^:dv 

i^ + £)\\b\\L^\\if-9){v)Un. 

We deduce that 

\\Q{9, 9) - Qif, f)hn < 2 (1 + £f' \\bh. 11/ - g\C 
which yields 

and also implies that Q is bounded on Pg-^ since we can choose 5 to be a maxwellian 
distribution. 

6.5. Proof of (A3). For any ki > 2, any energy 

rd+i > 

and any mean velocity 

(.ri,...,rd) G B^d {0,^/rZn) 
we define 

PGi,r ■= {f ^ Pgi ; {f,Vj) =rj, j = l,...,d, (/, Ivp) =rd+i} 

with the notation 

r := (ri,...,rrf+i). 

This corresponds to the vector of constraints 

m = {vi,...,Vd,\v\'^) 
in Definition 12.41 

We also denote by R the set of all admissible contraint vectors r G M^^: 

R = J (n, . . . , ra+i) G M'^+^ s. t. ra+i > and ^ rf < ra+i \ . 

Let us finally define the weight 

Ai(/):=M,,_2(/) = (/, (^')'^-'> 

(this means that we choose m'g (v) = {v)^^^'^ in assumption (A3)). 

We claim that there exists a constant Ck^ > (depending on ki) such that for any 
rj G (0, 1) and any function 

cI>G fl Clf(Pg„r;M), 
reRf 

we have 

(6.5) yVGR'"'', \G''(^o^,^)-iG^^){^^)\<C,Jsv^mc^^^^p^^^^ 
where we recall that 

TV 



K(y)--=lT.(-^)''- 



N 

i=l 



This would prove assumption (A3) with the rate 
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Consider V G M'^^ and define 
Then let us consider a given 



and let us set 



We compute 



:=!)$[/.{)]. 



N 



G"(*°fn = wE 



\Vj — Vn\ 



<^ U^A - ^ {li^) b{9,j)da 









2A^ 



Uj — Vj\ X 



=:/i(^)+/2(V^). 






1+17 
Ml 



da 



As for the first term Ii{V), we have 
1 ^ /• 



«J=i 



2Af2 



V - w| / 6(0) [<p{v*) + (t>{w*) - (l){v) - (l){w)] fi'^ (dv) ^{; (dw) da 



which is precisely the term we are searching for thanks to assumption ( A2) and Lemma r2.13i 
As for the second term l2{V), using that 









-2 



1 

N 



< - I I Y.{ve)'^-']+2{{v,f + {v,f)'^-' 



< 



2^1/2 / ^ 



E(-^)'^ 



-2 



^£=1 



CMi[_,{V) = CMr,,^^{fi^), 



KAC'S PROGRAM IN KINETIC THEORY 69 



we deduce 



\Hy)\ < ^<-2(^)['&]c-(P,^,.^)(^E 




We then use the elementary inequahty 

for a constant C/j^ depending on ki but noi on A^, which yields 

We conclude that (16. 5p holds by piling the two last estimates for /i and l2- 

6.6. Proof of (A4) with time growing bounds. Let us consider some 1-particle initial 
data 

/cfi'o e Pg^. 

Similary as in the previous section, we then define (under the assumption (16. ip on the 
collision kernel) the associated solutions ft and gt to the nonlinear Boltzmann equation 
(|l.ip . as well as 

ht := £f^ [/o] (50 - /o) 

the solution to the linearized Boltzmann equation around /<, given by 

dtft = Qift, ft), f\t=o = /o 

dtgt = Q{gt,gt), g\t=o = 9o 

, dth = 2 Q{ft, ht), h\t=Q = ho := go - fo- 

We also define as before 

wt := gt- ft- ht- 

We shall now again expand the limiting nonlinear semigroup in terms of the initial data, 
around fo- The goal is to prove assumption (A4) with the choice of indices 

{r,',ri") = {l,{l+v)/2). 

This imposes the choice of weight 



A2(/) = Ai(/)^ = ^M,,_2(/). 



Lemma 6.3. For any given energy £ > and any rj > there exists 

• some constant ki > 2 (depending on £ and 5), 

• some constant C (depending on £), 

such that for any 

r G Rf and fo,go e Pg^^r 
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we have for any t > 0; 



(6.6) \\gt - /ilUfi < e^^'+*^ ^max{Mk,-2ifo),Mk,-2igo)} ||/o - 5o|Imi . 

(6.7) \\ht\\,,^ < e^(i+*) ^M,,_2(/o) ll/o -5o|Im1 > 



(6.8) \\oJt\\M^ < e^^'^'^ ^max{Mfc,_2(/o),Mfc,„2(5o)} ||/o - 5o|lij • 

Proof of Lemma \6.3[ . We proceed in several steps and number the constant for clarity. 
Let us define 

yfGM\R% ||/||mi:= / {v)'\f\idv), \\f\\Ml--= [ {v)'{l + log{v)Y\f\{dv). 



k /„j "^"fc/ 



Step 1. The strategy. Existence and uniqueness for ft, gt and ht is a consequence of the 
following important stability argument that we use several times. This estimate is due to 
DiBlasio [22j in a L} framework, and it has been recently extended to a measure framework 
in [29. Lemma 3.2] (see also [32j and [Mj for other argument of uniqueness for measure 
solutions of the spatially homogeneous Boltzmann equation) . 
Let us sketch the argument for h. We first write 



(6.9) ^ l{vf\ht\{dv)< III \ht\{dv)ft{dv,)\u\b{9)\{vY + {v:f-{vy'-{v,y 



da 



+ 2 /// \ht\idv) ftidv,)\u\bie) {v,f da 



(this formal computation can be justified by a regularization proceedure, we refer to [29] 
for instance). Since the first term vanishes, we deduce that 

(6.10) jJhthq<Ci\\fhq\\hthq 

for some constant Ci > only depending on b. 
Then in the case when 

(6.11) ll/slU/i ^ ^ (0)*) o'^ some time interval s G [0,t] 

we may integrate this differential inequality and we deduce that h is unique. 
More precisely, we have established 

(6-12) sup ll/islUfi < ||5o - /oIImi exp ( Ci / ||/s||mi ds ] , 

se[o,i] ^ ^ \ Jo ■^ J 

and similar arguments imply 

(6.13) sup ll/s - 5^11^.1 < \\go - /oIImi exp ( Ci / ||/^ + gs\\jy.n ds ] . 

se[o,t] ^ ^ \ Jo ^ J 

It is worth mentioning that one cannot prove (|6.11|) under the sole assumption 

II/oIImj < °o 

on the initial data since it would contradict the non-uniqueness result of [39]. However, as 
we prove in ()6.15p below, one may show (thanks to the Povzner inequality, as developped 
in [HSliZ]) that (|6lT]) holds as soon as 

II/0IIM2V ^ °°- 
This will be a key step for establishing (j6.6p and (j6.7|) . 
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Now, our goal is to estimate the Mg norm of 

wt := gt- ft- ht 
in terms of \\gQ — /olU/i- The measure ojt satisfies the evokition equation: 

dt(^t = Q{9t,gt) -QiftJt) -QihJt) -Q{ft,ht), ioo = 
which can be rewritten as 

dtujt = Q{(^t, ft + gt) + Q{ht,gt - ft)- 
The same arguments as in (I6.9p - (l6.10p yield the following differential inequality 

^ ll^ilU/i < C2 ||wt|lA/i ll/t + athji + \\Q{ht,gt - ft)hq , llwolUfi = 

for some constant C2 > depending on b. 
We deduce 

/■* 

sup WuJsWm^ < 
se[o,t] 

Since 



/ WQihsJs- gs)\\M^ dsj expf C2 / Wfs + gshq dsj . 



/ \\Q{hs,fs- gs)\\Mi ds<C2 sup ||/is||mi ( / \\9s - fsWui ds] 
Jo ^ \se[o,t] ^J \Jo ^ J 

+ C2 I sup hs-fsWAq] ( / W^sWuids] 

we deduce from (j6.12p and (|6.13p 

(6.14) sup ll^slU/i < C2 \\go - /oIImi exp ( C2 / UfsWul + IbslU/i ) ds ) x 

J \\9s - fshqj exp (Ci J II/.IIM3I 1 

+ U II^^IIa/1 j exp (Ci J (\\fs\\Aq + \\g 



jsWAq 



Hence the problem now reduces to the obtaining of time integral controls over 

II/sIImi ' llfl'slU/i' \\fs-gs\\n,q and ll/isllMgi • 

Step 2. Time integral control of f and g in Mi . In this step we prove 



(6.15) 



/ WfsWM},, dt<Csi£)t + c,\\fohq^ ^ = 1,2, 



for the solution ft, where C3{£) > is a constant depending on the energy, and C4 > is 
a numerical constant. The same estimate obviously holds for the solution gt. 

The estimates (j6.15p are a consequence of the accurate version of the Povzner inequality 
which has been proved in |591 147) . Indeed it was shown in [59, Lemma 2.2] that for any 
function 

^ : M"^ -> M, ^{v) = V'd'yp) with ip convex, 
the solution ft to the hard spheres Boltzmann equation satisfies 

— - / ^{v)ft{dv)= / ftidv) ft{dv*)\v -v*\K,i,{v,v^) 

dt JlRd jRdjRd 
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with K^i, = G^ — H^, where the term G\^ "behaves mildly" (see below) and the term H\^ 
is given by (see [59l formula (2.7)]) 

rTr/2 



H^{v,v^,) = 271 



V^ {\v\^ cos^ 9 + \v^f sin^ 6) - cos^ d i; {\vf) - sin^ d i; {\v^\^) 



dO. 



Note that H^^ > since its integrand is nonnegative because of the convexity of ip. 
More precisely, in the cases that we are interested with, namely 

$(t,) = ^2,c(I^P) with 'il^k,c{r) = r''/'^ [log r)'^ and C = l,2, 
it is established in [59] that (with obvious notation) 

yv,v,e M^ \G.^^Jv,v,)\ < C^{v) (log(t;))^ {v,) (log(t;,))^ 

for some constant C5(C) > depending on (. 

On the other hand, in the case C = 1 we compute, with the help of the the notation 
X := cos^ 6 and u = |f*|/|f |, 



VxG [1/4,3/4], Vug [0,1/2], 

■02,1 (|'y| cos + lf*! sin ( 



cos^ 7/^2,1 (|t;|2) -sin2 002,i(l^'*l^) 



(l- x) V'2,1 {u^) + X V'2,1 (1) - V'2,1 {{I- X)U^ + x) 



> Ca \v\ 



for some numerical constant Cg > 0, which only depends on the strict convexity of the 
real function '02,1- We deduce that there exists a constant C7 > such that 

Similarly, in the case C = 2, we have 

Vxe [1/4,3/4], Vne [0,1/2], 

= 2|t;|2 logl^l^ {(l-x)02,i iu^) +x02,i(l) -V2,i((l-x)'u2 + x)} 



> Cs |u| log |f| 



+ |t;|^ (1 - x) 02,2 (■U^) + X 02,2(1) - ■02,2 ((1 - x) U^ + x) 

for some constant Cg > depending on the strict convexity of ^02,1 and 0^2,2- Hence we 
obtain for some constant Cg > 

H^^,^{V,V^) > Cgl^l^ log|u|^l|„|>2|„.|. 

Putting together the estimates obtained on G2,(j and H2^(^ we deduce 

|7;-'t;^,|i^2,c < Ciq{v)'^ {v^f log (v) log(t'*) - Cn \v - v^\ [til^ (log [til)^"^ l|i'|>2|,;.| 
for some constants Cio,Cii > 0. Since 

\v-v*\ |'t;|^(log|7;|)^"n|^|>2|^,.| >Cst.(v)^ (log(w))'^~^ l|„l>2|^,| - Cst. 

> Cst. {vf (log(t.))^-^ - Cst. {vf {v^f log (v) log (v^) 
we deduce 

(6.16) \v - v,\K2,c < Ci2 {{vf + {v,f) - Ci3 {vf (log(^))'^-i 

for some constants C12, C13 > 0, and we finally obtain the differential inequality 

j^WMImI^ < 2Ci2 {! + £)- Ci3M3,c-i, 
from which (|6.15p follows. 
Step 3. Exponential time integral control of f and g in M^ . 
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This step yields a proof of (|6.6p and ([6 
Let us first prove that 



(6.17) Vt>0, e^^^^^^^^^"("^^"-3^"""-3)'^^<e^"(^+^)* (max{M.(/o),M,(,o)})^ 

for some constant C14 > 0, for any fc > %, with % big enough depending on the energy 
S. 

We shall use the previous step and an interpolation argument. For any given probability 
measure 

/ e Pki^'') with M2U) < £, 

we have for any a > 2 



Ml. 



Ijvf ^i + M^^ (!(.).<. + l(„p>.) f{dv) 

< il + £)(l + i^) +- f {v)^ (1 + log(t;)) fidv) 

< (i + f)(i + ^) + i|l/llM. 

where we have used inequality logx <j; — lonx>lin the last step. 
By choosing 

«^=ll/llLi> 

we get 

(6.18) ||/|U,i^<2(l + £:) (l + log||/|U,i 

Remark 6.4. Observe here that it was absolutely crucial to be able to control the right 
hand side of ()6.15p in terms of the Mg ^ moment, that is only a logarithmic loss of moment 
as compared to M2 ■ This is what allows us to control this right hand side in terms of the 
logarithm of a higher moment of /, so that the exponential in ()6.17p can be controlled 
in terms of some polynomial moment of /, hence fulfilling the requirement on the loss 
of weight of the stability on the semigroup for the abstract method. Recall indeed that 
the moment associated with the weight Ai has to be propagated along the flow of the 
A'^-particle system. And we have been unable to control an exponential moment for such 
a high-dimension evolution. 

On the other hand, the following elementary Holder inequality holds 

(6.19) y^y E N, fc' < K v/ G Ml ii/iu,., < ||/||i;f /'= \\f\\%' < ||/||J;/^ 

Then estimate (16.170 follows from ()6.15p . ()6.18p and ()6.19p with k' = 5 and k > 5 large 
enough in such a way that 

(Ci + C2)C42(l+£:)|<i. 

We then deduce (|6.6p from (|6.13p . and (similarly) (|6.7p from (|6.12p . 

Step 4. Time integral control on d and h. Let us denote as before 

dt ■■= ft - 9t- 
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Let US prove 

ft 



(6.20) (^^ ||4|U,i ds^ and ^ 



I'^sIIm-'- 



Cl/o(ilWlMl+ll9=llMl) 



(is 



<Ci5 IMoIIm] e ^"V -3 -3; (^C3(f)t + C74 \\M\miJ + \\Mmi,- 

for some energy dependent constant C^ and some numerical constant C". Performing 
similar computations to those leading to ()6.9p . we obtain 

^ylll^tllMii < \ht\{dv)ft{dv^) \v -v^\ K2,i{v,v^) 

+C16 / / / |/it|(dt;)/i(dT;*) |w-i;*| (u*)^ (l + log(f*)) 

for some constant Ciq > depending on b. Thanks to the Povzner inequality (j6.16p (with 
C = 1), we deduce for some constants C17, Cig > 

Integrating this differential inequality yields 

ll^tllMi, +^"18 J WhsWul ds < Cn f sup WhsW^i J U II/sIImi, dsj + ||/io|Im1i • 
Using the previous pointwise control on ||/it||Afi and (j6.15p (with C = 2) we get 

/" III. II J ^ ^16 II , II C*! /oII/sIImI ^* //-» /C\i I /-< II -P II \ I II J II 

J \\hshq ds<-^ Wdohq e "3 l^CsiS) t + C4 ||/o|Imi_J + MoIImi, ■ 

Arguing similarly for dt, we deduce ()6.20p . 
Step 5. Conclusion. We first rewrite ()6.14p as 



sup WOJsllMi 

se[o,i] 



, II (Ci+C2)/„*(||/.!| i+iis.ll 1) 
MlMi e V 3 37 



ds 



< (^2 IMollA,fi e V 3 ^y ( / ( ||4||;,,^i + ||/ij|),fi ) ds 



^2 ||«0|lAfi e V 1/ l^ll"*llA/l "^ II"''^IIa/i 



Then we use the estimate (|6.2Up for the last term: 



sup \\UJs\\Mi 

se[o,t] 



(2Ci+C2)/o'f||/.l|„i+l|9.||„i)ds 

,1 e V 3 3/ 



< C2C15 ||do|U,i NolUfi, e ^" V "3 -37 (^C3(£:)t + C4 ||/o||m1 

Finally we use estimate (j6.17p for the exponential term with k = ki and we obtain 

sup ll^slUfi <C'2Ci5 ||(io|lA/i IMolUfi X 
se[o,t] ^ ^'^ 

X e^"(i+^)* (max{Mfc(/o),Mfe(5o)})^ (C3(£:)t + C4 ||/o|Imi,J • 

Then arguing as in the end of Step 3, for any r] £ (0, 1), using ()6.19p with fci large enough, 
we have 

IMoIImj, <(max{Mfc(/o),Mfc(go)})^ 11^011^7 
and 

II/oIImi, < (max{Mfc(/o),Mfc(5o)})^ • 



KAC'S PROGRAM IN KINETIC THEORY 75 



We conclude with 

sup \\ujs\\Mi ^ ^^^'''*^ x/max{Mfc(/o),M,,(5o)} ||/o - go& , 
se[o,i] 

from which estimate (16.81) foUows. D 



6.7. Proof of ( A4) uniformly in time. Let us start from an auxihary result from [60] 
Let us define the linearized Boltzmann collision operator at 7 

£^(/) = 2Q(7,/) 

where 



7 



is a maxwellian distribution with momentum u G M*^ and temperature T > 0. 

Theorem 6.5 (Theorem 1.2 in [60]). First the linearized Boltzmann semigroup e^* for 
hard spheres satisfies 

(6.21) l|e^^*llri/ -M <C,e-^^ 

wj/iere 

m^(w) :=e^l''l, z > 0, 

and A = A(T) is the optimal rate, given by the first non-zero eigenvalue of the linearized 
operator C^ in the smaller space L^(7^^). 

Second the nonlinear Boltzmann semigroup S*^^ satisfies 

(6.22) \\srifo) - 7|Li(„,i) < C. e-^* ||/o - 7lLi(^-) 
where 7 = 7^^ is the maxwellian equilibrium associated with /q: 



-f{v) 



{2'kTY 
with 

yi = l,...,d, Ui = foVidv and T = - fo\v\ dv 
jRd d JiRd 

and Cz is some constant possibly depending on z and the energy E of the solution consid- 
ered, and A = A(T) is the same rate function as before. 

Let us now prove uniform in time estimate for the expansion of the limiting semigroup 
in terms of the initial data. 

Lemma 6.6. For any given energy £ > and for any r] G (0, 1), there exists 

• some constant %^^ > 2 (depending on E and r\), 

• some constant C (depending on 8), 

such that for any 

fo,goePg,{R'') 
satisfying 

\/i = l,...,d, {fo,Vi) = {go,Vi) and (/o, |f P) = (50, bP) < <?, 
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and for any k > k^^^i we have 

(6.23) \\gt - /ilUn < Ce-t* ^max{Mfc(/o), Mk{go)} \\go - /olli^ , 



(6.24) \\ht\\Mi<Ce-2'^Mk{h) ||<7o - /o 



|1— r; 
'M. 



1 ) 



(6.25) Wuthn < Ce-t Vmax{Mfc(/o), Mk{go)} \\go - /oHL". 
Note that these estimates imply (A4) with T = oo, Pg^ = Pg^. 

Remark 6.7. Observe that in the following proof we shall use moment production hounds 
on the limiting equation. Indeed once stability estimates for small times have been secured 
(as in Lemma l6.3p . one can use, for t > Tq > 0, moments production whose bounds only 
depend on the energy of the solution. This, together with the linearized theory in L^ 
setting with exponential moment bounds of Theorem 16.51 will be the key to the following 
proof. 

Proof of Lemma \6.6l From [481 Theorem 1.2-(b)] (see also [Ij for a simpler different proof), 
there exists some constants z, Z (only depending on the collision kernel and the energy of 
the solutions) such that 

(6.26) sup\\ft + gt + ht\\Li < Z, m2z{v) := e^'\^\ 

t>l 2z 

(note that the proof in [48] applies to the solutions ft and gt, however it is straightforward 
to apply exactly the same proof to the linearized solution ht around ft, once exponential 
moment is known on ft). 

We also know from ()6.22p that (maybe by choosing a larger Z) 

(6.27) Vt>l ||/,-^||^^^^+||5,-^||^^^^ <2Ze-^*, 
where 



7 := 7/0 = 790 



{2-jrTy 
with 

\/i=l,...,d, Ui = {f,Vi) = {f,Vi) and T = -{f,\vf) = -{g,\v\^) 

stands for the normalized maxwellian associated to /o and gQ. 
We write 

dtift - gt) = Qift - gt, ft + gt) 

= ^jift - gt) + Qift - gt, ft-i) + Qift - gt,gt - 7) 

and, using also (|6.21|) on the linearized semigroup, we deduce for 

-"(*) := Wft-gtWui,^ 

the following differential inequality for t > Tq > 1 

uit)<Ce~^^'~^^Kin) 

+ C Te-^^*-^) (\\Qifs-gs,fs-i)\Wn^m-.) + \\Qifs-gs,gs-i)\\M^^.)) ds 

J To 

(this formal inequality and next ones can easily be justified rigorously by a regulariz- 
ing proceedure and using a uniqueness result for measure solutions such as [32l [291 SH])- 
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Therefore we obtain 
u{t)<Ce-^^'-^'''^u{To) 



+ C' / e ■" [\\fs-7\\M^{{v)m^) + \\9s-7\\M'^{{v)m^)) Wfs- 9s\\Mmv)m,) ds. 

J To 

We then use the control of M^{{v) m^) by M^{m2z) together with the controls (|6.26p - 
(|6.27p . the decay control (|6.2ip and the estimate 

g-As-A(t-s) ^ g-|t-|s^ 

We get 

u{t)<Ce-'^'-^^^K{To) + Ce-^' f e'^ ' \\fs - gs\\M^iv)m.) ds. 

J To 
We then use the following control for any a > 0: 



{v) mz J 

<af \f-9\e^\''\+e'^- f (f + g) 



e2z|i,| 
<au{t) + e-"'Z. 



Hence we get 



\\f-g\\Mi, < 

{v) mz 

and we deduce 



n + e ^ Z < [1 + Z)u when u>\, (choosing a := 1) 
- I log u\u + uZ when u <1 (choosing — za := log u) 

\\f-g\\M}, <Ku{l + {\ogu)^), K:=l + - + Z. 

(v) mz Z 

Then for any 5 G (0, 1), we have, by choosing Tq large enough, 

Vt>ro, e~t*<5e-^* 
and we conclude with the following integral inequality 

(6.28) u(t) <Ce-^(*-^")u(To)+5e-^* / e'^' Us {l + {log Us)_) ds. 

Jto 

Let us prove that this integral inequality implies 

(6.29) Vt>ro, u{t)<Ce-i*u{Tof-^. 
Consider the case of equality in (|6.28p . Then we have 

n(t)>e-^(*-^»)u(To) 

and therefore 

(l + (logni)_) < (l + (logu(To))_ + A(t-ro)). 
We then have 

U{t) := / e-2'us (l + (log'u,)_) ds 
Jto 

< / e~^ 'Us {I + {log u{To))_ + X{s- To)) ds 
JTo 

< (l + (logn(ro))_) / e-^n, 

JTn 
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By a Gronwall-like argument we can therefore obtain 

n(t)<Ce-^(*-' 
Then thanks to the inequahty 



u{t)<Ce-^^'~^°^u{To) + C5e-^' (l + (logu(To))_) u{To) 



Vxe(0,l], -(logx)x< 



x^-^ 



6 

we can prove (|6.29p when u(ro) < 1, and in the case when u(Tq) > 1, we can use (|6.26p 
again to get 

u{To)<{2Z)^u{To)^-^. 

This concludes the proof of the claimed inequahty (|6.29p . 

Then estimate (j6.23p follows by choosing 6 small enough (in relation to r/) and then 
connecting the last estimate (j6.29p from time Tq on together with the previous finite time 
estimate (|6.6p from time until time Tq. 

Then the estimate (j6.24p is proved exactly in the same way by using the equation 

dth = C^ift - gt) + Q{ht, ft - 7) 

(which is even simpler than the equation for ft — gt). 

Concerning the estimate (I6.25P we start from the equation 

dtUJt = '^C^tiht) + Q{(^t,ft -7) + Q{^t,gt -7) + Q{ht,dt). 
Then we establish on 

y{t) := W^iWmI^^ 
the following differential inequality 

y{t) < C e-^ (*-^«) yin) + C 5 e-^ (l + (log y(ro))_) y(To) 

which implies 

y{t)<Ce~-^' (y{To)'~' + WdrJll' IIHIIm') " 
Then estimate ()6.25p follows by choosing 5 small enough (in relation to rj) and then 
connecting the last estimate from time Tq on together with the previous finite time estimate 
(16. 8p from time until time Tq. D 



6.8. Proof of (A5) uniformly in time. Let us prove that for any ZjTWj £ (0, c«) there 
exists some continuous function 

e : M+ ^ M+, 9(0) = 0, 

such that for any 

/o,5oG^m.(M'^), m,-(t;):=e^"H 

such that 

||/o|lAf4_ <-^5, ||50|Im1,_ <-^5. 

there holds 

supH^i (5f^(/o), 5f^(<7o)) < e {Wi (/o, 50)) . 
t>o 

where Wt stands for the Kantorovich-Rubinstein distance. Let us denote 

Wt:=Wi[S^\fo),S^\go)). 
As we shall see, we may choose 

(6.30) e(u;):=min|e, Gei-^^+l'^s'"!))'^', — — + C2wA, 6(0) = 0, 

[ (1 + |logu;|)2if J 
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for some constants 0, Ci, C2 > (only depending on z and M.z)- 
We start off with the inequahty 

vt > Wt< ||(/t - 5t)bllUp < \ Il(/* + 5t)(«)lMi = 1 + ^ =: e. 

Let us now improve this inequahty for smah value of Wq. Therefore we assume without 
restriction that 

in the sequel. 

On the one hand, it has been proved in ^321 Theorem 2.2 and Corollary 2.3] that 

(6.31) Wt<W^^K I W,(\^ (log W,)^) ds, 

Jo 

for some constant K. 

One can then check that the function 

W{t) := e^--"' (lyo)^"" 

satisfies 

W'{t) = K {I- log W{t)) W{t), W{0) = Wq. 

Therefore it is a super-solution of the differential inequality ()6.3ip as long as Wt < 1. It 
is an easy computation that this super-solution satisfies 

log(l + |logP7o|) 



K 



W{t) < 1 as long as t < to := 

Observe also that W{t) is increasing on t e [0,to]- 
We then define 

...,^^ , to log (1 + 1 log ^ol) 

(6.32) t,:=- = 

and we deduce the following bound on the solution of (|6.3ip : 

(6.33) VtG[0,ti], Wt<Wt<Wt,=e^-^^+^^°^^''^^^'^\ 

On the other hand, from ()6.22p . there are constants A, Z > 0, z G (0, z) such that 

(6.34) Vt > ll/t -7/oIIl^, + \\9t-7,o\\Ll^^ < ^ e'^\ 

where 7/g and 7gg stand again for the normalized maxwellian associated to /o and g^. 
We denote by wj^, and Ug^ the momentum of /o and g^ 

Uh = {fo,v), Ug,^ = {go,v), 

by Tjq and Tg^ the temperature of /o and go 

Tfo = ^ (/o, \v - Ufgfj , Tg, = - (^go, \v - Ug^fj , 
and by Sf^ and Sg^ the energy of to /o and go 

£fo = {fo,\v\^) , £go = {9o,\v\^)- 
Then we can use (14. 5 p for some 

s > large enough and a ^ ( 0, — 

V 2s 
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and we compute 



(6.35) 



< 

< 

< 

< 
< 
< 




C {W2ifo,gof'' + Wiifo,gof'') 
C(W^i(/o,5or + m(/o,5o)'") 



where the final constant C > depends on s > as well as on upper and lower bounds 
on the temperatures of /o and qq. 

Remark 6.8. We also refer to [21J for more general estimates of the Wasserstein distance 
between two gaussians. Such refined estimates are however not needed in our study. 



Then we combine (I6.34p and (16.350 by triangular inequality to get 

(6.36) V t > Wt < Ci e"^ * + Ca VFo" 

for two constants Ci, C2 > 0. 

We then consider times t >ti and we deduce from (16.320 and (16.360 the following bound 
from above 



(6.37) 



Vt > ti, Wt< Ci e-^*i + C2 Wn" 



+ C2 w;^. 



(l + |logWo|)5^ 

It is then straightforward to conclude the proof of ( A5) uniformly in time for the function 
K^ by combining (lUIHHD and (lOTj) . 

By using Theorem 13.11 whose assumptions have been proved above, this concludes 
the proof of point (i) in Theorem 16.11 together with the the estimate on W^ (/) from 
Lemma 14.21 

One can also conclude the proof of point (ii) in Theorem 16.11 by using 

• Lemma 14.41 for the construction of the sequence initial data /g which satisfies the 
required integral and support moment bounds, 

• The previous in order to apply Theorem l3.lt 

• Lemma 14.71 in order to estimate 

6.9. Proof of infinite-dimensional Wasserstein chaos. Let us now prove Theorem l6.2[ 
Its proof is similar to Theorem 15. 3[ 

First the proof of (j6.3p as a consequence of point (i) in Theorem 16.11 and [371 Theo- 
rem 1.1] is strictly similar to the one of (j5.8p as a consequence of point (i) in Theorem 1 5. II 
and [571 Theorem 1.1]. 

Then the proof of (j6.4p is also similar to the one of (j5.9p , the only difference being that 
one needs the following result of lower bound (independent of N) on the spectral gap of 
the A^-particle system. 
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Theorem 6.9 ([llj). Consider the operator Lhs for the hard spheres N-particle model 
with collision kernel B{v — w) = \v — w\. Then there is a constant A > such that for any 
probability f on S one has 

where S is the hypersurface 

S^:=\V€E^;^J2^, = and ^EN|' = i|. 
I fc=l k=l ) 

and for some constant A > independent of N . 
Then using Theorem 16.91 we deduce that 

VA^ > 1, Vt > ||/l^ - l|L2(5iv,^iV) < e-^* \\K - l|lL2(5iv,^iV) , 

where h^ = df^/d'y^ is the Radon-Nikodym derivative of f^ with respect to the measure 
7^ and the end of proof of (j6.4p is then exactly similar to the one of (|5.9p . 

7. H-TYIEOREM AND ENTROPIC CHAOS 

This section is concerned with the i:f-theorem. We answer a question raised by Kac ^42j 
about the derivation of the i7-theorem. 

7.1. Statement of the results. Our main results of this section state as follows: 

Theorem 7.10. Consider the Boltzmann collision process for Maxwell molecules (with or 
without cutoff) or hard spheres, and some initial data 



/o G L°° (W^) s. t. I e^l^l dfo{v) < +00 



for some z > 0, and the sequence of N-particle initial data (/o^)Af>i constructed in 
Lemma Um and \4-T\ 
Then we have: 

(i) In the case of Maxwell molecules with cut-off and hard spheres, if the initial data 
is entropy- chaotic in the sense 



N 
with 



^ H{fi^\^^)^^^±^Hifoh). 



J5^ 



df^ 



N 



d-i^' 

then the solution is also entropy chaotic for any later time: 

Vt>0, ^H{f^\^^)^^^^±^H{f,\^). 

This proves the derivation of the H -theorem in this context. 
(ii) In the case of Maxwell molecules, and assuming moreover that the Fisher informa- 
tion of the initial data /o is finite: 

|V./o|' ^^ 

■ dv < +00, 



/o 

the following estimate on the relaxation induced by the H -theorem uniformly in the 
number of particles also holds: 

viv>i, l.H{fi'\^'')<m 

for some function f3{t) — )■ as t — )■ oo. 
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Remarks 7.11. (1) The assumptions on the initial data could be relaxed to just P4n-L°° 
as in point (iii) of Theorem 15. li However this stronger allows to apply the previous 
theorems for hard spheres and Maxwell molecules. We do not search for optimal 
statement here, but rather emphasize the ideas of our method. 

(2) A stronger notion of entropy chaoticity could be 

1 H {n[f^%^) ^^^^ 0. 

The propagation of such a stronger property is an interesting open question. 

(3) The point (ii) holds for the hard spheres conditionally to a bound on the Fisher 
information uniformly in time and in the number of particle. However at now, 
proving such a bound for the many-particle hard spheres jump process is an open 
problem. 

(4) In point (ii) one could expect to have 

ViVGN*, Vt>0 ^H{ft^\j^) <Ce-^^. 

7.2. Propagation of entropic chaos and derivation of the //-theorem. In this 
subsection we shall prove the point (i) of Theorem 17.101 Its proof relies on a convexity 
argument. 

Let us define h^ := df^/d'y^ and then compute 



I f ^ f .,v.... ,v...^, h^m 



■^ t E l_^ i^m) - ^nV)) log j^ Biv. - .,, .) d. ,-idV) 



where we recall that V*, was defined in (jS.ip . which implies 



ft 



(7.38) Vt>0, ^H{f,^\^^)+I^D^{f^)ds = ^H{f,^\j^). 

We also note that the same kind of equality is true at the limit (see e.g. [17] ) 

Vi>0, H{fth)+[D'^{f,)ds = H{Mj) 
Jo 

with 

D^ if) -.= 1 I iff: - ff*) log ^ B{v -v,,a) dv dv, da 

(be careful to the factor 1/2 in our definition of the collision operator p.2|) when computing 
the entropy production functional). 

We then have the following lower semi-continuity property on these functionals, as a 
consequence of their convexity property. 

Lemma 7.12. The many-particle relative entropy and entropy production functionals de- 
fined above are lower semi- continuous: if the sequence {f^)N>i is f -chaotic then 

liminfli/(/^|7^)>i/(/|7) 



and 



liminflD^(/^) >!)-(/). 
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Let us first explain how to conclude the proof of point (i) of Theorem 17.101 with this 
lemma at hand. We first deduce from (|7.38p and the entropic chaoticity of the initial data 
that 



V 



t > 0, j^H (//^|7^) + 1^ D^ (/f ) ds ^^^^ H (/o|7) 

= H{ft\-f)+ [ D^{f,)ds 
Jo 

Second we use Lemma 17.121 on the LHS to deduce that 

Vt > 0, linUnf (1 H {ff\^^) + J^ D^ (/,^)) > H {ft\^) + J^ D^ (/,) ds 

where each of the limit of the two non-negative terms on the LHS is greater that the 
corresponding non-negative term in the RHS. We deduce from the two last equations that 
necessarily 

and 

Vt > 0, fo^ (/f ) ds ^^^ f D^ (/,) ds 
Jo Jo 

which concludes the proof of point (i) of Theorem 17.101 



Proof of Lemma 7.12\ These inequalities are consequences of convexity properties. The 



lower continuity property on the relative entropy on the spheres was proved in |12t Theo- 
rem 12] (actually the proof in this reference is performed on the sphere S^~^, but extending 
it to the invariant subspaces of our jump processes S is straightforward). We refer to 
|16j for a detailed proof of the latter. 

Let us now prove the inequality for the entropy production functional D^ . Denoting 
Z = h^ {y*2)/h^ , we first rewrite thanks to the symmetry of f^ as 

^^(n=^^^|^/ / J{Z)B{v,-v,,a)fnv„v,)daJ^im_ 
2N^ JsN Jsd-i fi {n,V2) 

where J{z) := {z — 1) logz and /2^ denotes the 2-marginal. Since the function z i— )■ J{z) 
is convex, we can apply a Jensen inequality according to the variables V3,...,V]\[ with 
reference probability measure f^/f2^, which yields 

D"" (/^) > ^^f^~ ^^ / / J{Z)B{v,-V2,cj)fi'{v,,V2)dadv,dv2 

with 

f'iV) _f^{{V^2)l.(y^2)2) 



Z{vi,V2) ■■= Z ■ 



lv3,...,VN&S'^{v^,V2) fi^ {V1,V2) fi^ {vi,V2) 

where S^{vi,V2) := {vs, . . . ,V]\[ G E^~'^, {vi, . . . ,vn) G S^}. We therefore deduce a 
control from below of the A^-particle entropy production functional in terms of the 2- 
particle entropy production functional, denoting (/^ = /^ ((^1*2)1) (^1*2)2))- 

D''{n>^%^l I mr-fi')log^-^B{v,-V2,a)dadv,dv2 

Finally we take advantage of the convexity of the functional 

X 

h{x,y) = {x-y) log - 

y 
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which iinphes that the function 

(/2,fi'2)->/ / {f2- g2)log — B{vi-V2,a) 

is lower semi-continuous for the weak convergence of the 2-particle distributions /2 and 
(72 as proved in [241 Step 2 of the proof] . Hence we obtain thanks to the chaoticity of the 
second marginal 

liminfZ?^(/^) 



which concludes the proof. D 

7.3. Many-particle relaxation rate in the i^-theorem. In this subsection we shall 
prove point (ii) in Theorem 17.101 Its proof goes in two steps. First we shall prove that 
it follows from an estimate on the Fisher information thanks to the so-called "HWI" 
interpolation inequality [72]. Second we shall prove such a uniform bound on the Fisher 
information in the case of Maxwell molecules. Let us take the opportunity to thank 
Maxime Hauray who kindly communicated to us a proof for the latter step. 
Let us define the Fisher informations for the A^-particle distribution: 

jRdN J 



N 



and 

for a probability / having a density with respect to the Lebesgue measure in M and 
with respect to the measure 7 respectively. The gradient in this formula has to be 
understood as the usual Riemannian geometry gradient in the manifold S . The tangent 
space TSy (of dimension Nd — 2) at some given point V G S is given by 

TS{f = iw gR"^^ s. t. "^wi = and VF ± F I . 

For more informations and other results on the Fisher informations on S we refer 

to [5] . We shall prove the following lemma whose proof is inspired from [7T] ) . 

Lemma 7.13. Consider the N -particle jump process for Maxwell molecules as defined in 
Subsection 15. 1\ for a give N > 1, and some initial data /q^ whose Fisher information is 
finite /(/q I7 ) < +00 on S . Then one has the following uniform in time bound on the 
Fisher information of the solution 

Proof of Lemma \7.13\ We shall first consider the case of cutoff Maxwel molecules whose 



collision kernel b is integrable, and for a positive and smooth solution /^ on 5^. These 
assumptions can be relaxed by a mollification argument. 

It is possible to study directly the estimate to be proved on the manifold S'^ , however 
it means that one has to consider some local coordinates and local basis for the tangent 
space. Another simpler method is to take advantage of the fact that the dynamics leaves 
the energy unchanged. 
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Starting from an initial data /g' on 5^ we consider the flatened initial data 



N 



with £{V) 



1=1 



and m{V) 



^/W) J 



N ' ' N 

Observe that from the conservation of energy and momentum and the uniqueness of the 
solutions to the linear master A^-particle equation 



Vi>0, fl' ■.= ai£{V))l3{m{V))fl 



^ V-m{V) 



VW) 



where //^ denotes the solution in W^^ starting from f^. If the functions a and /3 are 

regular and compactly supported, as well as /g^, this produces a smooth solution on R'^^. 

Assume that the result on the Fisher information is true in 



I f^ 



cN 



v/f 



ir 



■dV < 



V/o^ 



fN 
JO 



dv = nfn- 



cN 



Then we have the orthogonal decomposition of the gradient locally in terms of radial and 
ortho-radial directions 

that we can plug into the Fisher information inequality: 



I ( ft ) := (Vf a + Vm^y + ( / a{£) f3{m) 

ie,m J Js^ 



< (V^Q + V^/3)^ + ( / a{£)^(m) 






r{dV) 



SN 



C 



i^\dv) = i[fn. 



pN 



Dropping the terms which do not depend on time we obtain the desired inequality on 5^. 
Let us now prove the inequality on W^^ . Let us first fix some notation: the A'^-particle 
solution //^ satisfies 



1 ^ 

- y 



^^f'" =Tt 1^ if i^^jAV)) 6(cos%) da - f^{V)) da 



-: NB (Q+'^(r) - r) 



where we use the following notations. We define 



Q-^'^n-^ E Qtfif)^ Qp''in--=[ f!^b (cose.,) da, 

COS 9ij := a ■ kij with kij = {vi — Vj)/\vi — Vj\ 
and where we assume that b is even and that 



/ b{a ■ k) da = Cb for any k, 



1. 



For any function g on M shall use the shorthand notation gf^ to denote the function 
V I—)- g{rij^(j[V)), which depends also implicitly on a. We shall make use of the measure 
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preserving involution 

' r X M" X S'^-^ -^ M" X M" X S^-i 

where a' = {vi — Vj)/\vi — Vj\ = kij. 

Finally as in |71) , we shall use the following endomorphism of M"^ 

Makix) = {k ■ a)x — {k ■ x)a 

Pak{x) = {a-x)k + M„k 

and we recall that ||P(jfc(3;)|| < ||x|| with equality only if x,(T, A; are coplanar. 

We claim that it is enough to prove that 

(7.39) /(Q+.A.(/A^))=/[ 1 ^ I f{n,,„{V))h{cosd,,)da\ <CBl{f) 

Indeed with this result at hand, we can write for e > 0: 

flis = e-^^-^ fl" + NCb r e^^-(-^) Q+'^ {fH,) ds 

Jo 

and therefore from the convexity of / 

/ ifiis) < e-"^- / an + (1 - e--^-) / (^[ Q-'- ifii.) ^fffX^ 

Observe that 

r NCBe^^Bis-e) 



ds . 



rfs = 1 



and then we can use the convexity of / again to get 

I {fits) < e-'"''-' I Un + r HQ^'"" iftls)) NCbc^^-^^-^Us. 

Jo 

Finally using the claimed result (j7.39p we obtain 

I{ftle)-I{fn < _ (^-^"'^''"') /(//V) + l r I [fits) iVC|e^^^(^-^)d5. 
e £ e Jo 

Then taking e — )• and using Lebesgue's theorem we deduce 



d 



j^i{fn<-NCBi{fn+NCBi{fn<o 

which concludes the proof. 

Let us now focus on the proof of the claim (j7.39p . Taking advantage of the convexity 
of /, it is enough to prove 

v./jG[ii,ivi], /(Qj'^(n)<c^/(n. 

Let us compute each partial derivative of Q.' (/^)- If ^ ^ {hj} then the derivative 
does not act on the kernel b and we obtain: 

v., (qJ'^I/^^)) = / V.,{f!J) b {cos e,j) da = f {V,J^) b{cose,,)da 
= 2/ (/F) (V,,v7^) b {cos 9,,) da. 
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If i G {hj}, then it is slightly more complicated. Without restriction we perform 
calculations in the case i = i. Let us first prove the formula 



(7.40) V^AQIfif 



)+,N ( fN 



"^vj )ii + (v.,r ),, + p.k (v.,r) .. - (v.^r )„ mcos^,,) da 



(the same equality obviously holds where i is replaced by j) 
Simple computations (see for instance [TTj) yield 



eN\ ^ f /-ri fN\ , 1^ fN\ \ , 1 / /vT fN\ 1^ fN\ 



v.(/^) = ^ (v.rL+(v.,nj+^ (v.n.,-(v„,nj-a 






V„. [^(cos^)] 



v.,n,,-(v.,,n 



b' {a -kij) Ilf.±a, 



where n^± is the projection on the hyperplane k-^. Using the first and third equality 
above, we get 



(7.41) v^AQtrif 



1+,Af / j:N 



6(cos%) 1{V,J) +{V,J) + nv,J\^-{V,j)A-a 



'ij 



b' {cos 9 ij) 



Vi - Vi I 



A; I da 
Ilj.±a da 



and we use the following formula of integration by part on the sphere S'^ ^ (see |71t Lemma 
2]) 

/ b' (cos 0ij) F{a) U^,±ada= b (cos %) M^k (V^F(cr)) da 

and the second equality above to rewrite the term involving b' in (j7.4ip into 

i ^^_^ 6(cos%) M^k ((V.J),^. - {V.J)^^) da 

Putting all together, we get formula (|7.4U|) . 

We deduce that for £ / i,j we have by Cauchy-Schwarz 



v.. (qS'^ (/^ 



< 4 



/ f p (cos Oji) da) if 



V \/7^ 



«j 



6 (cos^jj) da 



and therefore 



v.. (qJ^ (/ 



Q.^-'^ (/^) 



< 4 



v.,/F 



u 



6 (cos 9ij ) do". 



88 



S. MISCHLER AND C. MOUHOT 



Now integrating in V we obtain 

V -^ / J^dN Jgd-1 

< 4 



^ VfV J 






b [cos 6 ij) da dV 
b {cos 6 ij) da dV 



< ^Cb 



V 



Vl 



m dV =: I, (/ 



N\ 



where we have used Cauchy-Schwartz's inequahty and the change of variable @ij, and 
where Ii is defined from the last line (Fisher information restricted to the i-th derivative) . 
When i = i,j, we use (j7.40p to get 



v.. Qj^ (/ 



< 



f'^b {cos 9ji) da 



v.y/^ + V,V/ 



f-N 



+ Pak V.V/ 



N 



v.,V/ 



N 



b {cos 9 ij) da 



where we have used Cauchy-Schwartz's inequality and the change of variable Qij 
Since for fixed V, Pak is odd in a and b{cos6ij) is even in a, we have 



/ A-P,k{B)da = 



for any functions A, B independent of a. Using finally that Pq-^ has norm less than 1 (for 
the subordinated norm to the euclidean norm) we get 



v.. (Qtf if 



< 2 



f'^b {cos Oji) da 



'^U.Vf 



N 



+ 



v.,V/ 



N 



b {cos 9ij) da 



and therefore 



Finally we end up with 

N N 



£=1 



which concludes the proof. 



D 



Let us now conclude the proof of point (ii) in Theorem 17.101 We make use of the so- 
called "HWI" interpolation inequality on the manifold S . Observe that S has positive 
Ricci curvature since it has positive curvature. Then [.73', Theorem 30.21] implies that 



^(;^|^^)<lf2(/^7^)./nFb^ 



1 

N" y-' " J - ^ V N 

We can then use the uniform bound on the Fisher information provided by Lemma f7.13l (vi) 
and the bound on the initial data to get: 



N 



< 



N 



< C 
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for some constant C > independent of N. Moreover Lemma |4. II and the propagation of 
moments on the A^-particle system in Lemma 15.41 imply that 



1 -^ N 



for some constant C > and exponent q > independent of N. Then using Theorem 
(case (b)) we deduce that 



with j3{t) — 7- as t — )■ +oo, which implies that 

and concludes the proof. 

8. The BBGKY approach revisited 

The so-called BBGKY (Bogoliubov, Born, Green, Kirkwood and Yvon) and statistical 
solutions approach is very popular in physics and mathematics for studying many-particle 
systems: see for instance [3j where this approach is used for Kac's master equation for hard 
spheres, or see, among many other works, the recent impressive series of papers [U [26l [271 
[28] where this approach is used for the derivation of nonlinear Schrodinger equations in 
mean-field theory in quantum physics. The basic ideas underlying the BBGKY approach 
to mean-field could be summarized as: 

• Write a BBGKY hierarchy on marginals of the A^-particle system and prove that 
the A^-particle system solutions converge to the solutions of an "infinite hierarchy" 
when A^ goes to infinity. 

• Then prove that solutions to this infinite hierarchy are unique, which is the hardest 
part of this program. 

• Then deduce the propagation of chaos by exhibiting, for any chaotic initial data 
to the infinite hierarchy, a natural solution obtained by the infinite tensorization 
of the 1-particle solution to the limiting nonlinear mean-field equation. 

This section revisits the so-called BBGKY hierarchy and statistical solutions approach 
in the case of collision processes, by showing (1) how these notions are included in our 
functional framework (cf. the abstract semigroup T^ defined in Section [2]), and (2) how to 
give a proof of uniqueness and propagation of chaos based on them by using the functional 
tools we have introduced. 

8.1. The BBGKY hierarchy. Let us recall the master of the A^-particle system under- 
going a Boltzmann collision process (the notion of BBGKY hierarchy has wider application 
range, but we shall stick to this concrete case for clarity), see (|5.3p - (j5.4p : 

(8.1) 5,</,^,^> = (//^,G^c^> 

with 

{G^ip) iV) = ^ Y. r ih -v,\) [ 5(cos %) [(^:, - ip] da 

l<i<J<N J^"-' 

where 99*^- = ^ {V*^) and 99 = ^{V) e C^ {m.^'^^ . 
Then the BBGKY hierarchy writes as folows. Let us recall the notation 

Jvt+i,...,v^ 
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for the marginals. Then integrating the master equation ()8.ip against some test function 
ip = (^(ui, . . . ,V() depending only on the first £ variables leads to 



|(/r,^> = (/."^i,Gf+iM> 



where 



Gf+i(<^) := ^ Y. r {\v, - v,\) f b{cose,,) [99*^. - y,] da 

I\ /sd— 1 



Then denoting by 



Z^:=(f'',r{\v,-v,\) J^^_^b{cose,,) [^*j-ip] da^ 



we can futher decompose this sum as 



d_ 
di 



i,j<e i<i<j i,j<i 



by observing that Z^j = for i,j > £, and using the symmetry of /^ we finally deduce 

|(/,^.)^^(ea^.c(^). 

We thus end with a series of N coupled equations on the marginals /^ , where the 
^-equation (^ < A^ — 1) depends on the f^-^ marginal. 

8.2. The infinite hierarchy and statistical solutions. Assume now that 

y(> 1, /f ^vTf in P^™°'^ 



Starting from (18. 2p we obtain, for 93 = ip{vi, . . . ,V£) depending only on the first i 
variables, 

— (7r^,99) = (7r^+i,G^i(v9)) 
where Gi+i{ip) e Cb(M^(^+i)) is defined by 

Gr+,{^):=Ti\v,-v,\) f 6(cos%) [ip*,-^] da. 

In a more compact form, we have the following set of linear coupled evolution equations 

(8.3) y£ > 1, dtTTe = A,^i (7r,+i) with ^^^ := {GT+i)* • 

Since the family of ^-particle probabilities vr^ is symmetric and compatible in the sense 
that 

y£ > 1, n^ [vr^+i] =7Te 

(this follows from the construction), we can associate by Hewitt-Savage's Theorem [38] 
a unique tt G P(P(IR'^)) and this set of evolution equation translates into an evolution 
equation 

dtTr = A°°{7r) on P (p (rA] 
of statistical solutions and the corresponding dual evolution 

dt<^ = G°^<^ on Cb(P 



In order to make this heuristic rigorous at an abstract level, one needs at least some 
tightness on the sequence {f^)N>l for any i, and some convergence 

Gf+i((^)^G,-i(v.) 
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on compact subset of M'^(^+^-'. Both are satisfied for Boltzmann collision processes con- 
sidered in this paper (tightness follows from the moment estimates in Lemma 15.41 for 
instance). 

8.3. Uniqueness of statistical solutions and chaos. We now want, under appropriate 

abstract assumptions, to identify the limit evolution in P{P(M.'^)) or C(,(P(M'^)) obtained 
from the hierarchy, and show that it coincides with the the pushforward evolution semi- 
group Tf° introduced in Subsection 12.31 Meanwhile we shall prove that the statistical 
solutions to the infinite hierarchy are unique, and hence prove the propagation of chaos, 
without any rate, but also under weaker assumptions than previously. 
We shall make the following assumptions 

(Al') Assumptions on the A^-particle system. 

G^ and T/^ are well defined on Ci,{E^) and invariant under permuta- 
tion, and the associated solutions /^ satisfies: 

Wi > 1, the sequence {ff)N>t is tight in Pg^iE)^'^ 

where Gi is a Banach space and Pg^ {E) is defined in Definitions 12.4112. 51) , 
and is associated to a weight function mg^ and a constraint function mg^^, 
and endowed with the metric induced from Qi. 

(A2') Assumptions for the existence of the statistical and 
pushforw^ard semigroups. 

For some 5 £ (0, 1] and some a G (0, oo) we assume that for any a £ 
(a, do): 
(i) The equation (|2.ip generates a semigroup 

5f^ : BPg,,a ^ BPg,,a 

which is 5-Hoder continuous locally uniformly in time, in the sense 
that for any r G (0, cxd) there exists Cr G (0, oo) such that 

yf,g€BPg,,a, sup ||5fV-5f^9L, <Cr \\f-9\\g,- 

te[o,r] 

(ii) The application Q is bounded and (^-Holder continuous from BPg^^a 

into Qi. 
(iii) E \s a, locally compact Polish space and there is J-"i in duality with Qi 

such that J^i is dense in Cb{E) in the sense of uniform convergence on 

any compact set. 

(A3') Convergence of the generators. 

For any fixed £ G N* and any ip G Ch{E^), the sequence 

G;^i(v9) G Cb (^^+^) satisfies G^^^v? ^^^^ G^+i^ 

uniformly on compact sets, where Gf^^ip satisfies the following compatibility 
binary derivation structure: 

for any (^ = c/?i (g) ... (g) (pe e Cb{E)®^ and any V = {vi, ...,vi+i) G i?^+^ 

(8.4) Gt,^{p,){V) = J2\Wvi{^i)\ Q*{^i)iv^,ve+,) 

where Q* is defined from Q through the duality relation 
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Remark 8.1. The identity ()8.4p is called compatibility binary derivation structure for the 
following reasons: compatibility since it is a natural condition in order that any solution ft 
to the nonlinear Boltzmann provides a tensorized solution to the BBGKY hierarchy ()8.3p . 
Indeed, considering such a solution 

ft£C (M+; Pg, (E)) and ip = ipi C^ . . . (g) ipe G Cb{Ef^ 

we compute 

= J2\^{ft,vM {ft® fuQ* {v^)) = {ff^\Gn,^) . 

The word binary refers to the fact that Gf]-^ decomposes in function acting on one variable 
and adding one variable, which corresponds to the binary nature of the collisions. Finally 
the word derivation refers to the fact that the following distributivity property holds 

Let us mention that this distributivity property is at the basis of the original combinatorial 
proof of Kac [32] of propagation of chaos for the simplified Boltzmann-Kac equation. 

(A4') Differential stability of the limiting semigroup. 

We consider some Banach space Q2 D Qi (where Qi was defined in (A2)) 
and the corresponding probability space Pg^ (E) (see Definitions 12. 4112. 5p 
with the weight function mg,^ and the constraint function mg2, and en- 
dowed with the metric induced from Q2- 

We assume that the flow S^^ is UC\^{Pg-^^r,Pg2) foi^ any r G Rg^ in the 
sense that there exists C|? > such that 



IZ I' {^'^''^oi^i-'.^-^P,,, + [S."^cJ.,P,.,„P.,,) "' i Cf. 



Thanks to (A2'), we know from Lemma 12.131 that for any <1> G UCb{Pg^-,^) we may 
define the Co-semigroup T^[^ G UGb{Pg^) by 

and so that <I>j = Tf° [$] satisfies the equation 

with a generator G°° which is a closed operator on UCb(Pg^) and has domain Dom(G°°) 
which contains C^{Pg-^^), and is defined by 

The evolution corresponds to the following dual evolution equation 

(8.5) |(7r„cI>) = (vr„G-[cl>]). 

Our goal is to prove (1) that the evolution equations (18. Sp and (18. 5p (or in other words 
that G°° = G°°, and (2) most importantly that the solution to these equation is given by 
the characteristics equation 

(8.6) y<^€ UCb{Pg,;R) {Trt,<^) = (^0,^^^) • 
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(or in other words that the heuristic generator (5°° introduced for the hierarchy is well- 
defined and equal to G°°). 

Let us explain why the relation ()8.6p indeed defines uniquely a probability evolution 
7ft G P{Pgi)- For any £ G N* we define 



That is a positive linear form on J^®^. Thanks to (A2')-(iii), the Stone-Weierstrass 
density theorem and the Markov-Riesz representation theorem, we conclude that vrf is 
well defined as a element of Pg^{E)®^ . Since now the sequence (vrf) is symmetric and 
compatible, the Hewitt-Savage representation theorem implies that there exists a unique 
probability measure vfi G P{Pgi ) such that for any ip G Tf 

(8.7) (^t,0:=(7ro,rr4). 

Theorem 8.2. Under the asumptions (A1')-(A2')-(A3')-(A4'), for any initial datum 
ttq G P{Pg-^), the flow Ttt defined from (|8.7|) is the unique solution in C{[0,oo);P{Pg-^)) to 
the infinite hierarchy evolution (j8.3|) starting from ttq. 

Moreover, if ttq is fo-chaotic (that is if no = 6fg with /o G P{E)), then nt is S^^^fo- 
chaotic for any t > 0. As a consequence we deduce that if /q^ is fo-chaotic, then f^ is 
Si fo-chaotic. More generally, if f^ converges to ttq then f^ converge to 7tt the associated 
statistical solution. 

Proof of Theorem \8.B . We shall proceed in several steps. 

Step 1: Propagation of chaos. Let us recall that Hewitt-Savage's theorem [38] implies that 
for any vr G P{P{E)) there exists a unique sequence (vr^) G P{E^) such that 

V(^G(a(i?)r, (^^^) = (^,i^^ 

As a consequence, if ttq is /o-chaotic and tt satisfies (18. 7p . then 

nOO of \ rjiQO 



{7tt,e,^) = [7tt,%) = [7To,Tr%)=Tr[R^{fo) 

= K (S^'fo) = {Sr /o, ^i> . . . {Sr /o, ^e) , 

which means that vfj^^ = /® , or equivalently itt = Sf^, and the solution ttj defined by ([87 
is /t-chaotic. 

Step 2: Equivalence between (18. Sp and (18. 5p . 

First let us assume ()8.5p and prove (18. 3p . Consider / G Pg^iE) and <p G J"®^. Then we 
have R^ G C^''^{Pg^{E)) and we deduce from ([8^ that 

(r'+\G?^i'/') = {QUJ),DRUf)) = G^ [K] (/) 

which means 

R, 



-ft /-too ,^ — \J 



Then, using Hewitt-Savage's Theorem again, (j8.5p implies that 

(8.8) I (vr,,t, ^) = I (tt*, i?^) = (vr,, G^ [i?^] ) = (t^, R'+I^^^ ) = (vr,+i,t, G,- ^ M) 
which means that nt satisfies (18. 5p 



f??^iM 



Assume conversely that vr, satisfies (|8.3p and let us prove (I8.5p . One needs to prove that 
one can recover any <1> G UC^{Pgj^{E)) from the previous equation (|87 
Therefore consider $ G UC^{Pg^{E)) and let us define 

(^ = f^^<l>) (y) = <!> te) , y = (t;i,...,7;, 
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and let us write ()8.8p for this choice of ip: 



dt(^*'^'^* 



7rt,G^ 



R' 



■TV R^+1 



Then, on the one hand, for any / G Pg^[E) 

by the law of large numbers. 

On the other hand, for any / G Pg-^{E), we have 

[DR.^{f),Q{fjy 



Y, <^(/y) dQ{f,f)ivi)lldfivj 



For any given i = 1, ...,£, we define 



d-l 



L>$ 



y^v,) ' ^i '■= (^i>--- ,Vi_i,fj+i,...,f^) 



and we write 
Observing that 

and that (<5(/, /), 1) =0 from assumption (A2')-(ii), we find 






'df{v,) 






^— !>00 



{D^{f),Q{fJ)) 



by the law of large number again. This implies (j8.5p . We then conclude that (jS.Sp holds 
for any $ G C^(Pgj(£J)) by density of UC^{Pg^{E)) in this space and the fact that the 
domain of G°° contains C^{Pg^{E)). 



Step 3: Uniqueness. Let us prove that any solution of (I8.3p - ()8.5p satisfies the characteristics 
equation ()8.7p . or in other words that TVt = Ttf. This shall imply uniqueness since we have 
already seen that the vft satisfying (18. 7p is unique. 

The fundamental point here is that for any <l> G UC^{Pg^{E)) if we define <^t '■= T^^, 
thanks to Lemma 12.131 we have 



Then since 



Vi>0, <^t€UC^{Pg,{E))cBom{G° 
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is C^ from the fact that ^t-r G C!^{Pg^{E)) belongs to the domain of G°° for any r, we 
compute 

dT 
and we deduce that 

(7rj,$o> = (7ro,$i) 

which proves that tt^ = vfj satisfies (j8.7p . and concludes the proof. D 



8.4. A remark on stationary statistical solutions. As we have seen: 

• The chaoticity of a sequence of symmetric A^-particle distributions / G P{E ), 
A^ > 1 is equivalent to the fact that the associated vr G P{P{E)) is a Dirac at some 
/o G P{E): IT = 6fg. Hence, in view of Hewitt-Savage's theorem, non-chaoticity 
can be reframed as saying that tt is a superposition of several, instead of one, 
chaotic states. 

• We have recalled the result in [321 US] stating that a chaotic (tensorized) sequence is 
asymptotically concentrated on the energy sphere, which is an effect of the Central 
Limit Theorem. 

• Finally let us make the simple observation that the A^-particle dynamics leaves 
the energy spheres invariant and relaxes on each energy spheres to the uniform 
measure. This is a consequence of the energy conservation laws: at the level of the 
particle system, the dynamics is layered according to the value of this conservation 
law. 

One deduces from these considerations that there is room for non- chaotic stationary states 
of the A^-particle system, namely superposition of several stationary states on different 
energy spheres. Let us make this more precise. 

Lemma 8.3. There exists a non- chaotic stationary solutions to the statistical Boltzmann 
equation. In other words, there exists tt G P(P(M'^)) such that tt ^ 6p for some p G P(M'^) 
and Af^-^iiTi+i) = for any £ G N. 

Proof of Lemma \8.3[ . It is clear that any function on the form 

V G M-^^^+i) ^ 7Te+i{V) = (/>(|y|2) 

is a stationary solution for the equation (|8.3p . that is A£^i{tt£^i) = for any i >1. Now 
we define, with d = 1 for the sake of simplicity, the sequence 

V € R^ ^ ttAV) = , ^f,, ,,,, V£>1, 

where the sequence of positive constants C£ is inductively constructed in the following way. 

• First ci is chosen in a unique way so that vri is a probability measure. 

• Then, once ci, . . . ,q are constructed, q+i is constructed so that n^[7r£+i] = tti, 
which means 

VFgM, / ./„,-./„ dt'* ~ 



lv,m (1 + |y|2 + |^;,|2)™+^/2+l/2 (^ + |y|2)-+^/2- 

This is always possible since 



.,eM(l + |y|2 + |^,|2r+^/2+V2 

Q+1 



dv^ 



(l + |y|2)-+^/2+l/2 J^^^^ , ^^,.m+e/2+l/2 



du* 



ce+1 



(1 + \V\^)rn+i/2 J^^^^ ^^ ^ |^^|2)m+^/2+l/2 



df* 
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which concludes the induction. 

We then deduce that the sequence tti, £ > 1, satisfies (18. 3p since every terms only 
depends on the energy, and also satisfies the compatibility condition n^[7r£+i] = tt^. This 
concludes the proof. CH 
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